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FORTRAN PROGRAMS FOR CALCULATING WIND-TUNNEL 
BOUNDARY INTERFERENCE 


By Harry H. Hey son 
Langley Research Center 

SUMMARY 

Boundary-interference programs, developed in NASA TR R-302, are presented 
without comment. These programs should be utilized only after careful consideration of 
the assumptions and procedures given in that report. 

INTRODUCTION 

Reference 1 develops a systematic computer procedure for calculating the wind- 
tunnel interference factors for arbitrary configurations from the interference calcula- 
tions for a vanishingly small model. The method is not limited to any one tunnel con- 
figuration since it is necessary only to substitute a subroutine appropriate to the tunnel 
for that given herein. 

The underlying theory (subroutine DLTAS) in the present computer programs is that 
of reference 2. It is directly applicable to models which produce large wake deflections, 
such as V/STOL models. These programs may also be used directly for more conven- 
tional testing at moderate lift coefficients by means of the few simple modifications 
described in reference 1. 

No sample calculations or check cases are provided herewith. The numerical 
values provided in reference 1 should be adequate for this purpose. The reader is 
cautioned against using these programs without first carefully considering the assump- 
tions, limitations, and procedures given in reference 1. 



COMPUTER PROGRAMS 


The programs are given in the appendixes. The following table should aid in 
locating the program of interest: 


Appendix 

Model 

Interference 

Page 

A 

Small 

At point 

4 

B 

Wing 

Average 

6 

C 

Wing 

Span distribution 

9 

D 

Wing 

At tail 

12 

E 

Jet 

At wing 

16 

F 

Jet 

Wing distribution 

20 

G 

Jet 

At tail 

24 

H 

Rotor 

Average 

28 

I 

Rotor 

Lateral axis 

31 

J 

Rotor 

Longitudinal axis 

34 

K 

Rotor 

At tail 

37 

L 

Tandem rotors 

Average 

41 

M 

Unloaded rotors 

Average 

46 

N 

Unloaded rotors 

At tail 

52 

O 

Side-by-side 

rotors 

Average 

57 

P 

Side-by-side 

rotors 

At tail 

64 

Q 

Subroutine DLTAS 

70 
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appendix a 


FORTRAN PROGRAM FOR CALCULATING WIND-TUNNEL INTERFERENCE 
NEAR A VANISHINGLY SMALL MODEL 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON COC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE OLTAS WHICH IS GIVEN IN 
APPENDIX Q. 


THIS IS THE BASIC WIND-TUNNEL INTERFERENCE PROGRAM FOR WHICH THE DERIVA- 
TIONS ARE GIVEN IN NASA TR R-124 IREF 2). THE SUCEEDING PROGRAMS ARE ALL 
DEVELOPED FROM THIS ONE PROGRAM. INPUT WILL BE FOUND AT ADDRESS 1 (ONE CARD 
PER CASE) IN FORMAT 103. THE REQUIRED INPUT VARIABLES ARE 


ZETA SEMIHEIGHT OF WIND TUNNEL OIVIDED BY HEIGHT OF MODEL ABOVE FLOOR 

ETA DISTANCE FROM MODEL TO RIGHT-HAND WALL DIVIDEO BY WIND-TUNNEL 

SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

XCVERH LONGITUDINAL POSITION OF POINT AT WHICH INTERFERENCE IS NEEDED, 

NONDI MENSIONALIZED WITH RESPECT TO TUNNEL SEMIHEIGHT 

YOVERH LATERAL POSITION OF POINT AT WHICH INTERFERENCE IS NEEDED, 
NONDIMENSIONALI ZED WITH RESPECT TO TUNNEL SEMIHEIGHT 

ZOVERH VERTICAL POSITION OF POINT AT WHICH INTERFERENCE IS NEEOED, 
NONDIMENS IONAL IZED WITH RESPECT TO TUNNEL SEMIHEIGHT 


PROGRAM WINOTUNIINPUT, OUTPUT, T APE5= I NPUT , T APE6 =OUT PUT I (A l) 

COMMON ZETA, ETA, GAMMA, XOVERH, YOVERH, ZOVERH, DELTAI28) (A 2) 

DIMENSION C( 8) (A3) 

DATA (C( I ) , 1 = 1 ,8 )/ 20 . , 30 . , AO . ,50. ,60. ,70. ,80. ,90./ (A A) 

1 READ (5,103) ZETA, ETA, GAMMA, XOVERH, YOVERH, ZOVERH (A 5) 

IF (EOF, 5) 999,2 (A 6) 

2 WRITE (6,148) G AMM A , ZE T A ,E TA , XO VERH ,YOVE RH , ZOV ERH (A 7) 

WRITE (6,210) (A 8) 

WRITE ( 6,211 ) ( A 9) 

WRITE (6,212) ( A 10) 

WRITE (6,213) (A 11) 

WRITE (6,214) (A 12) 

WRITE (6,215) (A 13) 

WRITE (6,216) (A 14) 

WRITE (6,217) (A 15) 
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Appendix A - Concluded 


WRITE 16,218) (A 16) 

DO 41 K=l,8 ( A 17) 

CALL DLTAS (C(K>) (A 18) 

*****44************** SEE APPENDIX Q FOR SUBROUTINE DLTAS *********** *********** 
WRITE 16,149) CIK) ( A 19) 

WRITE 16,150) IDELTAU ) ,1=1 ,25,4) (A 20) 

WRITE 16,151) (DELTA! I ) , 1=2,26,4) (A 21) 

WRITE (6,152) (DELTA! I ) ,1=3,27,4) (A 22) 

WRITE (6,153) (DELTA! I), 1=4,28,4) (A 23) 

41 CONTINUE (A 24) 

GO TO l (A 25) 

103 FORMAT (6F10.3) (A 26) 

148 FORMAT ( 1H1 //35X*I NT ERFEREMCE FACTORS AT A POINT NEAR A VANISHING! (A 27) 

1Y SMALL M00EL*//35X*GAMM A =*F8.3,9X*ZETA =*F8.3 , 11X*ETA =* (A 28) 

2F8.3//35X*X/H =*F8.3,9X*Y/H =*F8. 3 , 11 X*Z /H =* F8 .3 H ) (A 29) 

1.49 FORMAT ( /5X6H CHI =F6.2/) (A 30) 

150 FORMAT ( 3 X5H( W,L )7( F 17 .4 ) ) (A 31) 

151 FORMAT (3X5H(U,L )7( F17.4 ) ) (A 32) 

152 FORMAT (3 X5H ( W ,D )7 ( F17 .4 ) ) (A 33) 

153 FORMAT ( 3X5H( U, D )7( F17.4 )/ ) (A 34) 

210 FORMAT (1X13K1H-)) (A 35) 

211 FORMAT (1X1HI11X IHI 31 X61 HCORRECT ION FACTORS FOR CORRECTING FROM A (A 36) 

1WIND TUNNEL WHICH IS25X1H!) (A 37) 

212 FORMAT ( 1 XIH 11 IX 1H I 1 17 ( 1H-) IHI ) (A 38) 

213 FORMAT ( 1X1H 11 IX 1H 1 16X1H I5X6HC L0SED5X1HI 16X IHI 2X12HCLOSED FL00R2X1 (A 39) 

1HI6X4H0PENSX1HI 16X1HI5X6HCL0SED4X1HI ) (A 40 

214 FORMAT ( 1X1HI3 X5HDEL T A3 X IHI 5X6HCL0S ED5XIHI 4X9H0N B0TT0M3X 1H I 6X4H0P (A 41) 

1EN6X1H I 6X4H0NL Y6X1HI 5X5HFL00R6X1H 1 5XSHCL OS ED5XI H I3X9H0N BOTTOM! XIH (A 42) 
21 ) ( A 43) 

215 FORMAT ( 1 XIH II IX IHI 1 6X1H 1 6X4H0NL Y6X 1H I 16X1 8H I (GR CUND EFFEC T ) I6X4H ( A 44 ) 

10NL Y6X 1H 1 1 6X1H 1 6 X4H0NL Y5 XIH I ) ( A 45 ) 

216 FORMAT ( 1 XI H 1 1 IX 1H I 84 ( 1H- ) l H 13 2( 1H- ) IHI ) (A 46) 

217 FORMAT ( 1 XI H 1 1 1 X 1H 1 36X1 1 HTO FREE AI R 37X1HI 8 XI 6HT C GROUND EFFECT 8X1 (A 47) 

IHI) (A 48) 

218 FORMAT ( 1X131 ( 1H-) / ) (A 49) 

999 STOP (A 50) 

END (A 51) 
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APPFNDIX B 


FORTRAN PROGRAM FOR CALCULATING THE AVERAGE WIND-TUNNEL 
INTERFERENCE OVER A SWEPT WING 


THIS PROGRAM WAS WRITTEN IN COC FORTRAN, VERSION 2.1, TO RUN ON CDC 5100 
SERIES COMPUTERS WITH THE SCOPF 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
APPENDIX g. 


INPUT WILL BE FOUND AT ADDRESS 1 (ONE CARD PER CASF) IN FORMAT 903. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE APEX OF THE SWEPT LIFTING LINE. 
THE REQUIRED INPUT VARIABLES ARE 


LI LOAD INDICATOR, LI=1 FOR JNIFORM LOADING, L=2 FOR ELLIPTIC 

LOADING 

ZETA 1 SEMI HP I GH T OF TUNNEL DIVIDED BY HEIGHT OF ORIGIN ABOVE FLOOR 

ETA I UISTANCF FROM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEM IW I DTK 

GAMM4 WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMA RATIO OF WING SPAN TO TUNNEL WIDTH 

LAMBDA WING SWEEP ANGLE, DEG 

AL°HA 4NGL c OF ATTACK OF WING, DEG 


PROGRAM W INDTUN( INPUT, OUTPUT ,TAPE5 = INPUT , T A PE 6 = 0U T PUT > IB 1) 

COMMON ZET A, FT A, GAMMA , XOVFRH , YOVERH , ZO VERH , DELTA ( 28 ) <R 2) 

DIMENSION XDFLTAI28) ,XLOAOC 10), XLE( 10 J,G<8) (B 3) 

REAL LAMBDA (BA) 

CATA (C( I ) , 1 = 1 , 8 >/2o ., 30 . ,t*r . , 5 r . ,50. , 70 . , pr. , 90 ./ (P 5 ) 

XLE(1) = XLE ( 43579 (R 6 ) 

XL E(2 ) = XLE (9 )=p.7142? (P 7) 

XLEI 3)=XLE< 8 ) =''.356"'3 (R P) 

XLE(4 ) = XL E ( 7 ) =^.05394 (B °> 

XLE(5)=XLE( 6)=0.9°499 (B l^t 

DO 8 03 L 1= 1, 2B (BID 

B-p 3 XDELTAl LI )=p. (B 12) 

PI=3.1B1592f5P5 n 979 ( B 13 ) 

RAO= .3174532 9251 99 (B 1 A- ) 

1 READ (5,903) LI, ZETAl.ETAl, GAMMA, GIGMA, LAMBDA, ALPHA (B 15) 

IF ( EOF, 5 ) 909,47 (6 16) 
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Appendix B - Continued 


47 IF ( L I • E Q . 1 ) GO TO 804 
IAL?HA=8HEtLl PTIC 
SUML=0. 01261^4 

00 808 M2 = I * l 0 

80 8 X L 04 D ( M 2 ) = X L E ( M ? ) 

GO TO 160 
80 4 SUML =0 • 0 1 

1 ALP HA=8HUN I FORM 
DO 809 M 2 = 1 , 1 9 

809 XL OAD( M2 )=1 .0 

160 WRITE (6,901) I AL PH A , G AMM A , ET 4 1 , S IG MA , Z E T A 1 , AL PH 4 , L AM BDA 
WRITE ( 6,210 ) 

WRITE l 6 , 2 L 1 ) 

WR ITE ( 6,212 ) 

WRITE (6,218 ) 

WR ITE (6,214 ) 

WRIT E (6,216 ) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218 ) 

C ONS T 1= 1. 

L A MB DA=L AM BD A *R AD 
ALPHA=4LPHA*RAC 
DO 41 K = 1 , 8 

IF ( S IGM4 .NF.O . ) GO TO 811 
M6 =M 7 = N6 =N 7 = 1 
X L 0 A Q ( 1 ) = 1 .0 
SUML = 1 . 

GO TO 812 

811 IF (ET41.NE.1.) GO TO 813 
N6 = i 

M6 =6 

N7=M 7= 1^ 

CONST 1=2 . 

GO TO 812 
813 M6=SI6 = 1 
M 7 = N 7 = 1 0 

812 DO 301 M 1 = M 6 , M 7 
00 802 N 1 = N6 , N 7 

XSTAR=( 11 . -2. AFLOAT (Ml) i / in. 

YSTAR=(?.*FL04T(N1 )-ll*)/in. 

ZSTAR=(11*-2# jSc FLQAT ( N 1 ) )/io. 

ET A=ETAl + YSTAR*SIGMA 


(B 17) 
(B 18) 
< B 19) 
( 9 20 ) 
( B 21 ) 
( B 22) 
( 8 2 3 1 
( 8 24) 
(B 26) 
(B 26) 
(B 27) 
(B 28) 
(B 29) 
(B 30) 
(8 31) 
( B 32 ) 
( B 33 ) 
(B 34) 

( B 35 ) 
(8 36) 
( B 37) 
(B 3 B ) 
(B 39) 
(R 40) 
( B 41) 
(B 42) 
(B 43) 
( B 44) 
(B 45) 
(B 46) 
(9 47 ) 
(B 48) 
(8 49) 
IB 5<M 
(B 51) 
( B 52 ) 
( 8 5 3) 
( R 5 4) 
(ft 55) 
(B 56) 
( B 57) 
( B 58 ) 
(8 59 ) 


ZETA=ZETA1/ (1 •-ABS(YSTAR)*SIGM4*G4MM4*ZETA1*T4N(L4M0DA ) *S IN ( ALPHA) ( B 60) 


1) ( B 61 ) 

XOVERH=S IGMA*GAMMA *TAN (L AMBDA) *COS ( AL P HA ) * ( AB S ( X ST AR ) - A B S ( Z S T AR ) ) ( B 62 ) 

Y9VERH=( FLOAT ( Ml ) -FLOAT ( N1 ) ) *S TGMA*G4MMA*( - .2 ) ( B 63) 

ZOVERH= SIGMA *GAMMA*T AN(LAM3DA)*$ TN(ALPHA)*(A8S(ZSTAR)~ABS(X$TAR > ) (8 64) 

CALL DLTAS (C(K) ) ( B 65) 

fc* ** ****** ****** ***** SEE APPENDIX Q FOR SUBROUTINE DLTAS ***** *** **** ******** ** 

DO 80 5 L 1 = 1, 2 8 ( B 66) 

805 XDELTA(L1 )=XJELTA(L1 )+DELTA(Ll)*XLOAD(Nl) (B 67) 

302 CONTINUE (B 68) 

30 1 CONT INUE ( B 69) 

DO 807 L 3=1, 28 ( B 70 ) 

807 DELTA(L3)=XDELTA<L3)*SUML*C0NST1 (B 71) 

WRITE (6,149) C(K) (872). 

WRITE (6,150) (DELTA(I), 1=1,25,4) (B 73) 

WRITE (6,151) (DELTA( I) , 1 =2 ,2 6,4 ) ( 8 74) 

WRITE (6,15?) (DELTA! I ),I=3,?7 ,4) ( 8 75) 

WRITE (6,153) (DELTA( I ), 1 = 4,28,4) (B76) 1 
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Appendix B - Concluded 


DO 8H 1.4=1,29 1 8 7 7 ) 

814 XDELT 4IL4 ) =0 • IB 78) 

41 CONTINUE (B 79) 

GO TO 1 (B 80) 

149 FORMAT (1X*CHI =*F7 .3/ ) |B 81) 

150 FORMAT I 3 X5H ( W , L ) 7 { F 1 7 . 4 ) ) ( 8 82) 

151 FORMAT ( 3X5H(U,L )7(F1 7.4 ) ) (B 83) 

152 FORMAT ( 3X5H ( W , 0 ) 7( F 1 7. 4 ) ) (8 84) 

153 FORMAT I 3X5HIU ,0 )7 ( FI 7.4) / /) IB 85) 

210 FORMAT 11X131 I IF- ) ) IB 86) 

211 FORMAT I1X1HI11X1H 131X61 (-CORRECTION FACTORS FOR CORRECTING FROM 4 IB 87) 

1WIND TUNNEL WHICH IS25X1HI) IB 88) 

212 FORMAT I IX 1HI 11X1HI11 71 1H-)IHI ) IB 89) 

213 FORMAT 1 1 X 1H II 1 X1HI16 X 1HI 5X6HCL0SE0 5X 1H I 1 6X1H I 2 XI 2HCL0 SEO FLOOR2X1 IB 90) 

1 H I 6X4H0PEN6 X1H II 6X1HI 5X6HCLQSED4X 1HI ) IB 91) 

214 FORMAT I 1 X 1 H 1 3 X5HOEL T A3X 1 HI 5X 6HCL0 SE05 XIH I 4X9H0N BOT TOM3X 1H I 6X4 HOP IB 92) 

1 EN6X1HI6X4HQNL Y6X1 HI 5X5HFLOQR6 X1HI5X6HCL0SE05XHI3 X9H0N BOTTOM3 X 1H !R 93) 
21 ) (B 94) 

215 FORMAT I IX l H 1 1 1 XI H 1 16X1HI 6X4H0NLY6X 1H 1 1 6X1 8H I I GROUND EFFECT) I6X4H IB 95) 

10NLY6XIHI16X1HI6X4H0NL Y5X1HI ) IB 96) 

216 FORMAT I IX 1HI11X1H! 841 1H-) 1HI32I 1H- ) 1H I ) IB 97) 

217 FORMAT I 1 XI H 11 1 XI H 1 36 XI 1HTO FREE A I R 37X 1H I 8X1 6HT0 GROUND EFFECT 8X1 IB 98) 

1 H I ) ( B 9 9 ) 

216 FORMAT 11X131 I 1 H- ) / ) IB 100) 

900 FORMAT I II ,F9.3,5F10.3) (B 101) 

901 FORMAT. I 1H1 / / t 42 X *AV ER AGE INTERFERENCE OF SWEPT WING OF FINITE SPA IB 102) 

IN* // 58XA 8 , * L0ADING*//36X*G4MM4 =*F6. 3, 10X*£ TA =*F7.3,19X 18 103) 

2*S IGMA =*F7 . 3//36X*ZET4 =*F6. 3 , 10X#ALPH4 =*F7.3,10X IB 104) 

3*L AM9DA =*F7.3//J IB 10 5) 

999 STOP IB 106) 

END IB 107) 
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APPEND! X C 


FORTRAN PROGRAM FOR CALCULATING THE DISTRIBUTION OF 
WIND-TUNNEL INTERFERENCE OVER THE 
SPAN OF A SWEPT WING 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC SDO" 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
APPENDIX Q. 


INPUT WILL BE FOUND AT ADDRESS 1 (ONE CARD PER CASE) IN FORMAT 103. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE APEX OF THE SWEPT LIFTING LINE. 
THE REQUIRED INPUT VARIABLES ARE 


LI LOAD INDICATOR, LI = 1 FOR UNIFORM LOADING, L=2 FOR ELLIPTIC 

LOADING 

ZETA1 SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMA RATIO OF WING SPAN TO TUNNEL WIDTH 

LAMBDA WING SWEEP ANGLE, DEG 

ALPHA ANGLE OF ATTACK OF WING, DEG 

C EFFECTIVE WAKE SKEW-ANGLE, DEG 


IN SYMMETRICAL CASES THIS PROGRAM COMPUTES THE INTERFERENCE DISTRIBUTION 
OVER ONE SEMISPAN ONLY. THIS PROGRAM REJECTS CASES OF ZERO SPAN. FDR SUCH 
CASES, THE INTERFERENCE IS UNIFORM AND THE VALUES ARE IDENTICAL TO THOSE PRO- 
VIDED BY THE PROGRAM OF APPENDIX B. 


PROGRAM WINDTUN< INPUT, OUTPUT, T APE 5= INPUT , TA PE6=0UTPUT) (C 1) 

COMMON ZETA, ETA, GAMMA, XO VE RH ,Y0 VERH ,ZOVERH, DELT A (2 8 ) (C 2) 

DIMENSION XDELTAI28),XL0AD( 10),XLE( 10) IC 3) 

REAL LAMBDA <C 4) 

XLE(1)=XLE(10)=. 43579 IC 5) 

XLE(2)=XLE(9)=.71422 (C 6) 
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Appendix C - Continued 



XLE(3) = XLE(8)=. 86603 


(C 7) 


XLE(4)=XLE(7)=. 95394 


(C 8) 


XLE( 5)=XLE<6)=. 99499 


(C 9) 


CO 805 N 2 = 1 , 2 8 


( C 

10 ) 

805 

XDELT A ( N2 I = 0 . 


(C 

11 ) 


RAO=. 0174532925199 


(C 

12 ) 

I 

READ (5,103) LI, Z ET Al , E TA1 , G AMMA ,S I GM4 , LAMBDA, ALPHA, C 


( c 

13) 


IF (EOF, 5) 999,814 


(C 

14) 

8 14 

IF (LI.EC.l ) GO TO 806 


(C 

15) 


IALPHA=8FELIIPTIC 


(C 

16) 


SUMl= • 1 26 104 


(C 

17) 


DO 808 M2 = l , 10 


(C 

18) 

808 

XL0AC(M2»=XLE(M2 1 


(C 

19 ) 


GO TO 803 


(C 

20) 

8C6 

SUML= • 1 


(C 

21 ) 


I ALPHA = 8HUN IFORM 


(C 

22 ) 


DO 809 M2=l, 10 


(C 

23) 

809 

XLOAO( M2 1 = 1 .0 


(C 

24) 

803 

XOVERH=ZO VERH=0. 


(C 

25) 


H1=0 


(C 

26) 


IF (ETA1.E0.1. ) Ml =6 


(C 

27) 


WRITE (6,930) GAMM A, ZET A 1 , I ALPHA , ALPHA , S IGM A , ET Al, L AM8DA ,C 


(C 

28) 


WRITE <6, 210) 


( c 

29) 


WRITE (6,21U 


(C 

30 ) 


WRITE (6,212) 


( c 

31 ) 


WRITE (6,213 ) 


( c 

32) 


WRITE (6,214) 


(C 

33) 


WRITE (6,215) 


( c 

34 ) 


WRITE (6,216) 


( c 

35) 


WRITE (6,217) 


(C 

36) 


WRITE (6,218) 


(C 

37) 


IF (SIGMA. NE.O.) GO TO 813 


(C 

38) 


WRITE (6,901 ) 


(C 

39) 


GO TO 1 


(C 

40) 

813 

ALPHA* ALPhA*RAD 


(C 

41) 


L AMBD A=L AMBDA*R AD 


(C 

42) 

804 

X ST AR= ( 1 1 .-2 .*FLOAT ( Ml) )/13. 


(C 

43 ) 


CO 800 N 1=1 , 10 


(C 

44) 


YSTAR=(2.*FL0AT(Ni)-ll. 1/19. 


(C 

45) 


ZST AR= ( 11.- -2. AFLOAT (Nil )/10. 


(C 

46) 


ETA=ETAl*YSTAR*SIGMA 


(C 

47) 


ZETA=ZETAl/( l .-ABS ( YST AR ) *S IGMA*GAMMA*ZE TA1 *T AN I LAMBDA )*S IN ( ALPHA ) 

(C 

48) 


11 


(C 

49) 


XOV£RH= 5 I GMA *GAMMA *T AN( L AMBD A ) *COS ( AL PH A ) * ( ASS ( X STAR ) -AB S ( 1 ST AR 

) ) 

(C 

50) 


YCVERH=(FL0AT(M1)-FL0AT (Nl) ) *S IGMA*G AMM A* ( - . 2 ) 


(C 

51) 


ZOVERH=S IGMA*G4MMA*TAN( LAM3DA) *SIN( ALPHA )*( ABS ( Z ST AR ) -ABS ( XST4R 

) ) 

( c 

52 ) 


CALL DLTAS (C ) 


(C 

53) 

4 * 4 ***************** SEE APPENDIX Q FOR SUBROUTINE DLTAS 4444 **** 44 * 

*********** 


DO 801 N2=l , 28 


(C 

54) 

801 

XOEL TA ( N2 ) =X0ELT4( N2 I*0ELTA( N2 I*XL0A0(N1 ) 


(C 

55) 

800 

CONTINUE 


(C 

56) 


DO 802 N2=l,28 


(C 

57) 

802 

CELT A ( N2 )= XDELT A( N2 J4SUML 


(C 

58) 


WRITE (6,149) XSTAR 


(C 

59) 


WRITE (6,150) (DELTA( I ) , 1 = 1, 25,4) 


(C 

60) 


WRITE (6,151) (DELTA( I ), 1=2,26,4) 


(C 

61 ) 


WRITE (6,152) (DELTA! I 1, 1=3, 27,4) 


(C 

62) 


WRITE (6,153) (DELTA! I ), 1 = 4, 28,4) 


(C 

63) 


CO 810 N2 = l , 2 8 


(C 

64) 

810 

XDELTA( N2 ) =0 .0 


(C 

65) 


M1=M1 + 1 


(C 

66) 


10 



Appendix C - Concluded 


IF (M1.LTU2 I GO TO 804 (C 

GO TO 1 (C 

103 FORMAT ( 1 1 , F9.3, 6F 10 . 3) (C 

149 FORMAT ( 10X12HY/ SEMI SPAN =F4.1/> (C 

150 FORMAT ( 3X5H ( W ,U 7( FI 7.4 ) ) (C 

151 FORMAT (3X5H (U • L )7 (F17.4) I (C 

152 FORMAT (3X5H(W,D I7CF17.4 )> (C 

153 FORMAT ( 3X5H(U ,D >T( F 17. 4 >// ) (C 

210 FORMAT C 1 XI 31 t 1 H- ) > (C 

211 FORMAT ( 1X1HI11X1HI31X61HC0RRECTI0N FACTORS FOR CORRECTING FROM A IC 

1WIND TUNNEL WHICH IS25X1HI) (C 

212 FORMAT ( 1 XIH 11 1 XIH I 117( 1H- ) 1HI) (C 

213 FORMAT I IX 1H II 1 X1HI 1 6X1HI5X6HCL0 SED5X 1HI 1 6X 1HI 2X 12HCL0SFD FL00R2X1 (C 

1HI6X4H0PEN6X1HI 16X1HI 5X6HCL 0SED4X1 H I) (C 

214 FORMAT (1XIHI3X5F0ELT A3X1HI 5X6HCLOSED5X1HI 4X9H0N B0TT0M3X1HI6X4H0P (C 

1EN6X1HI 6X4H0NL Y6X1HI 5X5HFL00R6X1HI 5X6HCL0S ED5XI H I3X9H0N BDTT0M3 XIH ( C 
2 I ) (C 

215 FORMAT (1X1HI11X1HI 16X1HI6X4H0NLY 6X1HI16X18HI (GRCUND EFFECT) I6X4H (C 

10NLY6X1HI16X1HI6X4H0NLY5X1HI ) ( C 

216 FORMAT I 1X1H 11 1 XIH I 84 1 1H-) 1H 13 2( 1H- ) 1HI ) (C 

217 FORMAT ( 1 XIH II IX IH I 36X 1 1HT0 FREE A IR 37XIHI 8X1 6HT0 GROUND EFFECT 8X1 (C 

1HI) (C 

218 FORMAT (1X131* 1H— >/l <C 

900 FORMAT ( 1 HI / //38X*I NT ERFERENCE DISTRIBUTION OVER SWEPT WING OF FIN (C 

1ITE SPAN*///15X*GAMMA =*F7.3*15X*ZET4 =* F7 . 3 , 1 5X , A8* LOADING* (C 

21 5X*ALPHA =*F7.3//15X*SIGMA = *F7 . 3» 15X*ET A =*F7 .3, 15X*L AMBDA = * (C 
3F8.3. 15X*CHI =*F7.3//) (C 

901 FORMAT (/ ///40X*S IGMA ECUALS ZERO - JSE AVERAGE INTERFERENCE PRDGR (C 

1AM*///) (C 

999 STOP (C 

END (C 


67) 

68 ) 

69) 

70) 

71) 

72) 

73) 

74) 

75) 

76) 

77) 
78 ) 

79) 

80) 
81) 
82) 
83) 
84 ) 

85) 

86 ) 

87) 

88 ) 
89 ) 

90 ) 

91) 

92) 

93) 

94) 
95 ) 

96) 

97) 
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APPENDIX D 


FORTRAN PROGRAM FOR CALCULATING THE AVERAGE WIND-TUNNEL 
INTERFERENCE OVER A TAIL 
BEHIND A SWEPT WING 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CQC 6090 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH GARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
APPENDIX Q. 


INPUT WILL BE FOUND AT ACCRESS 1 (TWO CARDS PER CASE! IN FORMAT 900. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE APEX CF THE SWEPT LIFTING LINE. 
THE REQUIRED INPUT VARIABLES FOR THE WING, GIVEN ON THE FIRST CARD, ARE 


LOAD INDICATOR, LI*l FOR UNIFORM LOADING, L=2 FOR ELLIPTIC 
LOADING 

SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT OF ORIGIN ABOVE FLOOR 

DISTANCE FRCM ORIGIN TO RIGHT-HANO WALL DIVIDED BY TUNNEL 
SEMIWIDTH 

WIDTH-HEIGHT RATIO OF WIND TUNNEL 
RATIO OF WING SPAN TO TUNNEL WIDTH 
WING SWEEP ANGLE, DEG 
ANGLE OF ATTACK OF WING, DEG 


THE REQUIRED INPUT VARIABLES FOR THE TAIL, GIVEN ON THE SECOND CARD, ARE 


SIGMAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL TAIL LENGTH BEHIND ORIGIN AT ZERO ANGLE OF ATTACK, NOND IMENS ION- 

AL1ZED WITH RESPECT TO TUNNEL SEMIHEIGHT 

TL TAIL LENGTH BEHINO ORIGIN AT ZERO ANGLE OF ATTACK, NINO I ME NS I ON- 

AL I ZED WITH RESPECT TO TUNNEL SEMIHEIGHT 


(D 1) 
(D 2) 


12 


PROGRAM WINDTUNt INPUT, OUTPUT *TAPE5 3 INPUT , TAPE 6= OUT PUT) 
CCMMON ZETA,ET A »G AMM A » XO VER H , YOVE RH, ZOVE RH , DEL T A ( 28 » 


LI 

ZETA1 

ETA1 

GAMMA 

SIGMAW 

LAMBDA 

ALPHA 



Appendix D — Continued 


DIMENSION XDELTAI28) ,XLOAO( 10 I ,XLE( 10) ,C(8) ID 3) 

REAL LAMBDA I D 4) 

DATA (C(I), 1=1, 81/20. ,30. ,40. ,50. ,60. ,70. ,80. ,90./ (0 5) 

XLE( 1)=XLE(10)=C .43579 (D 6) 

XLE(2)=XLE(9)=0.7IA22 (D 7) 

XLE(3 )=XLE(8 )=0 .86603 (D 8) 

XLE(4)=XLE(7)=0. 95394 (D 9) 

XLE(5)=XLE(6)=0. 99499 (0 10) 

DO 803 Ll=l ,28 (D 11) 

803 XDELTAILl ) = 0. (D 12) 

BA D=. 0174532 92 5159 (D 13) 

Pl=3. 14155265358979 (D 14) 

1 BEAD (5,900) L I , ZE T Al , ETA1 .GAMMA ,S IGMAW .LAMBDA , AL PH A, (D 15) 

1 SIGMAT ,TL*TH (D 16) 

IF (EOF, 5 ) 999, A7 (D 17) 

A 7 IF ILI.EQ.l) GO TO 80A (D 18) 

I ALPHA=8 FELL IPT IC (D 19) 

SUML=0. 031526 (D 20) 

DO 808 M2=l , 10 (D 21) 

808 XL0AD(M2)=XLE(M2 ) (D 22) 

GC TO A 8 (D 23) 

80A SUML=0 .025 ID 2A) 

I ALPHA=8HUNIF0RM (D 25) 

DO 809 M2 = 1 , 10 (D 26) 

809 XLDAD(M2 )=1 .0 (D 27) 

A8 WRITE (6,901) S IGMAW ,TH, IALPHA .GAMMA , ZET A1 , S IGMA T, TL, LAMBDA , ALPHA , (D 28) 

1ETA1 (D 29) 

WRITE (6,210) (D 30) 

WRITE (6,211) (D 31) 

WRITE (6,212) (D 32) 

WRITE (6,213) (D 33) 

WRITE ( 6 , 21 A ) (D 3A ) 

WRITE (6,215) (D 35) 

WRITE (6,216) (D 36) 

WRITE (6,217) (D 37) 

WRITE (6,218) (D 38) 

L AMBDA=L AM8DA*RAD (D 39) 

ALPHA* ALP FA*R AD (D AO) 

DO A 1 K= 1 , 8 ( D Al ) 

IF (SIGMAW. EO.O..ANO. SIGMAT. EQ.O.) GO TO 850 (D A2) 

IF (SIGMAW. EQ.O.. AND. SIGMAT. NE.O.) GO TO 855 (D A3) 

IF (SIGMAW. NE.O. .AND. SIGMAT. EQ.O.) GO TO 860 (0 AA) 

M 7=4 ( D A5 ) 

N7= 10 (D A6) 

CCNST1=1 .0 ID A7> 

IF (ETA1.NE.1.) GO TO 812 <D A8 ) 

M7=2 (D 49) 

C0NST1=2.0 (D 50) 

GO TO 812 (D 51 ) 

850 M7=N7= 1 (D 52) 

XLO AC (1 )= 1.0 (D 53) 

SUML=0.025 (D 54) 

CONST 1=40 .0 (D 55) 

GO TO 812 (D 56) 

855 M7=4 ID 57) 

N7=l (D 58) 

XLO ADI 1 ) = 1 .0 (D 59) 

SUML=0.025 (D 60) 

C0NST1=10 .0 (D 61) 

IF (ETA1.NE.1.) GO TO 812 ID 62) 

M7=2 ( D 63 ) 


Appendix D - Continued 



CONST 1 = 20 .0 


( 0 

64) 


GO TO 812 


( 0 

65 ) 

860 

M7= 1 


(D 

66) 


N7=10 


( 0 

67) 


C ONST 1 = 4.0 


( 0 

68 ) 


IF (ETAl.NE.l.) GO TO 812 


iO 

69) 


N7=5 


( D 

70) 


CCNST1=8 .0 


( 0 

71) 

812 

DC 801 M 1= 1 » M7 


{ 0 

72 ) 


DO 802 K1=1,N7 


( 0 

73) 


XSTAR=l 2.* FLO AT ( M )-ll . ) /iO. 


( 0 

74) 


YSTAR=C ll.-2.*FL0AT(Nl)) /IO. 


(D 

75) 


ZSTAR = I 5.-2.*FL0AT(Ml))/4. 


( D 

76) 


ETA=ETA1*XSTAR*S IGMAW 


( D 

77) 


ZETA=ZETAl/< l.-ABS(YSTAR )*S IGM AW*GAMMA*ZETA1 *T AN (L AMBDA)*SIN( AL PH A 

(0 

78) 


1) ) 


(0 

79) 


XCVERH=TL*COS( ALPHA) *TH*S INI ALPHA )-S IGM AW*GAMMA* TAN (L AMBDAI*COS< AL 

( 0 

80) 


1PHA)*ABS(YSTAR> 


( 0 

81) 


YOVERH= ZST AR*S IGMAT*GAMM A-Y ST AR*S I GMAW*GAMMA 


(0 

82) 


ZOVERH= TH*COS I AL PH A ) -TL* S I N ( AL PH A ) +S IGM AW*G AMMA*T AN( L AMBO A ) *S IN ( AL 

( 0 

83) 


IPHA)*ABS( YSTAR > 


(D 

84) 


CALL OLTAS (C(KII 


( D 

85) 

*4* ************ ****** SEE APPENDIX 0 FOR SUBROUTINE DLTAS ********************** 


DO 805 Ll=l ,28 


(0 

86) 

805 

X DELTA l LI )=XDELTA( LI ) + ( OELT A ( L 1 ) *XLO AD < Nl) I 


( 0 

87) 

802 

CONTINUE 


( D 

88) 

801 

CONTINUE 


( D 

89) 


DO 807 L3 = l » 2 8 


(0 

90) 

807 

CELTA(L3)=XDELTA(L3i*SUML*CONSTl 


( D 

91 > 


WRITE (6,149) C(K) 


(0 

92) 


WRITE (6,150 (OELTA(I) ,1=1,25,4) 


(D 

93 ) 


WRITE (6,151 1 (DELTA(I>, 1 = 2, 26,41 


( D 

94) 


WRITE (6,152) (DELTA(I), 1=3,27,4) 


(D 

95 ) 


WRITE (6,153) (DELTA! I ), 1=4,28,4) 


(0 

96) 


DO 814 L4=l ,28 


( 0 

97) 

814 

XDELTA(L4)=0. 


< 0 

98 > 

41 

CONTINUE 


(0 

99) 


GO TO 1 


(D 

ICO) 

149 

FORMAT (5X5HCHI =F7.3/> 


(0 

1C1) 

150 

FORMAT (3X5H(W,L)7(F17.4)) 


( D 

102) 

151 

FORMAT (3X5H(U,L )7(F17.4)) 


(0 

103) 

152 

FORMAT (3X5H(W,D)7(F17.4)) 


(0 

1C4 ) 

153 

FORMAT (3X5H(U,D)7(F17.4)//) 


(D 

105) 

210 

FORMAT (1 XI 31(1 H-)) 


(0 

106) 

211 

FORMAT ( 1X1HI11X1HI31 X6 1 HCORREC TI ON FACTORS FOR 

CORRECTING FROM A 

(0 

107) 


1WIND TUNNEL WHICH IS25X1HI) 


(D 

108) 

212 

FORMAT (LX1HU1X1H1117( IH-I1HI 1 


(0 

109 ) 

213 

FORMAT (1X1HI11X1HI1 6X1H I5X6HCL0SE05X 1HI 16X1HI 2X12HCLOSED FL00R2X1 

(0 

110) 


1HI6X4H0PEN5X1H 1 16X1HI 5X6HCL0SED4X1HI ) 


(D 

111 ) 

214 

FORMAT (1X1HI3X5HDELTA3X1HI5X6HCL0SED5X1HI 4X9H0N 

B0TT0M3X1HI6X4H0P 

(0 

112) 


1EN6X 1HI 6X4H0NL Y6 X1H I 5 X5HFL00R6 X 1H 1 5 X6HCL0S ED5X1 H 13 X9H0N B0TT0M3 X 1 H 

(D 

113) 


21 ) 


(0 

114) 

215 

FORMAT ( 1X1HI11X1HI 1 6X1H I6X4H0NL Y6X 1HI 16X1 8H I ( GRCUND EFFECT) I6X4H 

(D 

115 ) 


10NLY6X 1HI16X1HI6X4 HONL V5 XlH I ) 


(0 

116 ) 

216 

FORMAT ( 1XIHIUX1HI84(1H-)1HI3 2< 1H-I1HI ) 


(0 

117) 

217 

FORMAT ( 1 XI H 1 1 1X1HI 36X1 1HT0 FREE A I R 37X 1HI 8 XI 6HT C GROUNO EFFECT 8X1 

(D 

118 ) 


1HI) 


(0 

119) 

218 

FORMAT < IX 1 3 1 ( 1H-)/) 


(D 

120 ) 

900 

FORMAT ( 11 ,F9.3,5F10.3/3F10*3> 


(0 

121) 

901 

FORMAT (1H1///39X*INTERFERENCE AT TAIL BEHIND SWEPT WING OF FINITE 

(D 

122) 


1 $PAN*///10X*SIGMA (WING) = *F7 .3, 7X*T AIL HEIGHT 

=*F7.3,8XA8,* LOAD 

(D 

123) 
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Appendix D - Concluded 


2ING ! *7X*GAMMA =*F8.3,7X*ZETA=*F7.3//10X*SIGM A (TAIL) =* F7. 3 » 7X*T A I L (D 129) 
3 LENGTH =*F7 .3 , 8X*LAMB0A =*F8 . 3,7X* ALPHA =*F8.3 , 7X*ETA =*F7.3///) (D 125) 
999 STOP (0 126) 

END (0 1271 
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APPENDIX E 


FORTRAN PROGRAM FOR CALCULATING THE AVERAGE WIND-TUNNEL 
INTERFERENCE OVER A SWEPT WING 
CAUSED BY THE PRESENCE OF LIFTING JETS 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2,1, TO RUN ON COC S300 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


, THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 

APPENDIX 0. 


INPUT WILL BE FOUND AT, AND ABOVE, ADDRESS l. NOTE THAT THE REFERENCE 
ORIGIN HAS BEEN CHOSEN AT THE APEX OF THE SWEPT LIFTING LINE. ONLY ONE CON- 
FIGURATION OF JETS CAN BE TREATED PER RUN. ANY NUMBER OF WING CONFIGURATIONS, 
MAY, HOWEVER, BE TREATED IN ONE RUN FOR THIS ONE JET CONFIGURATION. THE FIRST 
VARIABLE REQUIRED (IN FORMAT 1C3) IS 


NJ TOTAL NUMBER OF JETS IN CONFIGURATION 


AS MANY AS 13 JETS CAN BE CONSIDERED BY THE PROGRAM AS LISTED HEREIN. IF 
MORE JETS ARE REQUIRED, AS MANY AS 99 CAN BE OBTAINED BY SUITABLE INCREASES IN 
XH1, YH1 , ZH1 , AND XLOAO IN THE DIMENSION STATEMENT. FURTHER INCREASES 
REQUIRE ALTERATION OF FORMAT 103. 


THE NEXT VARIABLES REQUIRED ARE THE LOCATIONS AND THE RELATIVE STRENGTHS 
OF THE JETS. AS MANY CARDS MAY BE USED AS NEEDED, HOWEVER, THE TOTAL NUMBER OF 
SETS OF DATA MUST AGREE WITH NJ. INPUT VARIABLES FOR THE JETS (IN FORMAT 904) 
ARE 


XH1 


X-C00RDIN4TE OF NOZZLE EXIT AT ZERO 
S IONAL I ZED WITH RESPECT TO THE 


ANGLE OF ATTACK, NONDIMEN- 
TUNNEL SEMIHEIGHT 


YH1 


Y-COORDINATE OF NOZZLE EXIT AT ZERO 
SIONALIZED WITH RESPECT TO THE 


ANGLE OF ATTACK, NONDIMEN- 
TUNAEL SEMIHEIGHT 


ZHl 


Z-COORDINATE OF NOZZLE EXIT AT ZERO 
SIONALIZED WITH RESPECT TO THE 


ANGLE OF ATTACK, NONDIMEN- 
TUNNEL SEMIHEIGHT 


XLOAD RELATIVE PORTION OF LOAD CARRIED BY JET 


SUBSEQUENT TO SPECIFICATION OF THE 
(ONE CARD PER CASE) IN FORMAT 903. THE 


JETS, MODEL CONFIGURATIONS ARE GIVEN 
REQUIRED INPUT VARIABLES ARE 


16 



Appendix E - Continued 


ZETA1 SEN IHE IGHT OF TUNNEL 01 VIDEO 8Y HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMA RATIO OF WING SPAN TO TUNNEL WIDTH 

LAM80A WING SWEEP ANGLE, DEG 

ALPHA ANGLE OF ATTACK OF WING, DEG 


PROGRAM WINDTUNI INPUT, OUTPUT ,TAPE5= INPUT ,TAPE6=0UTPUT) 

COMMON Z ETA, ETA, GAMMA, XOVERH ,YOVERH .ZOVERH, DELTA <28 ) 

CIM6NSICN XD£LTA<28)»XLOAD< 10 », XHll 101 ,YHI( iO> ,ZHlt 10) ,CI 8) 

REAL LAM8DA 

CAT A (C( It, 1 = 1, 81/ 20., 30., 40. ,50. ,60. ,70., 80. ,90./ 

DO 803 Ll=l ,28 
803 XDELTA(L1 1=0. 

RAD=. 0174532925199 
IFIRST=0 

READ <5,103) NJ 

READ <5,904) < XH1< I ) , YH1 < I ) , ZH1 < I ) , XLOAD< I ) , 1=1 , NJ ) 

1 READ <5,903) ZET A1 » ET A1 , GAMMA, SIGMA»LAMBDA» ALPHA 
IF < EOF , 5 ) 999,48 

48 IF < IFIRST .NE.O ) GO TO 47 

WRITE 16,903) < XHl < I ) , YH1< I ) , ZH1 < I ) , XLOAD < I ) , I = 1, NJ ) 

47 WRITE <6,901) GAMMA, ZET Al, LAMBDA , S IGMA , ETAl .ALPHA 
WRITE <6,210) 

WRITE <6,211 ) 

WRITE <6, 212) 

WRITE 16,213) 

WRITE <6, 214) 

WRITE <6,215) 

WRITE <6,216) 

WRITE <6,217) 

WRITE <6,218) 

I F IRST = 1 

L AMBCA=LAMBDA*RAD 
ALPHA® ALPHA*R AO 
SUML=0. 

CO 820 M2=1,NJ 
SUML=SUML ♦XL0AD<M2) 

820 CONTINUE 

SUML =1 ./< 10.*SUML) 

CO 41 K=1 ,8 

IF IETA1.NE.1.) GO TO 813 
M3=0 

DO 815 M 4=1 , N J 
IF <YHHM4).EQ.O. ) GO TO 817 
DO 816 M5=l , NJ 
IF IM4.E0.M5) GO TO 816 

IF < XHKM4) .EQ.XH1IM5) . AND.ZHKM4) .EQ.ZH1IM5) . AND. YH1IM4) . EQ.-YH1 < 
1M5 ) ) GO TO 821 
GO TO 816 

821 M3=M3+1 
816 CONTINUE 


( E 

1 ) 

CE 

2) 

(E 

3) 

CE 

4) 

CE 

5) 

C E 

6) 

(E 

7 ) 

C E 

8) 

(E 

9) 

(E 

10) 

( E 

in 

C E 

12) 

C E 

13) 

( E 

14) 

(E 

15) 

CE 

16) 

C E 

17) 

C E 

18) 

CE 

19) 

(E 

20 ) 

C E 

21 ) 

CE 

22) 

(E 

23) 

(E 

24) 

(E 

25) 

(E 

26) 

C E 

27) 

( E 

28 ) 

(E 

29) 

(E 

30) 

(E 

31) 

(E 

32) 

CE 

33) 

( E 

34) 

(E 

35) 

( E 

36) 

C E 

37) 

CE 

38) 

CE 

39) 

(E 

40) 

(E 

41) 

(E 

42) 

( E 

43) 

(E 

44) 

(E 

45) 
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Appendix E - Continued 



GO TO 815 * 


( E 

46) 

817 

M3=M3 +1 


(E 

47) 

815 

CONTINUE 


(E 

48) 


I F (M3.NE.NJ ) GO TO 813 


(E 

49) 


I ALPHA- 10 H SYMMETRIC 


( E 

50) 


M6=N6= 1 


(E 

51 ) 


M7=5 


( E 

52) 


N 7=NJ 


(E 

53 ) 


CONST 1 = 2 . 


( E 

54) 


IF (SIGMA. NE.O.) GO TO 812 


( E 

55 ) 


M7= 1 


(E 

56) 


C0NST1=IC. 


( E 

57) 


GO TO 812 


( E 

58) 

813 

M6= N6=l 


(E 

59) 


M 7=10 


(E 

60) 


N7=NJ 


( E 

61) 


CONST 1 = 1 • 


(E 

62) 


I ALPHA= 10 HA SYMMETRIC 


( E 

63 ) 


IF (SIGMA. NE. 0. ) GO TO 812 


( E 

64) 


M7=l 


(E 

65) 


CCNST1 = 10 • 


( E 

66) 

812 

CO 801 M 1 =M6 ,M 7 


(E 

67) 


DO 802 N1=N6,N7 


(E 

68) 


YST AR=( 11 .-2 • *FL OA T( Ml ) )/13. 


( E 

69) 


ZET A=ZET Al/ (1 .*ZETA1*( ZHU N1 ) *COS ( ALPHA ) -X HI ( N1 )*S IN ( ALPHA) ) ) 

( E 

70) 


ETA=ETAl-< (1 ./GAMMA) *YHl (Nl > ) 


(E 

71) 


XOVERH= ABS ( YST AR ) *S I GM A*G AMMA*T AN ( LAMBDA ) *30 SC ALPH4J- 

( XHHNIMCOSI 

(E 

72) 


1 ALPHA J)-< ZHU N1 MS INI ALPHA)) 


{ E 

73) 


YOVERH=YST AR*S I GMA*GAMMA-YH1I Nil 


(E 

74) 


ZOVERF=-ABS(YSTAR)*SIGMA*GAMMA*T4NC L4MBD4 1 * S IN ( ALPHA) 

-1ZH1 (Nl MCOS 

( E 

75 ) 


1 ( ALPHA) ) MXHK Nl )*S IN) ALPHA)) 


( e 

76) 


CALL OLTAS <C (K ) ) 


( E 

77) 

+**+* 4 *++******+*+*** $££ APPENDIX Q FOR SUBROUTINE OLTAS ********************** 


CO 805 L 1 = 1 ? 2 8 


( E 

78) 

805 

X0ELTA(L1 )=XD£LT A( LI ) +DE LTA ( LI ) 


(E 

79) 

8C2 

CONTINUE 


(E 

89) 

801 

CONTINUE 


( 6 

81) 


CO 807 L3 =1 » 28 


(E 

82) 

807 

DELTA! L 3 ) = XDEL T4 ( L 3 ) *SUML*CONST 1 


( E 

83) 


WRITE ( 6,149) C(K ),I ALPHA 


( E 

84) 


WRITE 16,150) (DELTA! I) , 1 = 1 ,25,4) 


( E 

85) 


WRITE (6,151) (DELTACI) ,1=2,26,4) 


(E 

86) 


WRITE (6,152) (DELTACI ), 1=3, 27,4) 


(E 

87) 


WRITE (6,153) (DELTACI), 1=4, 28,4) 


( E 

88) 


DC 814 L4=i ,28 


(E 

89) 

814 

X DELTA ( L4 )=0 • 


(E 

90) 

41 

CONTINUE 


( E 

91) 


GO TO 1 


( E 

92) 

103 

FORMAT (12) 


(E 

93) 

149 

FORMAT (1X*CHI = *F 7. 3, 6XA10 ,* JET CONFIGURATION*/) 


(E 

94) 

150 

FORMAT (3X5H(W,L ) 7 ( F 1 7 • 4 ) ) 


(E 

95) 

151 

FORMAT ( 3X5H!U,L )7CF17.4)) 


( E 

96 ) 

152 

FORMAT (3X5H(W,D )7(F17.4)) 


( E 

97) 

153 

FORMAT (3X5H(U,0)7(F17.4 )//) 


( E 

98 ) 

210 

FORMAT (1X131! 1H-) ) 


(E 

99) 

211 

FORMAT (1X1HI11X1HI31X61HC0RRECTI ON FACTORS FOR CORRECTING FROM A 

(6 

100) 


1WIND TUNNEL WHICH IS25X1HI) 


(E 

101) 

.212 

FORMAT UX1HI11X1HI 117! 1H-)IHI) 


< E 

102) 

213 

FORMAT (1X1HI11X1HI 1 6X1HI 5X6HCL0SED5 X1H I 16XIHI2X12HCL0SE0 FL00R2X1 

(E 

103) 


1HI6X4H0PEN6X1H116X1HI5X6FCL0SED4X1HI ) 


(E 

104) 

214 

FORMAT II XIHI3X5HOELTA3X1HI 5X6HCLOSE35XIHUX9HON B0TT0M3X I H 16X4 HOP 

(E 

105) 
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Appendix E - Concluded 


1 EN6 XiH I6X4H0NL Y6X1FI I 5X5HFL00R6X1HI5X6HCL0SED5X1HI3 X9HON BOTTOM3X1H ( E 106) 
21 ) (E 107) 

215 FORMAT ( l XI H 11 IX 1HI 1 6X 1H I 6X4H0NL Y6X 1HI 1 6X1 8H I ( GR CUND EFFECT) I6X4H (E 108) 

10NL Y6X 1 H I 16X1 H 16 X4H0NL Y5 XIH I ) <E 109) 

216 FORMAT I 1 XIH 1 1 IX 1H I 84( 1 H- ) 1 H 132 ( 1H- ) 1 HI ) (E 1 10 ) 

217 FORMAT 1 1 XIH 1 1 1 X 1H 1 3 6X1 1FT0 FREE A I R 37X1HI 8 XI 6H TO GROUND EFFECT 8X1 (E 111) 

1HI) ( E 112 ) 

218 FORMAT 1 1X131 1 1H- )/ ) (E 113) 

900 FORMAT 16F10.3) (E 114) 

901 FORMAT ( 1 HI/ // 37X*AV ERAGE INTERFERENCE OF SEVERAL JETS ON A FINITE (E 115) 

1 SWEPT WING*/// 3 CX*G AMMA =*F7.3,16X*ZETA =*F7. 3 , 15X*L AMBD4 =* ( E 1 16 ) 

2F7.3//30X*SIGMA =*F7 . 3, 16X*ETA =*F 7. 3, 1 5X* AL PH A =*F7.3//) (E 117) 

903 FORMAT < / / ID IX *REL AT I VE */25 X*XOVERH *20X* YO VERH*2 CX*Z0VERH*1 8X (E 118) 

1*STRENGTH*//I24XF7.3,19XF7.3,19XF7.3,19XF6.3/) ) (E 119) 

904 FORMAT I 3F7 .3 , F 5 .3 ,3F7. 3 ,F5 . 3, 3F7 . 3 , F 5. 3 ) (E 120) 

999 STOP (E 1211 

END (E 122) 
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APPENOIX F 


FORTRAN PROGRAM FOR CALCULATING THE DISTRIBUTION OF 
WIND-TUNNEL INTERFERENCE OVER THE 
SPAN OF A SWEPT WING 
CAUSED BY THE PRESENCE OF LIFTING JETS 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN. VERSION 2.1, TO RJN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
.MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT CF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQLAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
APPENOIX Q. 


INPUT WILL BE FOUND AT, AND ABOVE, ADDRESS 1. NOTE THAT THE REFERENCE 
ORIGIN HAS BEEN CHOSEN AT THE APEX OF THE SWEPT LIFTING LINE. ONLY ONE CON- 
FIGURATION OF JETS CAN BE TREATED PER RUN. ANY NUMBER OF WING CONFIGURATIONS, 
MAY, HOWEVER, BE TREATED IN ONE RUN FOR THIS ONE JET CONFIGURATION. THE FIRST 
VARIABLE REQUIRED (IN FORMAT 103) IS 


NJ TOTAL NUMBER OF JETS IN CONFIGURATION 


AS MANY AS 10 JETS CAN BE CONSIDERED BY THE PROGRAM AS LISTED HEREIN. IF 
MORE JETS ARE REQUIRED, AS MANY AS 99 CAN BE OBTAINED BY SUITABLE INCREASES IN 
XH1 , YH1 , ZH1 , AND XLOAD IN THE DIMENSION STATEMENT. FURTHER INCREASES 
REQUIRE ALTERATION OF FORMAT 103. 


THE NEXT VARIABLES REQUIRED ARE THE LOCATIONS AND THE RELATIVE STRENGTHS 
OF THE JETS. AS MANY CARDS MAY BE USED AS NEEDED, HOWEVER, THE TOTAL NUMBER OF 
SETS OF CATA MUST AGREE WITH NJ. INPUT VARIABLES FOR THE JETS (IN FORMAT 904) 
ARE 


XH1 X-COORDINATE OF NOZZLE EXIT AT ZERO ANGLE OF ATTACK, NONDIMEN- 

SIONALIZED WITH RESPECT TO THE TUNNEL SEMIHEIGHT 

YH1 Y-COORDINATE OF NOZZLE EXIT AT ZERO ANGLE OF ATTACK, NONDI MEN- 

S I ONAL I ZED WITH RESPECT TO THE TUNNEL SEMIHEIGHT 

ZH1 Z-COORDINATE OF N3ZZLE EXIT AT ZERO ANGLE OF ATTACK, NONDI MEN- 

SIONAL I ZED WITH RESPECT TO THE TUNNEL SEMIHEIGHT 

XLCAD RELATIVE PORTION OF LOAD CARRIED BY JET 
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Appendix F — Continued 


SUBSEQUENT TO SPECIFICATION OF THE JETS, MODEL CONFIGURATIONS ARE GIVEN 
(ONE CORD PER CASE) IN FORMAT 90D . THE REQUIRED INPUT VARIABLES ARE 

ZETA1 SEMIHEIGHT OF TUNNEL DIVIDED BV HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL DIVIDED BV TUNNEL 

SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMA RATIO OF WING SPAN TO TUNNEL WIDTH 

LAMBDA WING SWEEP ANGLE, DEG 

ALPHA ANGLE OF ATTACK OF WING, DEG 

C EFFECTIVE WAKE SKEW-ANGLE, DEG 

IN SYMMETRICAL CASES THIS PROGRAM COMPUTES THE INTERFERENCE DISTRIBUTION 
OVER ONE SEMISPAN ONLY. THIS PROGRAM REJECTS CASES OF ZFRO SPAN. FDR SUCH 


CASES, 

THE INTERFERENCE IS UNIFORM AND THE VALUES ARE IDENTICAL TO THOSE 

PRO- 


VIDEO 

EY THE PROGRAM Of APPENDIX E. 




PROGRAM WINDTUN1 INPL T ,0L TPUT ,T APE5= I NPUT , T A PE6= OUTPUT ) 

( F 

1 ) 


COMMON ZET A, ET A, GAMMA, XOVERH.YOVERH , ZD VERH, DELTA (28) 

(F 

2) 


CIMENSICN XDEL T A ( 2 8 ) ,XLOAD( 101 ,XH1( lf'I.YHl (l -I.ZHl (10) ,0ATE(2) 

( F 

3 ) 


REAL LAKBCA 

( F 

4) 


DO 80 3 L 1 = 1 ,28 

( F 

5 ) 

, 803 

XDELTAI LI ) =0 . 

(F 

6) 

l 

RAO =.0174532925199 

i F 

7) 


IFIRST=Q 

(F 

3 ) 


READ (5,103) NJ 

(F 

9 ) 


READ (5,904) 1 XHH I) , YH 1 ( I ) , ZH l ( I ) , XLO AD ( I ) , I =1 ,N J ) 

( F 

n ) 

1 

READ (5,900 ) Z ET A1 , ET A1 , GAMMA , S I GM A , L AM BD A , AL PH A, C 

( F 

ii ) 


IF (EOF, 5) 999,48 

( F 

12 ) 

48 

IF ( IFIRST .NE.O ) GO TO 47 

( F 

13 ) 


WRITE 1 6,933) I XH 1 ( I ) , YH1 ( I ) , ZH 1 ( I ) , XL 0 AD i I ) , i = 1 , NJ ) 

( F 

14) 

47 

WRITE (6,901) GAMMA, ZETA1, LAMBDA, SIGMA, ETA! .ALPHA, C 

(F 

IS ) 


WRITE (6,210) 

( F 

16) 


WRITE (6,211) 

( F 

17) 


WRITE 16,212 ) 

(F 

18) 


WRITE 16,213) 

( F 

19) 


WRITE 16,214) 

( F 

2 r ) 


WRITE 16,215) 

(F 

21 ) 


WRITE 16,216) 

{ F 

22 ) 


WRITE (6,217) 

( F 

23) 


WRITE 16,218 ) 

(F 

24) 


I F IRST = i 

( F 

26 ) 


L AMBCA = L AM3DA*RAD 

(F 

26 ) 


ALPHA= ALPHA *R AD 

(F 

27) 


SUML = 0 . 

(F 

28 ) 


CG 8 20 M2 = l »N J 

( F 

29 ) 


SUML=SUML+XL0AD(M2) 

( F 

VM 

8 20 

CONTINUE 

( F 

31 ) 


SUML = 1 . / S L ML 

( F 

32 > 


IF ( SIGMA. NE.C . ) GO TO fill 

( F 

33) 


WRITE (6.5C5) 

< F 

34 ) 


C7LL CAYTIM (DATE) 

( F 

35 ) 
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Appendix F - Continued 


811 


DATE 


) GO TO 813 


GO TO 817 


821 

816 

817 

815 


813 

812 

801 


***** 

8C5 

802 

8C7 


WRITE ( 6 , 906 > 

GO TO 1 
IF (ETA1.NE.1 
H 3=0 

CO 815 M4 = 1 , N J 
IF (YH1 IMA) .E0.0.1 
DO 816 M5=1,NJ 
IF (M4.EQ.M5) GO TO 816 

IF (XH1 (M4> .EQ.XHKM5) .AND.ZH1 ( M4 J.EQ.ZH1 (M5)* AND.YH1 1 M4 ) .EQ.-Y HI ( 
1M5) ) GO TO 821 
GO TO 816 
M3 = M3H 
CONTINUE 
GO TO 815 
M3=M3* 1 
CONTINUE 

IF (M3*NE#NJ ) GO TO 813 
I 4LPH4= 10H SYMMETRIC 
M1CHEK=5 
GO TO 812 
M 1 CH EK = 1 1 

I ALPHA= 1 0H4 SYMME TR I C 
M 1=0 

DO 802 Nl = l » NJ 

Y ST AR= (11 *-2 • *FLOA T (Ml I )/10. 

ZETA=ZET A 1/ ( 1 . ♦ ZET A1 * ( Z HI ( N 1 ) *CQS ( ALPHA )-XHl( N1 >*S I N ( AL PH A )) ) 
ET4=ET4 1-( Cl. /GAMMA )*YH1(N1 ) ) 

XOVERH=ABSl YSTAR) *SIGMA*GAMMA*T AN( LAMBDA )*COS( AL PH A) - ( XH1 ( N 1 ) *C OS ( 

1 ALPHA M-<ZH1(N1 ) *S IN ( ALPHA ) ) 

Y GVERH= YST AR*S I GM A*G AMMA-YH 1 ( N 1 ) 

ZOVERH=-4BS( YSTAR >*SCGM4*G4MM4*T4N( LAMBDA ) * S I N ( ALPH A) - C ZH1 < NL )*C0S 
1( ALPHA HMXHKN1 )*S IN (ALPHA) I 
CALL OLTAS (C) 

4*************** SEE APPENDIX 0 FOR SUBROUTINE DLTAS **************** 
DO 805 L 1 = 1 ? 28 

XDEL TA(L1)=X DELTA! L l ) + DE LTA ( LI) *XLO AD (N 1 ) 

CONTINUE 

CO 807 L3 = l , 2 8 

CELTAU3 )=XDELTA(L3 )*SUML 


WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
CO 814 


(6, 149) 

( 6, 150 ) 
(6,151 ) 
(6,152 ) 
(6,153 ) 
L4=l , 28 


YSTAR , I ALPHA 
(DELTA! I ) , 1=1 ,25,4) 
(DELTA! I), 1=2,26,4) 
(DELTA! I ) , 1=3, 27,4) 
(DELTA(I ), 1=4,28,4) 


814 


103 

149 

150 

151 

152 

153 
210 
211 

212 

213 


JET CONFIGURATION*/) 


XDELT A ( L4 ) =0 • 

M l=Ml*i 

IF (M1.LE.M1CHEK ) GO TO 801 
GO TO 1 
FORMAT (12) 

FORMAT ( //1X*Y/SEMISPAN **F4. 1 , 10X 41 3* 

FORMAT (3X5H!W,L ) 7 ( F 1 7 *4 ) ) 

FORMAT (3X5H(U,L)7(F17.4)) 

FORMAT (3X5H(W,D )7€F1 7.4 ) ) 

FORMAT (3X5H(U, D)7(F17.4)) 

FORMAT ( 1X13 1 ( 1H- ) ) 

FORMAT ( 1X1HI11X1HI31X61HCORRECTION FACTORS FOR CORRECTING FROM A 
1WIND TUNNEL WHICH IS25X1HI) 

FORMAT (1X1HI11X1HI117(1H-)1HI ) 

FORMAT ( 1X1HI11X1HI16X1HI5X6HCL0SED5X1HI 1 6X 1 H I 2 X 12HC L OSE D FL00R2X1 
1 H 1 6X4H0 PEN6X 1 HI16X1HI5X6FCLQSED4X1HI ) 


(F 

36) 

( F 

37) 

( F 

38) 

(F 

39) 

(F 

40 ) 

( F 

41 ) 

(F 

42) 

! F 

43) 

(F 

44) 

(F 

45) 

( F 

46 ) 

!F 

47) 

(F 

48) 

( F 

49 > 

! F 

50) 

(F 

51 ) 

(F 

52) 

(F 

53) 

(F 

54) 

(F 

55) 

(F 

56) 

t F 

57 ) 

( F 

58) 

(F 

59) 

! F 

60 ) 

( F 

61) 

(F 

62) 

(F 

63) 

( F 

64) 

(F 

65) 

(F 

66 ) 

( F 

67) 

( F 

68) 

r ****** 

( F 

69) 

(F 

70 ) 

(F 

71) 

( F 

72) 

( F 

73) 

(F 

74) 

(F 

75) 

(F 

76) 

(F 

77) 

(F 

78) 

(F 

79) 

<F 

80) 

(F 

81 ) 

( F 

82) 

(F 

83) 

(F 

84) 

( F 

85) 

!F 

86) 

! F 

87) 

(F 

88) 

( F 

89 ) 

(F 

90) 

(F 

91 ) 

(F 

92) 

(F 

93) 

( F 

94) 

(F 

95) 


22 


Appendix F - Concluded 


214 FORMAT (1XIHI3X5HDELT43X1HI 5X6HCLOSED5X1 HI 4X9H0N BOTTOM3X 1 H 16X4 HOP (F 96) 

1EN6X1HI6X4H0NLY6XIHI5X5HFL3OR6X1HI5X6HCLOSED5XIH I3X9H0N BOTTOM3X1H (F 97 ) 
21 ) ( F 98) 

215 FORMAT (1X1HI11X1HI16X1HI 6X4H0NL Y6X 1H I 16X1 8H I ( GROUND EFFECT) I6X4H (F 99) 

1CNLY6X1HI 16X1H16X4H0NLY5X1H1 ) CF 100) 

216 FORMAT I IX 1 H 11 IX 1H 1 841 1H- > 1H 1 32 I 1H- ) 1HI ) IF 101) 

217 FORMAT (1X1HI11X1H 136X11 HTO FREE AIR 37X1 HI 8X1 6HT0 GROUND EFFECT 8X1 (F 102) 

1 H I I (F 103) 

218 FORMAT (1X131C1H-)) (F 104) 

900 FORMAT (7F10.3J (F 105) 

901 FORMAT 1 1 HI // / 37X*D I STR I BUT I ON OVER FINITE SWEPT WING CAUSED BY SE (F 106) 

1VERAL JET S*/ /30X*GAMMA = *F7 . 3 , 16X*ZET A = *F 7. 3, 1 5 X*L AMBDA =*F7.3// IF 107) 

230X*S IGMA =*F7.3, 16X*ETA = *F7. 3, 15X*ALPHA = *F7 .3//60X*CHI =* IF 108) 

3F7.3 //) (F 109) 

903 FORMAT 1 10 IX *R ELAT IVE*/25X*X0VERH*20X*Y0 VERH*20X*ZOVERH*18X IF 110) 

l*STRENGTH*//(24XF7.3,19XF7.3,19XF7.3,l9XF6.3//)> (F 111) 

904 FORMAT I 3F7. 3, F5 .3 , 3F7. 3 ,F5 . 3, 3F7 . 3 , F5.3 ) IF 112) 

9C5 FORMAT I////40X*S IGMA EQUALS ZERO - USE AVERAGE INTERFERENCE PROGR (F 113) 

1AM*// ) (F 114) 

999 STOP IF 115) 

END IF 116) 


I 
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APPENDIX G 


FORTRAN PROGRAM FOR CALCULATING THE AVERAGE W IND-T JNNEL 
INTERFERENCE OVER A TAIL 
CAUSED BY THE PRESENCE OF LIFTING JETS 


THIS PROGRAM WAS WRITTEN IN COC FORTRAN, VERSION 2.1, TO RUN ON COC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFORE ME NT I ONE C COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
APPENCIX Q. 


INPUT WILL BE FOUNO AT, ANO ABOVE, ADDRESS 1. NOTE THAT THE REFERENCE 
ORIGIN HAS BEEN CHOSEN AT THE APEX OF THE SWEPT LIFTING LINE. ONLY ONE CON- 
FIGURATION OF JETS CAN BE TREATED PER RUN. ANY NUMBER OF WING CONFIGURATIONS, 
MAY, HOWEVER, EE TREATED IN ONE RUN FOR THIS ONE JET CONFIGURATION. THE FIRST 
VARIABLE REQUIRED (IN FORMAT 103) IS 


N J TOTAL NUMBER OF JETS IN CONFIGURATION 


AS MANY AS 10 JETS CAN BE CONSIDERED BY THE PROGRAM AS LISTEO HEREIN. IF 
MORE JETS ARE REQUIRED, AS MANY AS 99 CAN BE OBTAINED BY SUITABLE INCREASES IN 
XH1, YH1, ZH1, AND XLOAD IN THE DIMENSION STATEMENT. FURTHER INCREASES 
REQUIRE ALTERATION OF FORMAT 103. 


THE NEXT VARIABLES REQUIRED ARE THE LOCATIONS AND THE RELATIVE STRENGTHS 
OF THE JETS. AS MANY CARDS MAY 3E USED AS NEEDED, HOWEVER, THE TOTAL NUMBER OF 
SETS OF DATA MUST AGREE WITH NJ. INPUT VARIABLES FOR THE JETS (IN FORMAT 90A) 
ARE 


XH1 


X-COORDINATE OF 
SIONAL IZED 


NOZZLE EXIT AT ZERO 
WITH RESPECT TO THE 


ANGLE OF ATTACK, NONDIMEN- 
TUNNEL SEMIHEIGHT 


YH1 


Y-COORDINATE OF 
SICNALIZED 


NDZZLE EXIT AT ZERO 
WITH RESPECT TO THE 


ANGLE OF ATTACK, NONDIMEN- 
TUNNEL SEMIHEIGHT 


ZH1 


Z-COORDINATE OF 
SICNALIZED 


NDZZLE EXIT AT ZERO 
WITH RESPECT TO THE 


ANGLE OF ATTACK, NONDIMEN- 
TUNNEL SEMIHEIGHT 


XLOAD RELATIVE PORTION OF LOAD CARRIED BY JET 


SUBSEQUENT TO SPECIFICATION OF THE 
(ONE CARD PER CASE) IN FORMAT 903. THE 


JETS, MODEL CONFIGURATIONS ARE GIVEN 
REQUIRED INPUT VARIABLES ARE 
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Appendix G - Continued 


ZETA1 SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL TAIL LENGTH BEHIND ORIGIN AT ZERO ANGLE OF ATTACK, NON D I ME NS I ON- 

ALIZED WITH RESPECT TO TUNNEL SEMIHEIGHT 

TH TAIL HEIGHT ABOVE ORIGIN AT ZERO ANGLE CF ATTACK, NOND I HENS I ON- 

ALIZED WITH RESPECT TO TUNNEL SEMI HE IGHT 

ALPHA ANGLE OF ATTACK OF WING, DEG 



PROGRAM WINDTUNI INPUT, OUTPUT ,T APE5= I NPUT ,T A PE 6= OUTPUT I 

CGI) 


COMMON Z ETA, ETA, GAMMA, XOVERH .YOVERH .ZOVERH , DELTA 12 8) 

(G 2) 


DIMENSION XDELTAI28) , XLOADI 10) ,XH1( 10 » ,YHI ( 10) , ZHK 10 ) , C I 8 ) 

(G 3) 


DATA ( C(I), 1=1. 8 1/20., 30., AO., 50., 60. ,70., 80. ,90./ 

<G 4) 


DO 803 L 1=1 , 28 

CG 5) 

803 

XCELTAILl )=0 . 

(G 6) 


PI =3.141592653 5 8979 

CG 7) 


R AD= .0174532925199 

(G 8) 


I FIRST=0 

CG 9) 


READ 15,103) NJ 

(G 10) 


READ 15,904) ( XH1 (I) , YH 1 (II ,ZH 11 U , XLOAO(I) , I = 1 ,N J ) 

(G 11) 

1 

READ (5,900) ZE T A 1 , ET A1 ,GAMM A , SI GM AT ,TL ,TH , ALPH A 

(G 12) 


IF (EOF, 5) 999,48 

CG 13) 

48 

I F ( IFIRST.NE.O) GO TO 47 

(G 14) 


WRITE (6,903 ) ( XH1 ( I ) , YH1 ( I ) , ZH1 ( I ) ,XLOAD (I) , I = 1, NJ) 

(G 15) 

47 

WRITE (6,901) GAMMA, ZETAl.TL, SIGMAT, ETAl.TH, ALPHA 

CG 16) 


WRITE (6,213) 

(G 17) 


WRITE (6,211) 

(G 18) 


WRITE (6,212) 

(G 19) 


WR (TE (6,213) 

(G 20) 


WRITE (6,214) 

(G 21 ) 


WRITE (6,215) 

( G 22 ) 


WRITE (6,216) 

(G 23) 


WRITE (6,2171 

CG 24) 


WRITE (6,218) 

CG 25) 


I FI RST = 1 

CG 26) 


LAMBDA= LAMBDA* RAO 

( G 27) 


ALPHA* ALPHA *RAD 

(G 28) 


SUML=0. 

(G 29) 


DC 820 M2=1,NJ 

(G 30) 


SUML= SUML + XL0ADIM2 ) 

(G 31) 

8 20 

CONT INUE 

CG 32) 


SUML= l./(4.*SUML) 

(G 33) 


DO 41 K=1 ,8 

CG 34) 


IF (ETA1.N6.1. ) GO TO 813 

(G 35) 


M3=0 

CG 36) 


DO 815 M4=l , NJ 

(G 37) 


IF (YH1IM4) .EQ.O. ) GO TO 317 

CG 38) 


DO 816 M5=1,NJ 

(G 39) 


IF ( M4.EQ.M5) GO TO 816 

CG 40) 


IF (XH1 (M4l.EQ.XHl (M5) . AND. ZH1 ( M4) . EQ. ZH1 ( M5 ) . AND. YH1 ( M4) . EQ.-YHK 

CG 41) 
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Appendix G - Continued 


1*5) ) GO TO 821 ( G 42 ) 

GO TO 816 (G 43) 

821 M3=M3*1 (G 44) 

816 CONTINUE I G 45 ) 

GO TO 815 (G 46) 

817 M 3=M3+1 (G 47) 

815 CONTINUE (G 48) 

IF (M3.NE.NJ) GO TO 813 <G 49) 

I ALPHA=10H SYMMETRIC (G 50) 

M 6=N6 = 1 IG 51) 

*7=2 (G 52) 

N7=N J (G 53) 

CONST 1 = 2 • (G 54) 

IF ISIGM4T.NE.0. ) GO TO 812 (G 55) 

*7=1 CG 56) 

C0NST1=4. (G 57) 

GO TO 812 IG 58) 

813 M 6=N6= 1 I G 59 ) 

M7=4 (G 60) 

N 7=N J (G 61) 

CONST 1=1. IG 62 ) 

I ALPH4=10HASYMMETRIC I G 63 ) 

IF (SIGMAT.NE.O. ) GO TO 812 (G 64) 

*7=1 (G 65) 

CON ST 1 = 4 . (G 66) 

812 DO 801 M1=M6,M7 IG 67) 

CO 802 N1 =N6 * N7 IG 68) 

YSTAR=(5.-2.*FL0AT(M1) )/4. IG 69) 

ZET4=ZET41/ 1 1 .♦ ZET A 1* I Z HI I Nl ) *COS ( 4L PHA )-XHl INI )*SIN( ALPHA ) ) ) (G 70) 

ETA=ETAl-( ( 1 ./GAMMA) *YHHN1) ) IG 71) 

XOVERH=I TL-XHl (Nl ) )*COS< ALPH4) + (TH-ZHl (Nl )) *SINI ALPHA) IG 72) 

YOVERH= YST AR*S I GM AT *GAMMA-YH 1 1 Nl ) IG 73) 

Z0VERH=ITH-ZH1 (Nl ) )*COS( ALPHA) -(TL-XHl I Nl )) *SIN( ALPH4) IG 74) 

CALL DLTAS ICIK )) ( G 75 ) 

*****4*************** SEE APPENDIX Q FOR SUBROUTINE DLTAS ********************** 
DO 805 L 1 = 1 » 28 IG 76) 

805 XDELTA ( L 1 ) = XDE LTAIL1 )+(DELTA(Ll)*XLQAO(Nl) ) l G 77 ) 

802 CONTINUE I G 78 I 

801 CONTINUE IG 79) 

CO 807 L3 = l *28 (G 80) 

807 CELTAIL3 )=XDELTAIL3 )*SUML*C0NST1 IG 81) 

WRITE 16,149) C I K ) , I ALPHA IG 82) 

WRITE 16,150) (DELTAII ) ,1 = 1,25,4) IG 83) 

WRITE 16,151) (DELTAII) ,1=2,26,4) (G 84) 

WRITE (6,152) (DELTA! I ), 1 = 3,27,4) IG 85) 

WRITE (6,153) (DELTA(I) ,1 = 4,28,4) (G 86) 

DO 814 L4=l ,28 (G 87) 

814 XDEL TA( L4 ) = 0 . (G 58) 

41 CONTINUE (G 89) 

GO TO 1 (G 90 ) 

103 FORMAT (12) I G 91 ) 

149 FORMAT I//1X*CHI =*F 7 .3 , 6 X4 1 0 , * JET CONFIGURATION*/) (G 92) 

150 FORMAT (3X5H (W ,L )7( F 17 .4 ) ) (G 93) 

151 FORMAT (3X5H(U,L )7(F17.4)) (G 94) 

152 FORMAT I 3X5 M W , 0 )7( FI 7.4 1 ) (G 95) 

153 FORMAT ( 3X5HIU , D ) 7( FI 7 .4 ) ) (G 96) 

•210 FORMAT (1X1311 1H-)) (G 97) 

211 FORMAT ( 1X1HI11X1HI31 X6 1 HCORRECT I ON FACTORS FOR CORRECTING FROM A (G 98) 

1WIND TUNNEL WHICH IS25X1HI) (G 99) 

212 FORMAT 1 1 X1H 1 1 IX 1H 1 1 17 1 1 H- > 1H I ) (G 100) 


*213 FORMAT ( 1X1 H 1 1 1 X1H 1 1 6X1 H 1 5 X6HCLQ SE05X 1H 1 1 6X 1 HI 2 X12HCL0SED FLOOR 2X 1 I G 1 01 ) 
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Appendix G - Concluded 


1HI6X4H0PEN6X1HII 6X 1H 1 5 X6 HCL 3SE D4X 1H I ) 

214 FORMAT (1X1HI3X5HDELTA3X1HI 5X6HC L0SE05X1 HI4X9H0N BOTTOM3X 1 H 1 6X4 HOP 
1 EN6X1 H16X4H0NLY6XIH I5X5HFLOOR6X1HI5X6HCLOSED5X1H I3X9HON BOTTOM 3X1H 
21 > 

215 FORMAT ( 1 X 1HI 1 IX 1HI 16X1HI 6X4H0NLY6X 1HI 1 6X1 8HI ( GROUND EFFECT) I6X4H 
10NLY 6X1HI16XIHI6X4H0NLY5X1HI I 

216 FORMAT ( 1X1 H 11 IX 1HI 841 1H- >1H 1 32( 1H- ) 1HI ) 

217 FORMAT ( 1 XI H II IX 1H 136X1 1 HTO FREE AIR37X1HI8X16HTO GROUND EFFECT8X1 
1H I > 


218 FORMAT C1X131C1H-I) 

900 FORMAT I7F13.3) 

901 FORMAT (1H1///37X*AVERAGE INTERFERENCE AT FINITE TAIL CAUSED BY SE 
IVERAL J ETS*//31 X^GAMMA =*F7 . 3 f 10X*ZETA = *F7 .3 »1 OX*TAI L LENGTH/H 
2 =*F6.3//31X*SIGMA(T) =*F7 . 3 , 10X*ET A =*F7. 3. 10X*TA IL HEIGHT/H =* 
2F 6 • 3/ / 57 X*ALPH A =*F9.3//» 

903 FORMAT ( 101X*REL AT IVE*/25X*XOVERH«-20X*YOVERH*20X*ZOVERH*18X 
1*STRENGTH*//(24XF7.3,19XF7.3,19XF7.3,19XF6.3//) ) 

9C4 FORMAT ( 3F7 .3 , F5 .3 . 3F7 . 3 ,F5. 3, 3F7.3 ,F5. 3 ) 

999 STOP 
END 


CG 

102) 

CG 

103 ) 

(G 

104) 

(G 

105) 

(G 

106) 

(G 

107) 

CG 

108) 

(G 

109) 

( G 

110) 

(G 

111) 

CG 

112) 

(G 

113) 

(G 

114) 

<G 

115) 

<G 

116) 

f G 

117) 

(G 

118) 

(G 
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CG 

120) 

(G 
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APPENDIX H 


FORTRAN PROGRAM FOR CALCULATING THE AVERAGE WIND-TUNNEL 
INTERFERENCE OVER A SINGLE ROTCR 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC 6310 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
, APPENCI X Q. 


INPUT WILL BE FOUND AT ADDRESS 1 (ONE CARD PER CASEI IN FORMAT 900. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE CENTER OF THE ROTOR. THE RE- 
QUIRED INPUT VARIABLES ARE 


LI LOAD INDICATOR, LI=1 FOR UNIFORM OISK-LOAD DISTRIBUTION, LI = 2 

FOR TRIANGULAR DISK-LOAD DISTRIBUTION 
ZETA1 SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMA RATIO OF ROTOR DIAMETER TO WIDTH 

ALPHA ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE, DEG 


PROGRAM WINDTUNI INPLT, OUTPUT ,T APE5= I NPUT * T A PE6=0UTPUT ) (H II 

COMMON ZETA, ETA, GAMMA, XO VERH ,Y0 VERH ,ZOVERH , DELTA (28) (H2) 

DIMENSION XDELTA(28I ,PSI(20) ,XLOAD( 20) ,RUNIF(20) ,RTRI A ( 20 » , C ( 8) (H 3) 

CATA (RUNIF(I), 1=1, 20/4*0.2981, 8*0. 6255, 8*0.8921/ (H A) 

DATA (RTR IA( I), I *1,20) /A *0.4386, 8*0. 7 29 6, 8* 0.92 62/ (H 5) 

CATA (C( I), 1=1,8 )/ 20., 30., 40. ,50. ,50. ,70. ,80. , 90./ (H 6) 

P 1 = 3.14159265358979 (H 7) 

R AD=0 .0174532 92 5199 (H8) 

DO 803 L 1 = 1 , 28 ( H 9) 

9C3 XDELT A ( L 1 ) =0 . (H 10) 

PS I (1 )= ( P 1/4. ) IH 11) 

PSII2 ) = 3.*PSI( 1 ) (H 12) 

PSI 13) = 5.*PS I( 1 ) (H 13) 

PSI(4)=7.*PSI ( 1) (H 14) 

PSI(5)=PSI(13)=(PI/8.) IH 15) 

PSI (6) =PS If 14)=3.*PSI(5) ( H 16) 

PSII7)=PSI(15)=5.*PSI(5) (H 17) 

PSI(8) = PSI( 16) = 7.*PSI (5) (H 18) 

PSI(9) = PSK17) = 9.*PSI(5) (H 19) 
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Appendix H - Continued 



PSI ( 10)=PSI< 18) = 11 .*P$I( 51 


(H 

20) 


PSIU1 > = P SI( 19 ) = 13. *PSI ( 5) 


(H 

21 ) 


PSI ( 12 ) =PS IC 20 )=15.*PSI (5) 


( H 

22) 

1 

REAO < 5 * 900 ) LI,ZETA1,ETA1, GAMMA , 

SIGMA, ALPHA 

(H 

23) 


SUML= .002 5 


(H 

24) 


IP ( EOP f 5) 999,47 


(H 

25) 

47 

IF (LI. EC. 11 GO TO 804 


(H 

26) 


I ALPHA=1CFTR IANGULAR 


(H 

27) 


CO 808 M 2 = 1 , 20 


(H 

28) 

80 8 

XLOAD<M2)=RTRIA(M2) 


(H 

29) 


GO TO 160 


(H 

30) 

804 

I ALPHA= 10 H UNIFORM 


( H 

31 ) . 


DO 8C9 M2= 1 , 20 


(H 

32) 

809 

XLOADt M2 )=RUNI F I M2 ) 


(H 

33) 

160 

WRITE (6,901 > SIGMA, 1 4L PH4 , ZE T A1 , ET 41 ,G4MM4 , 4LPH4 

(H 

34) 


WRITE (6,210) 


(H 

35) 


WRITE (6,211) 


( H 

36) 


WRITE (6,212) 


( H 

37) 


WRITE (6,213) 


(H 

38) 


WRITE (6,214) 


(H 

39) 


WRITE (6,215) 


( H 

40) 


WRITE (6,216) 


(H 

41 ) 


WRITE (6,217) 


(H 

42) 


WRITE (6,218) 


( H 

43) 


ALPHA=ALPHA*RAD 


(H 

44) 


00 41 K=1 * 8 


(H 

45) 


M7=N7= 20 


( H 

46 ) 


IF ( SIGMA. NE.O.) GO TO 815 


( H 

47 ) 


M7=N7= 1 


(H 

48 ) 


CONST 1 = 400 • 


( H 

49) 


GO TO 812 


(H 

50) 

815 

IF ( ET A 1 . NE .1 . )G0 TO 813 


(H 

51 ) 


C0NST1 =2 • 


(H 

52) 


GO TO 812 


(H 

53) 

813 

CONS T 1 = 1 . 


(H 

54) 

812 

DO 801 Ml=l , M7 


(H 

55) 


DO 802 N 1=1 , N7 


< H 

56) 


IF ( ET A 1 .NE.l. ) GO TO 811 


(H 

57) 


IF (PSI(Nl).GT.PI) GC TO 802 


(H 

58) 

811 

ETA=ETA1— (XLOADINl >*SIGMA *Sl N( PS l ( Nl >) ) 

( H 

59) 


ZETA=l./( (1./ZETA1 )-(XLOAD( N1I*SIGM4 * SI Nl ALPHA) *COS(PSI( Nl M*GAMM 

(H 

60) 


1A)) 


(H 

61 ) 


XGVERH= SIGMA *GAMMA*COS( ALPH A ) * ( XLO A0( Ml ) *COS( PS II Ml ) ) -XLO AOINl )*C 

( H 

62) 


10SI PSI ( Nl ) 1 ) 


(H 

63) 


YOVERH=S IGMA *GAMMA* (XLO AO I M l ) *S IN l P S l ( Ml ))- XL D ACINI) *SINIPSI(N1)) 

(H 

64) 


1) 


(H 

65 ) 


ZOVERH=-SIGMA *GAMMA*S IN I AL PH A ) * ( XLO A0( M 1 ) *C0 S I P SI I Ml ) )- XLOADINl )* 

(H 

66) 


1 COS I PS 1 1 Nl ) ) ) 


( H 

67) 


CALL DLTAS (C(K>) 


(H 

68 ) 

* 4 *** 4 *************** SEE appendix Q FOR 

SUBROUTINE DLTAS ********************** 


DO 805 Ll=l ,28 


(H 

69) 

805 

XOELT A I L 1 )=XOELT A I LI )+DELTA(Ll ) 


(H 

70) 

802 

CONTINUE 


(H 

71) 

801 

CONTINUE 


(H 

72) 


DO 807 L3=l , 28 


(H 

73) 

807 

DELTA(L3)=XDELTAIL3 ) *SUML*C ONS T l 


(H 

74) . 


WRITE (6,1491 Cl K) 


(H 

75) 


WRITE 16,150) I DELTA 1 1 ) , 1 = 1, 25, A) 


(H 

76) 


WRITE (6,151) (DELTAU), 1-2 ,26,4) 


( H 

77) 


WRITE (6,152) (DELTA! I 1, 1 = 3,27,4) 


(H 

78) ' 


WRITE (6,153) (DELTA! I ) ,1=4,28,4) 


(H 

79) 
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Appendix H - Concluded 


CO 814 L 4 = 1 » 28 (H 80) 

814 XDELTAI L4 )=0. (H 81) 

41 CONTINUE (H 82) 

GC TO 1 (H 83) 

149 FORMAT C//1X*CHI =* F7.3/) <H 84) 

150 FORMAT I3X5H C W , L )7 ( F 17 .4 ) ) (H 85) 

151 FORMAT (3X5HCU.L I7IF17.4)) (H 86) 

152 FORMAT (3X5 H I W , 0 )7 I F 17 .4 ) ) (H 87) 

153 FORMAT ( 3 X5H( U , D ) 7 { F 1 7 .4 ) ) (H 88) 

210 FORMAT (1X131I1H-)) (H 89) 

211 FORMAT (1X1HI11X1HI31X61 HCORREC TI ON FACTORS FOR CORRECTING FROM A (H 90) 

1WIND TUNNEL WHICH IS25X1HI) (H 91) 

212 FORMAT ( 1 XI H 1 1 1 X1HI 1 17 ( 1H- ) 1HI ) (H 92) 

213 FORMAT ( 1 XI H II IX 1H 1 1 6X 1H 1 5X6HC L0SED5 X 1HI 16X 1HI 2X 12 HCLO SEO FL00R2X1 (H 93) 

1HI6X4H0PEN6X1H116X1HI5X6FCL0SED4X1HI) (H 94) 

214 FORMAT ( 1 XI H 13 X5HDELT43 X 1HI 5X6HCL0S E05X 1 H I4X9H0N 80TT0M3XI HI6X4H0P (H 95) 

1EN6XI HI6X4H0NL Y6X1H 1 5X 5HFLOOR6 X 1H 1 5X6HCL 0SE05X1HI3X9H0N 80TT0M3XIH (H 96) 
21) ( H 97 ) 

215 FORMAT 1 1 XI H 1 1 1 X 1H 1 1 6X 1H I6X 4H0NL Y6X 1H 1 1 6X1 8H I ( GR CUND EFFECT) I6X4H (H 98) 

10NL Y6X 1HI16X1HI6X4H0NLY5X1HI ) <H 99) 

216 FORMAT ( 1 X 1 H II 1 X 1H 1 841 1H- ) 1H I 3 2 i 1H- ) 1H I ) (H 100) 

217 FORMAT 1 1 XI H 1 1 1 X 1H 1 36X1 1 HTO FREE A I R 37X IHI 8X1 6HT C GROUNO EFF6CT8X1 <H 101) 

1HI) (H 102) 

218 FORMAT (IX13K1H-)) (H 103) 

900 FORMAT ( 1 1 , F9 . 3 , 6F 10 . 3 ) (H 104) 

901 FORMAT ( 1 HI / // /40X*AV ERAGE INTERFERENCE OVER FINITE SPAN ROTOR*// (H 105) 

140X*SIGM4 =*F6.3, 12XA10,* LOAD ING*//40X*ZET A =*F6 . 3 , 19X *ET A =*F6. (H 106) 

23//40X*GAMM4 =*F6.3» 18X* ALPHA =*F5.1//) (H 107) 

999 STOP (H 108) 

END (H 109) 
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APPENDIX I 


FORTRAN PROGRAM FOR CALCULATING THE DISTRIBUTION OF 
WINO-TUNNEL INTERFERENCE OVER THE 
LATERAL AXIS OF A SINGLE ROTOR 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON COC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REOUI RED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY RECUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE OLTAS WHICH IS GIVEN IN 
APPENDIX Q. 


INPUT WILL BE FOUND AT ADDRESS 1 (ONE CARD PER CASE) IN FORMAT 103. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE CENTER OF THE ROTOR. THE RE- 
QUIRED INPUT VARIABLES ARE 


LI LOAD INDICATOR, LI=l FOR UNIFORM OISK-LOAD DISTRIBUTION, LI=2 

FOR TRIANGULAR DISK-LOAD DISTRIBUTION 

ZETA1 SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMA RATIO OF ROTOR DIAMETER TO TUNNEL WIDTH 

ALPHA ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE, DEG 

C EFFECTIVE WAKE SKEW-ANGLE, DEG 

IN SYMMETRICAL CASES THIS PROGRAM COMPUTES THE INTERFERENCE DISTRIBUTION 
OVER CNE SEMISPAN ONLY. THIS PROGRAM REJECTS CASES OF ZERO SPAN. FOR SUCH 
CASES, THE INTERFERENCE IS UNIFORM AND THE VALUES ARE IDENTICAL TO THOSE PRO- 
VIDED BY THE PROGRAM OF APPENDIX H. 


PROGRAM WINDTUN< INPUT, OUTPUT, TAPE 5 = 1 NPUT ,T APE6=0UTPUT ) (I II 

COMMON ZETA, ETA, GAMMA, XOVERH .YOVERH, ZOVERH, DEL TA (28) (I 2) 

DIMENSION XDELTA(28 ),PSI (2D) ,XL0AD(20) ,RUNIF(20) ,RTRIA(20) (I 3) 

DATA (RUNIF( I) , 1=1, 20)/4*0. 2981, 8*0. 6255, 8*0. 8921/ (I A) 

DATA (RTRIA(I) ,1=1, 20)/4*0. A386,8*0.7296, 8*0. 9262/ (I 5) 

PI=3. 14159265358979 (I 6) 

RAD=0. 0174532925199 (I 7) 

PSI ( 1 ) = { P 1/4.) (18) 
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Appendix I - Continued 



PSI(2) = 3.*PSI< 1 > 

( I 

9) 


PSI(3I = 5.*PSI( 1 ) 

( I 

10 ) 


PSI <4 )= 7* *PS 1(11 

(I 

11) 


PSI<5)-PSI(13)=(PI/8*) 

( I 

12 ) 


PSI(6)=PSI(14)=3#*PSI(5) 

( I 

13) 


P SI ( 7 1 = PSI ( 15 ) =5**PS 1(51 

( I 

14) 


PSim«PSH16l«7.*PSI (51 

( I 

15) 


PSI(9) = PSI( 17)=9.*PSI (5) 

( I 

16 ) 


PSI CIO ) =PS I C 18 ) = 1 1 **PSI 1 5) 

( I 

17) 


PSI(11)=PSIC 19 )=13#*PSI ( 5) 

( I 

18) 


PSI(12)=PSI(20 >=15.*P$I(5) 

( I 

19) 


$UML = .05 

( l 

2 n ) 


DO 805 N2 = l ? 28 

( I 

21) 

/ 805 

XDELTAI N 2 ) = 0. 

( I 

22 ) 

1 

READ (5,103) LI. ZETA1,ETA1, GAMMA, SIGMA, ALPHA, C 

( I 

23 ) 


IF (EOF, 5) 999,700 

( I 

24) 

700 

IF (LI.EQ.l ) GO TO 806 

( I 

25) 

* 

I AL PHA= 10HTR I A NGUL AR 

( I 

26) 


CO 808 M2 = l , 23 

( I 

27) 

8C8 

XL0AD(M2)=RTRI A(M2 ) 

( I 

28) 


GO TO 702 

(I 

29) 

806 

I ALPHA= 10F UNIFORM 

( l 

30) 


DO 809 M 2 = 1 , 20 

( I 

31) 

809 

X LO AO (M2 ) = RUNIF (M2 ) 

(I 

32 ) 

702 

IF (ETA1.NE.1.) GO TO 813 

( I 

33) 


M 1CHEK= 6 

( I 

34) 


GO TO 47 

( I 

35) 

813 

M 1CHEK= 12 

( I 

36) 

47 

WRITE (6,900 SIGMA, I ALPHA ,ZET41 ,ET A1 , GAMMA, ALPHA, C 

( I 

37) 


WRITE (6,210) 

(I 

38) 


WRITE (6,211) 

( I 

39) 


WRITE (6,212 ) 

( I 

40 ) 


WRITE (6,213 ) 

( I 

41 ) 


WRITE (6,214) 

( I 

42 I 


WRITE ( 6,215) 

( I 

43 ) 


WRITE (6,216) 

( I 

44) 


WRITE (6,217) 

( I 

45) 


WRITE ( 6,218) 

( I 

46) 


IF ( S I GM A *N F .0 • ) GO TO 803 

( I 

47) 


WRITE (6,903) 

( I 

48) 


GO TO 1 

( I 

49) 

803 

ALPHA = ALP H4*R A D 

( I 

50 ) 


Ml =0 

(I 

51) 

804 

YSTAR=1 .2-0.2*FL0AT(Ml) 

( I 

52 ) 


CO 800 N 1=1 , 20 

( I 

53) 


ETA=ETA1-(XL0AD(N1 )* SIGMA *SIN( PS I (Nil)) 

( I 

54) 


ZETA=1 •/( (1 */ZET A1 ) — (XLOAD( N1)*SIGMA * SI N< A LPHA ) *COS ( PS I (Ml ) ) *G AMM 

( I 

55 ) 


1 A) ) 

( I 

56 J 


XOVERH= -S I GM A*GAMMA*COS ( ALPHA)* XL0AD(N1 )*COS(PS UNI) ) 

( I 

57) 


YOVERH= SIGMA *GAMM A* (YSTAR-( XLOAO( N1 )*S IN ( PSI (N l) )) ) 

( I 

58) 


ZOVER H= SIGMA ♦GAMM 4*S IN ( AL PHA )* XLOAO( N1 ) *COS ( P SI (N1 ) > 

( I 

59) 


CALL DLTAS (C ) 

( I 

60 ) 

********************* SEE APPENDIX Q FOR SUBROUTINE DLTAS ********************** 


CO 801 N2=l,28 

( I 

61 ) 

80 1 

XDELTA(N2)=XDELTAC N2 I + DELTMN2 ) 

( I 

62 ) 

800 

CONTINUE 

(I 

63) 


DO 8C2 N2=l » 28 

( I 

64 ) 

* 802 

OELT A ( N2 ) = X0ELTA(N2 )*SUML 

( I 

65) 


WRITE U.1A9) YSTAR 

( I 

66 ) 


WRITE 16.150) (DELTA ( 1 1 * 1=1 .25» A) 

( I 

67) 

4 

WRITE 16,1511 (OELTA(I), 1=2,26,4) 

( I 

68 ) 
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Appendix I - Concluded 


WRITE 16,152) (DELTA! I ) ? 1=3,27,4) 

WRITE (6,153) (DELTA ( I ) , 1=4, 28, A ) 

DO 810 N2 = l , 2 8 
810 XDELT A( N2 )=0 .0 
M 1=M1 +1 

IF (M1.LE.M1CHEK ) GO TO 804 
GO TO 1 

103 FORMAT (I1,F9.3,5F10.3) 

149 FORMAT (//10X*Y/R =*F4.1/» 

150 FORMAT ( 3X5H ( W , L )7 ( F 17 .4 ) » 

151 FORMAT (3X5H(U,L )7(F17.4)) 

152 FORMAT ( 3X5H ( W , D )7 ( F 17 . 4 ) ) 

153 FORMAT ( 3X5H( U , 0 )7 ( F17.4 ) ) 

210 FORMAT ( 1X13K 1H- ) ) 

211 FORMAT (1X1HI11X1HI31X61HC0RRECTI ON FACTORS FOR CORRECTING FROM A 
1WIND TUNNEL WHICH IS25X1HI ) 

212 FORMAT ( IX 1 H II 1 X 1H 1 1 1 7( 1 H- ) 1HI ) 

213 FORMAT 1 1 X1H 1 1 IX 1HI 1 6X 1 H I5X6HC L0SED5XIH I 16X 1H I 2X 12HCL0SED FL00R2X1 
1H I6X4H0PEN6X1H I16X1HI5X6 HCL3SED4X1H I ) 

214 FORMAT ( IX 1H 1 3X5HDELTA3 X1HI 5 X6HCL 0SED5X1 H I 4X9H0N B0TT0M3X1 H 16X4 HOP 
1EN6X1HI 6X4H0NLY6X1H 1 5X5HFLOOR6X1HI5 X6HCL0SED5X 1H I3X 9H0N B0TT0M3 XIH 
21 ) 

215 FORMAT (1X1HI11X1HI 16X1H I6X4H0NLY6X 1HI 1 6X1 SHI ( GROUND EFFECT) I6X4H 
10NLY 6X 1H 1 16X1HI 6X4HCNL Y5 XIH I ) 

216 FORMAT ( IX 1H 11 IX 1H I 84( 1H- ) 1H 1 32( 1H- ) 1HI ) 

217 FORMAT ( 1 XIH II IX 1H I 36X1 1 HTO FREE AI R 37X 1H I 8X1 6HT C GROUND EFFECT 8X1 
1 F I ) 

218 FORMAT (1X131I1H-) ) 

900 FORMAT ( 1 HI/ // /31 X *1 NTERFERENCE DISTRIBUTION OVER LATERAL AXIS OF 
1FINITE SPAN ROTOR*// 

240X*S IGMA =*F6.3, 15XA10,* L DAD I NG */ /40<* ZET A =* F6 . 3 , 21 X *ET A =*F7. 
33//40X*GAMMA =*F6. 3 , 19X* ALP HA =*F7.3//56X*CHI =*F8.3//» 

903 FORMAT ( ///40X54HSIGMA EQUALS ZERO USE AVERAGE INTERFERENCE P 

1RCGRAM ) 

999 STOP 
END 




69) 



70) 



71 ) 



72) 



73 ) 



74 ) 



75) 



76) 



77) 



78 ) 



79) 



89 ) 



81 ) 



82) 



83) 



84) 



85) 



86) 



87) 



88 ) 



89) 



90) 



91 ) 



92 ) 



93) 



94) 



95 ) 



96) 



97) 



98) 



99 ) 



100) 



101 ) 



102) 



103) 



104) 


t 
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APPENDIX J 


FORTRAN PROGRAM FOR CALCULATING THE DISTRIBUTION OF 
WIND-TUNNEL INTERFERENCE OVER T FE 
LONGITUDINAL AXIS OF A SINGLE ROTOR 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN* VERSION 2.1. TO RUN ON COC 5300 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 

* 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
^APPENDIX Q. 


INPUT WILL BE FOUND AT ADDRESS 1 (ONE CARD PER CASE) IN FORMAT 1"3. NOTE 
THAT TFE REFERENCE ORIGIN IS CHOSEN TO BE AT THE CENTER OF THE ROTOR. THE RE- 
QUIRED INPUT VARIABLES ARE 


LI LOAD INDICATOR. L I = 1 FOR UNIFORM DISK-LOAD DISTRIBUTION, L 1=2 

FOR TRIANGULAR DISK-LOAD DISTRIBUTION 

ZETAI SEMIHEIGHT CF TUNNEL DIVIDED BY HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FRCM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWI DTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMA RATIO OF ROTOR DIAMETER TO TUNNEL WIDTH 

ALPHA ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE, DEG 

C EFFECTIVE WAKE SKEW-ANGLE, DEG 


INPUT WILL BE FOUND AT ADDRESS 1 (ONE CARD PER CASE) IN FORMAT 103. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE CENTER OF THE ROTOR. THE RE- 
QUIRED INPUT VARIABLES ARE 


LI LOAD INDICATOR, LI=1 FOR UNIFORM DISK-LOAD DISTRIBUTION, L I =2 

FOR TRIANGULAR DISK-LOAD DISTRIBUTION 

ZETAI SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL OIVIDED BY TUNNEL 

SEMIWIDTH 

C-AMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 


« 
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Appendix J - Continued 


SIGMA RATIO OF ROTOR DIAMETER TO TUNNEL WIDTH 

ALPHA ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE. DEG 

C EFFECTIVE WAKE SKEW-ANGLE, DEG 


THIS PROGRAM REJECTS CASES OF ZERO SPAN. FOR SUCH CASES, THE INTERFERENCE 
IS UNIFORM AND THE VALUES ARE IDENTICAL TO THOSE PROVIDEO BY THE PROGRAM OF 
APPENCI X H. 



PROGRAM W I NDTUNI INPUT, OUTPUT ,T APE5= I NPUT ,T APE6= OUTPUT ) 

( J 

1 ) 


COMMON ZETA, ETA, GAMMA, XOVERH, YOVERH ,Z3VERH, DELTA (28) 

( J 

2) 


DIMENSION XDEL TA (28 ) .PSII20) ,XLOAD» 20) ,RUNI F (201 »RTRI A ( 2D I 

( J 

3) 


CATA (RUN IF C I > , 1=1 , 20 > /4 *0. 298 1 » 8*3 . 6255, 8* 0. 892 1/ 

( J 

4) 


DATA (RTR I A ( l 1 , 1*1 , 2 C > /4*0. 4386 ,8*3 . 7296, 8* 0. 9262/ 

( J 

51 


PI = 3. 14159265358979 

( J 

6) 


RAD=0. 0174532925199 

( J 

7) 


PSI(1I=(P 1/4.) 

( J 

8 > 


PSII2 )=3.*PSI(1 ) 

< J 

9) 


PSI(3) = 5.*PSim 

( J 



PSI (4)=7.*PSI( 1 ) 

( J 

111 


PSI(5)*PSI(13)=(PI/8.) 

( J 

12) 


PSI(6) = PS 1(14 )=3.*PSI (5) 

( J 

13) 


PSI(7I=PSI(L5)=5.*PSI(5I 

( J 

14) 


PSI (8>=PSI( 16 l=7.*PS 1(5) 

( J 

15) 


PSK 9) = PSI( 17)=9.*PSI (5) 

( J 

16) 


PSI(10I=PS I( 18 )=11 .*PSI (5) 

( J 

17) 


PSI ( ll)=PSI( 19)=13.*PSI(5) 

( J 

18) 


PSI(12I = PSI(20 1= 15 »*PSl ( 5) 

< J 

19) 


DO 805 N 2=1 , 28 

( J 

20) 

805 

XDEL TA( N2 )=0. 

(J 

21) 

1 

READ (5,103) LI, ZETA1.ETA1, GAMMA, SIGMA, ALPHA, C 

( J 

22) 


IF (EOF, 5) 999,700 

( J 

23) 

700 

SUML= .05 

( J 

24) 


CCNST=1. 

(J 

25) 


IF (LI.EC.l ) GO TO 806 

( J 

26 ) 


I ALPHA= 10HTR I ANGULAR 

t J 

27) 


DO 808 M2=l , 20 

( J 

28) 

80 8 

XL0AD(M2)=RTRIA(M2) 

( J 

29) 


GO TO 47 

( J 

30 ) 

806 

I ALPHA=10H UNIFORM 

( J 

31) 


DO 809 M2=l , 20 

( J 

32) 

809 

XLOAD(M2)=RUNIF(R2 ) 

( J 

33) 

47 

WRITE (6,900) SIGMA, IALPHA, ZETA1 ,ETA1, GAMMA, ALPHA, C 

( J 

34) 


WRITE 16,210) 

( J 

35) 


WRITE (6,211 ) 

( J 

36) 


WRITE (6,212) 

C J 

37) 


WRITE (6,213) 

( J 

38) 


WRITE ( 6,214 t 

( J 

39 ) 


WRITE (6,215) 

( J 

40 ) 


WRITE (6,216) 

( J 

41 ) 


WRITE (6,217) 

( j 

42) 


WRITE (6,218) 

( J 

43) 


IF (SIGMA. NE. 0. ) GO TO 803 

( J 

44) 


WRITE (6,503) 

( J 

45) 


GO TO l 

( J 

46) 

803 

ALPHA=ALPHA*RAD 

( J 

47) 


P1=0 

(J 

48) 
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Appendix J - Concluded 


804 YSTAR = .2*FL34T(Mi)-l .2 (J 49 ) 

CO 80 0 Nl = 1 » 20 (J 50) 

IF ( ET A 1 *NE *1 * ) GO TO 813 (J 51) 

CONST=2 • (J 52) 

IF (PSKN1KGT.PI) GO TO 800 (J 53) 

813 ET4=ETA l- ( XL04D( Nl )*SIGM4 *S IN ( PS I ( N1 ) ) ) (J 54) 

ZETA=1./U1./ZETA1 )-<XLOAD( N1)*SIGMA *SIN( ALPHA ) *COS( PSKN1 ) )*GAMM (J 55) 

14)) IJ 56) 

XOVERH=S IGMA *GAMM A*C0S ( ALPH A ) * ( YST Aft- ( XLOAD! N 1 ) *CQ S (P SI I M 1 ) ) ) ) (J 57) 

YOVERH=-SIGM4 *G4MM4 ♦XL040 ( Nl ) *S IN ( PS II Nl ) ) (J 58) 

ZCVERH= — S IGMA *G AMMA *S IN 1 AL PHA ) * < YS T AR- ( XL 3 ADI N 1 )*COS< PS It Nl ) ) ) ) (J 59) 

C4LL OLT4S (C ) (J 60) 

*4* 4* ******* 44 i ****** SEE APPENDIX Q FOR SU8R0UT I NE DLTAS 4444444444444444444444 
DO 801 N 2 = 1 1 2 8 (J 61 ) 

801 XDELTAIN2 )=XDELT A ( N2 ) +DELT A ( N2 ) 

800 CONTINUE 

CO 802 N2 = i » 28 

802 DELTA! N2 ) = XDELT A( N2 ) *SUML*CONST 
WRITE (6,149) Y STAR 

( DELTA ( I ) , 1 = 1 ,25,4) 

(DELTA!! ), 1=2, 26,4) 

(DELTA! I ) , 1=3,27,4) 

(DELTAU ) ,1=4,28,4) 


WRITE 
WRITE 
WRITE 
WRITE 
CC 810 


(6.150) 

(6.151 ) 
(6,152 ) 
(6,153 ) 

N2=l ,28 

.0 


810 XDELTA! N2 ) =0 - 
Ml = Ml*l 

IF (M1.LE.12 ) GO TO 804 
GO TO 1 

FORMAT ( 1 1 , F9 • 3 , 5F1 C« 3) 


103 

149 

150 

151 

152 

153 
210 
211 

212 

213 

214 


215 

216 

217 

218 
900 


FORMAT ( // 10X* X/R =*F4.1/) 

FORMAT (3X5H(W,L )7(F17.4I) 

FORMAT (3X5H(U#L )7(F17*4)) 

FORMAT (3X5H(W,D 171 F17.4 ) 1 
FORMAT <3X5H(U,0 )7 ( FI 7.4 ) ) 

FORMAT ( 1X131 ( 1 H- ) ) 

FORMAT (1X1HI11X1HI3 1X61HC0RRECTI0N FACTORS FOR CORRECTING FROM A 
1WIND TUNNEL WHICH IS25X1HI) 

FORMAT (1X1HI11X1HU17( 1H-)1HI ) 

FORMAT (1X1HU1X1HI 16X 1H I 5X 6HC L0SED5 X1H I 16X1HI 2X 12HCL0SED FLOOR 2 XI 
1H I6X4H0PEN6 X1H 1 1 6X1HI 5X6HCL 3 SED4X 1H I ) 

FORMAT ( 1X1H 13X5 HD ELTA3X1H I 5X6HCL0SED5X INI 4X9H0N BOTTOM3 X1HI6X4H0P 
1EN6X1HI 6X4H0NL Y6XIHI 5X 5HFL00R6X1 HI5 X6HCL0SED5X 1H I3X9H0N B0TT0M3 X1H 
21 ) 

FORMAT (1X1HI11X1HI 16X1H I6X4H0NL Y6X1HI 16X1 8HI ( GROUND EFFECT) I5X4H 
1CNLY6X1 HI 16X 1H I6X4H0NLY5 X1H I ) 

FORMAT C1X1HU1X1HK84( 1H-) 1HI32C 1H- )1HI ) 

FORMAT ( 1X1HI11X1HI36X11HT0 FREE A IR37X 1HI 8X1 6HTC GROUND EFFECT 8X1 
1H 1 1 
FORMAT 
FORMAT 


(1X13H1H-) ) 

I1H1////28X*! NTERFERENCE DISTRIBUTION OVER LONGITUDINAL AXI 
IS OF FINITE SPAN ROTOR*// 

240X*S IGMA =*F6.3,15XA10,* L3ADING*//40X*ZETA =*F 7. 3 ,2 3X*E TA =*F5.3 
3//40X*GAMMA =*F6 .3 , 1 9X*A LPH A =*F7. 3/ / 5 6X*CH I =*F8.3//) 

903 FORMAT ( /// 41 X54HS I GMA EQUALS ZERO USE AVERAGE INTERFERENCE P 

1RCGRAM) 

999 STOP 
END 


( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 


( 

( 

( 

( 

( 

( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
( J 
t J 
( J 
( J 
( J 
( J 
(J 
( J 


62) 

63) 

64) 

65) 

66 ) 

67) 

68 ) 

69) 

70) 
71 ) 

72) 

73) 

74) 

75) 

76) 

77) 
78 ) 
79) 
80 ) 
81 ) 
82 ) 

83) 

84) 
85 > 
86 ) 

87) 

88 ) 

89) 

90 ) 

91) 

92) 

93) 

94) 

95) 

96) 

97) 

98) 

99) 
100 ) 
101 ) 
102 ) 

103) 

104) 


j 
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APPENOIX K 


FORTRAN PROGRAM FOR CALCULATING THE AVERAGE WIND-TUNNEL 
INTERFERENCE OVER A TAIL 
BEHIND A SINGLE ROTOR 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN* VERSION 2.1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREME N T 10 NE C COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
APPENDIX Q. 


INPUT WILL BE FOUND AT ACDRESS 1 (TWO CARDS PER CASE) IN FORMAT 900. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE CENTER OF THE ROTOR. THE RE- 
QUIRED INPUT VARIABLES FOR THE ROTOR, ON THE FIRST CARD, ARE 


LI LOAD INDICATOR, LI=1 FOR UNIFORM DISK-LOAD DISTRIBUTION, L I =2 

FOR TRIANGULAR DISK-LOAO DISTRIBUTION 

ZETAi SEMIHEIGHT OF TUNNEL CIVIDED BY HEIGHT OF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL OlVIDED BY TUNNEL 

SEMIWICTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMAR RATIO OF ROTOR DIAMETER TO TUNNEL WIOTH 

ALPHA1 ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE, OEG 

THE REQUIRED INPUT VARIABLES FOR THE TAIL, ON THE SECOND CARD, ARE 


SIGMAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL TAIL LENGTH BEHINO ORIGIN AT ZERO ANGLE OF ATTACK, NON D IMENS ION- 

ALIZED WITH RESPECT TO ROTOR RADIUS 

TH TAIL HEIGHT ABOVE ORIGIN AT ZERO ANGLE CF ATTACK, NONO I MENS I ON- 

ALIZED WITH RESPECT TO ROTOR RADIUS 

ALPHA2 ANGLE OF ATTACK OF BODY CARRYING TAIL, OEG 


THIS PROGRAM REJECTS CASES OF ZERO SPAN. 
IN TERMS OF ROTOR RADIUS, SUCH CASES REPRESENT 


SINCE THE EQUATIONS ARE FORMED 
INPUT ERRORS. THE PROGRAM OF 
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Appendix K - Continued 


APPENDIX D CAN BE USED EOR SUCH CASES SINCE THE REPRESENTATION OF THE LIFTING 
SYSTEMS ARE IDENTICAL WHEN THE SPAN IS VANISHINGLY SMALL. 


PROGRAM WINDTUNI INPUT, OUTPUT , T APE 5 = I NPUT , T A PE 6=OUT PUT ) 

COMMON ZETA,ET4,G4MM4,XCVERH,YOVERH,ZOVERH,DELT4(28 ) 

DIMENSION XDELTA(28)*PSI(20)»XLOAD(20),RUNIF(20) ,R TRI A I 20 ) ,C I 81 
DATA IC( I), I=i,8)/2C.,3C.,40. ,50. ,60. ,70.,80.,90./ 

DATA < RUN IF ( I) *I=L, 201/4 *0.29 81, 8*0. 62 55, 8* 0. 8921/ 

DATA IRTRIA( II ,1=1, 201/4*0.4386,8*0. 7296,8*0.926 2/ 

RAD=0. 0174532925199 
PI=3. 14159265356979 
PSK1 1= (PI/4. I 
PSI(2I = 3.*PSI( 1 ) 

PSI (3) = 5.*PS I( 1 I 
PSI (4 ) = 7 . *PS K 1 > 

PSI15)=PS I( 13 )=( PI/8.1 
PSI (6)=PSI( 14)=3.*PSI(5I 
PSI(7)=PS1(151=5.*PSI<5) 

PSK8 l=PSI(16 1=7 .*PSI (51 
PSI (9)=PS K 17) = 9.*PSI (5) 

PSI (10)=PSI ( 18 )=li.*PSI( 5) 

PSI(11)=PSI ( 19)=13.*PSI ( 5» 

PSI(12)=PSI(20 I=15.*PSI (5) 

DO 803 L 1=1 ,28 

803 XDELTA(L1 )=0. 

1 READ (5,9CC) L I ,Z ET A1 ,E T A1 , G AMM A ,S I GMA R , AL PH 41 , SI GM AT , TL , TH , 

1 ALPHA2 
IF (EOF, 5) 999,47 

47 SUML=0.0125 

IF (LI.EQ.l) GO TO 804 
IALPHA=10HTRI ANGULAR 
DO 808 M2=l ,20 
808 XL0AD(M2)=RTRIA(N2) 

GO TO 48 

804 I ALPHA=10HUNIFORM 
CO 809 M2 = l » 20 

8C9 XLOADI M2 )=RUN I F { M2 ) 

48 WRITE (6,901) I AL PH A, ZE TA1 , S I GMAR , TL, AL PHA 1 , E T A 1 , SI GM4T , TH , 

1 ALPHA2, GAMMA 

WRITE (6,213) 

WRITE (6,211 ) 

WRITE ( 6,212 ) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215 ) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218) 

IF ( SIGMAR.NE.O. ) GO TO 800 
WRITE (6,101 ) 

GO TO 1 

800 ALPH1 =ALPHA1*RAD 
4LPH2=4LPH42*RAD 
DO 41 K= 1,8 

IF (SIGMAT.NE.O. ) GO TO 811 

N6=M6=M7=1 

N7=20 

C0NST1=4. 

GO TO 812 

811 IF (ETA1.NE.1.) GO TO 813 


(K 1) 
(K 2) 
(K 3) 
(K 4) 
(K 5) 
(K 6) 
( K 7) 
(K 8) 
(K 9) 
(K 10) 
(K 11) 
( K 12) 
( K 13 ) 
(K 14) 
(K 15) 
(K 16) 
(K 17) 
(K 18) 
(K 19) 
(K 20) 
(K 211 
<K 22) 
(K 231 
IK 24) 
(K 25) 
(K 26) 
(K 27) 
(K 2 R ) 
(K 29) 
(K 30 ) 
(K 31) 
( K 32 ) 
IK 331 
( K 34) 
<K 35) 
(K 361 
(K 37) 
(K 38) 
IK 39) 
(K 40) 
(K 411 
( K 42 ) 
(K 43) 
(K 44) 
( K 45) 
(K 46) 
(K 47) 
(K 48) 
(K 49) 
(X 50) 
IK 51) 
(K 52 ) 
(K 53) 
(K 54) 
( K 55) 
(K 56) 
(K 57) 
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Appendix K 


Continued 


811 IF (ETAl.NE.l.) GO TO 813 
N6=M6 = 1 

M7=2 
N7 = 20 
CONST1 =2 • 

GO TO 812 
813 M6=N6= 1 
M7=4 
N7=20 
CONST 1=1 

812 CO 801 M1=M6,M7 
DO 802 N1=N6,N7 

ETA=ETA1-(XL3A0(N1 ) *SIGM AR*S IN CPS I( N1 M ) 

ZETA=l./( (1./ZETA1 )-XLQAO(Nl > *SI GM AR*GAMMA* S IN ( ALPH1 ) *COS< PSI (N1 ) > 

1 ) 

XGVERH=SIGMAR*GAMMA*( (TL*C3S(ALPH2) ) +( TH*S I Hi AL PH2 ) )-(XLOAD(Nl) *CQ 
1S( ALPHI )*COS (PSI(Nl)ll ) 

XM1 =FLOAT (Ml ) 

YOVERH= S IGMAR*GAMM A*( - ( ( 2.*XMl-5. )/ 4. ) * ( S I GMAT / S IGM AR ) - ( XLO AD( N 1 ) * 
1SIN(PS I(N1 ) ) ) ) 

ZOVERH=SIGMAR*GAMMA*( I TH*COS ( ALPH2 ) ) -(TL *S I N ( AL PH2 > ) + < XLOADINl ) *SI 
IN ( ALPHI I ♦COSIPS I (N 1 M ) ) 

CALL OLTAS IC( K 1 1 

#** 4* 44 ************** SEE APPENOI X Q FOR SUBROUTINE OLTAS *************** 


( K 

58) 

(< 

59) 

IK 

60) 

( K 

61 ) 

(K 

62) 

(K 

63 ) 

I K 

64) 

I K 

65) 

IK 

66) 

IK 

67) 

( K 

68) 

(K 

69) 

(K 

70) 

(K 

71) 

( X 

72) 

(K 

73) 

( K 

74) 

(K 

75) 

(K 

76) 

I K 

77 ) 

(K 

78) 

IK 

79) 

( K 

80) 

******* 


DC 805 LI =1 ,28 

80 5 X DELTA (Li )=X0ELTA(L1 I + DELTA ( L 1 ) 
802 CONTINUE 
801 CONTINUE 

DO 807 L3 = 1 » 28 

807 CEL TAIL 3 ) =XDELTA < L3 ) *SUML*C CNST1 


1C1 

149 

150 

151 

152 

153 
210 
211 


FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 


EQUALS ZERO USE BASIC PROGRAM*) 


IK 
( K 
( K 
(K 
IK 
(K 


81) 

82) 

83) 

84) 

85) 
86 ) 



WRITE 

(6,149) 

C ( K ) 



(K 

87) 


WRITE 

(6,150) 

(DELTACI ) 

,1=1, 

,25,41 

(K 

88) 


WRITE 

(6,151 ) 

( DELTA ( I ) 

» 1=2, 

,26 ,4) 

« K 

89 ) 


WRITE 

(6,1521 

( OELT A ( l 1 

* 1=3 i 

,27,4, 

(K 

90) 


WRITE 

(6,153) 

(DELTAU ) 

,1=4, 

,28,4) 

(K 

91 ) 


CO 814 

LA-1,28 




I K 

92) 

814 

XDELTM L4 1 =0 . 




(K 

93) 

41 

CONT INUE 




(K 

94) 


GO TO 

1 




(K 

95) 


FACTORS FOR CORRECTING FROM A 


(//40X*SIGMA( ROTOR) 

I//5X+CHI =*F 5 • 2/ ) 

(3X5H(W,L )7( F17.4 ) I 
I 3X 5H ( U »L >7(F17*4M 
(3X5H(W,D )7(F17.4>) 

(3X5H(U,D)7(F17*4)) 

(1X131I1H-) ) 

( 1X1 HI 1 IX 1H 1 3 1X61 ECO RR EC T ION 
1WIND TUNNEL WHICH IS25X1HI) 

212 FORMAT I 1 X 1H 1 1 1 X 1H 1 1 1 71 1 h- ) 1HI) 

213 FORMAT (1X1HI11X1HI1 6X1HI 5X6HC L0SED5 X 1HI 16X 1HI 2X 12 HC LOSEO FL00R2X1 
1H I6X4H0PEN6 X1H 116X1HI5X6FCL0SED4X1HI ) 

214 FORMAT ( 1X1H 13 X5H0EL TA3X 1HI 5X6HCL0SED5X1H 14 X9H0N B0TT0M3XI H 1 6X4 HOP 
1EN6X1HI 6X4H3NLY6X1H I 5X5HFL00R6X1H I 5X6HCL0SE D5X 1H 1 3 X9H0N B0TT0M3 X1H 
2 I ) 

215 FORMAT I 1 X1H 1 1 IX 1H I 16X 1H I6X 4H0NL Y 6X 1H 1 1 6X1 8HI I GROUND EFFECT) I6X4H 
10NL Y6X1 H I16X1HI6 X4H0NL Y5 X1H I ) 

216 FORMAT ( 1X1HI1 1X1HI84I IF- )1HI32( 1H- ) 1HI ) 

217 FORMAT ( 1 X1H 11 1 X1H 13 6X1 1HT0 FREE A I R3 7X1 HI 8 XI 6HTC GROUND EFFECT 8X1 
1HI) 

218 FORMAT (1X131I1H-)) 

900 FORMAT ( 1 1 , F9 . 3 , 4F 10 . 3/ 4F 10 . 3) 


(K 
(K 
(K 
(K 
(K 
(K 
( K 
(K 
(K 
(K 
IK 
(K 
(K 
(K 
« 
IK 
IK 
(K 
(K 
(K 
IK 
(K 


96) 

97) 
98 ) 
99) 

ion) 

101 ) 

102 ) 

103) 

104) 

105) 

106) 
107) 
108 ) 
1C 9 ) 
110 ) 
111 ) 
112 ) 

113) 

114) 

115) 

116) 
117) 
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Appendix K - Concluded 


901 FORMAT (1 HI/ // 38X*AV ERAGE INTERFERENCE OVER TAIL BEHIND FINITE-SPA IK 118) 

900 FORMAT 1 1 1 »F9. 3 , AF 10 . 3/A FIO . 3 > IK 119) 

901 FORMAT 1 1 HI / // 38 X* AV ER AG E INTERFERENCE OVER TAIL BEHIND FINITE-SPA IK 120) 

IN ROTOR*/ / 56 XA 10 t* LOAD I NG*/ / 1 5X*ZET A =* F6 . 3 , 9X*SI GMAI ROTOR ) =* IK 121) 
2F6.3 » 9X*T AI L LENGTH/R =*F6. 3 , 9 X* ALPH A ( ROTOR ) =*F 7. 3// 1 5X*E TA =* IK 122) 

3F6.3»9X*SIGMA( TAIL ) =*F6. 3 , 9X*TA IL HE I GHT/ R =*F t. 3 , 9X* AL PH AI BODY ) IK 123) 

A =*F7.3//58X*GAMMA =*F7.3//» IK 12A) 

9 STOP IK 125) 

FND IK 126) 


APPENDIX L 


FORTRAN PROGRAM FOR CALCULATING THE AVERAGE WIND-TUNNEL 
INTERFERENCE OVER TANDEM ROTORS 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RJN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
APPENDIX 0. 


INPUT WILL BE FOUND AT ACCRESS 1 (TWO CARDS PER CASE! IN FORMAT 900. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE CENTER OF THE FRONT ROTOR. THE 
DIAMETERS AND LOAD DISTRIBUTIONS OF THE TWO ROTORS ARE ASSUMED TO BE IDENTICAL. 
THE REQUIRED INPUT VARIABLES FCR THE FRONT ROTOR, ON THE FIRST CARO ARE 

LI LOAD INDICATOR, LI =1 FOR UNIFORM DISK-LCAO DISTRIBUTION, L 1=2 

FOR TRIANGULAR DISK-LOAD DISTRIBUTION 

ZETA1 SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WINO TUNNEL 

SIGMA RATIO OF ROTOR DIAMETER TO TUNNEL WIDTH 

BETA SIDE-SLIP ANGLE, DEG 

ALPHAF ANGLE OF ATTACK OF TIP-PATH PLANE OF FRCNT ROTOR, DEG 

ALPHAS ANGLE OF ATTACK OF BODY CARRYING REAR ROTOR, DEG 

THE REQUIRED INPUT VARIABLES FOR THE REAR ROTOR, CN THE SECOND CARD, ARE 

LRR DISTANCE CF REAR ROTOR BEHIND ORIGIN AT ALPHAB = 0, NONDIM ENS ION- 

ALIZED WITH RESPECT TO ROTOR RADIUS 

HRR HEIGHT OF REAR ROTOR ABOVE ORIGIN AT ALPHAB = 0, NONDI MENS ION- 

ALIZED WITH RESPECT TO ROTOR RADIUS 

ALPHAR ANGLE OF ATTACK OF TIP-PATH PLANE OF REAR ROTOR, DEG 

THIS PROGRAM COMPUTES INDEPENDENTLY THE INTERFERENCE AT EACH ROTOR DUE TO 
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Appendix L - Continued 


ITS OWN PRESENCE AND THE INTERFERENCE AT EACH ROTOR DUE TO THE PRESENCE OF THE 
OTHER ROTOR. IN THE LIMIT, WHEN SIGMA IS ZERO, THE TWO ROTORS ARE COINCIDENT, 
ALL FCUR INTERFERENCES ARE IDENTICAL, AND ONLY ONE SET OF INTERFERENCES IS CAL- 
CULATED. 


PROGRAM W INDTUNI INPUT, OUTPUT ,TAPE5=INPUT ,T A PE 6=0UTPUT ) (L II 

COMMON ZETA, ETA, GAMMA, XO VET H ,YO VERH ,ZOVE RH , DE LT A ( 28 ) (L 2) 

DIMENSION XDELTA(28),PSI(29) ,XLOAO( 20) ,RUNIF(20) ,RTRI A(20) ,C(8) (L 3) 

REAL LRR ( L 

0AT4 I RUN IF ( 1 1 ,1=1, 201/4*0. 2981, 8 *0.6255, 8*0. 8921/ (L 5) 

CATA (RTRIA(I) ,1=1, 20/4*0. 4386, 8*0. 7296, 8*0. 9262/ (L 6) 

DATA (CIII, 1=1, 81/20. ,30. ,40. ,50. ,60. ,70. ,80. ,90./ (L 7> 

PI= 3. 1415926 53 58979 (L 81 

R4D=0. 0174532925199 (L 91 

CO 803 L 1=1 ,28 ( L 10) 

803 X DELTA! LI 1=0. (L 11) 

PSIU ) = (PI/4.) (L 12) 

PSI (2 )=3.*PS 1(1) (L 13) 

PSI (3 ) = 5 .*PS I( 1 ) (L 14) 

PSI(4)=7.*PSI(1) (L 15) 

PSI(5)=PSI(13)=(PI/8.) (L 16) 

PSI(6)=PSI(14)=3.*PSI(5) (L 17) 

PSI(7»=PSI(15)=5.*PSII5> (L 18) 

PSI(8)=PSI(16)=7.*PSI(5) (L 19) 

PSI(9)=PSI(17)=9.*PSI(5) (L 20) 

PSI ( 10 ) =PSI ( 1 8 » = 1 1 .*PS I ( 5 ) (L 21) 

PSI(11)=PSI (19)=13.*PSI ( 5) (L 22) 

PSI(12)=PSI<20 )=15.*PSI(5> (L 23) 

1 READ (5,900) L I , ZETA1 , ET Al , GAMMA , S I GMA ,BET A, ALPH AF , AL PHAB, (L 24) 

1 LRR,HRR, ALPHAR (L 25) 

IF (EOF, 5) 999,47 (L 26) 

47 A ALPF = ALPHAF (L 27) 

A ALPR= ALPHAR (L 28) 

A ALP8=4LPH48 (L 29) 

ABETA=8ET A (L 30) 

ALPHAF=4LPHAF*R AD (L 31) 

ALPHAR=ALPHAR*RAO ( L 32 ) 

ALPHAB = ALPHAB*RAD ( L 33 ) 

BETA=BETA*RAD (L 34) 

WRITE ( 6 , 100 ) (L 35) 

SUML= .0025 ( L 36) 

IF (LI.EC.l) GO TO 804 (L 37) 

IALPHA= 10HTR I ANGUL AR (L 38) 

DO 808 M2 = l , 20 ( L 39 ) 

808 XL0AD(M2)=RTRI A(M2) (L 40) 

GO TO 160 ( L 41 ) 

804 I ALPHA= 10HUNIF0RM (L 42) 

DO 809 M2 = l , 20 (L 43) 

809 XLOAD(M2)=RUNIF(M2 ) (L 44) 

160 NR0T=4 (L 45) 

IF (SIGMA. EQ.O.) NR0T=1 (L 46) 

00 42 I RTR=1 ,NROT (L 47) 

WRITE (6,701) (L 48) 

GO TO (601,602,603,604), IRTR (L 49) 

601 IF (SIGMA. EQ. 0. ) WRITE (6,707) (L 50) 

IF (SIGMA. EQ. 0. ) GO TO 610 (L 51) 

WRITE ( 6,702 ) (L 52) 

GO TO 6 10 (L 53) 

602 WRITE (6,703) (L 54) 

GO TO 610 (L 55) 
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603 WRITE 16,704 ) 

GO TO 610 

604 WRITE (6,705) 

610 WRITE ( 6,706) IALPHA, GA )- KA , S IGMA , LRR , A4LPF, 441 PB , ZET 41 , E T41 , 

1 HRR , A ALP R, ABET A 

WRITE (6,213) 

WRITE (6,211) 

WRITE (6,212) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,213) 

DO 41 K= 1 , 8 
M7=N7=20 

IF (SIGMA. NE.O.) GO TO 815 

M7=N7=1 

C0NST1=400 • 

GO TO 812 

815 IF (ETA1.NE.1..0R.ABETA.NE.O.) GO TO 813 
C0NST1=2. 

GO TO 812 
813 C0NST1=1. 

812 CO 801 Ml = 1 , M7 
DO 802 N1=1,N7 

IF ( ETA1 .NE .1 • .OR. ABET A. NE.O. ) GO TO 811 
IF (PSI1N1 ) .GT .PI ) GO TO 802 
811 GO TO (621,622,623,624), IRTR 

621 ETA=ETA1-XL0AD(N1)*SIGMA*(SIN(PSKN1» )*COS( BETA) + 

1 COSIPS KN1 ))*COS(ALPFAF)*SIN(BETA) ) 

ZETA=ZETAl/ (1.0-XL0AC(N1)*SIGMA*GAMMA*ZETA1*C0S( PSIINl ) ) 

1 *SIN( ALPHAF) ) 

XOVERH=S I GMA*GAMMA*I XLO AC ( M 1 1* ( COS ( PS HM1) >*COS( ALPH AF ) *C3 S < BET A ) 

1 -SINIPSI (Ml ) l*SIN(BET41 ) -XL04D( Nl > * (COS ( PS I( M ) )*COS( ALPHAF ) 

2 *COS ( BET A)— SINCPS I(N1))*SIN(BETA) ) ) 
Y0VERH=SIGMA*GAMMA*(XL04D(Ml)*(SIN(PSI( Ml ) ) *COS ( BE T A) +C0 S ( P S I 

1 (Ml ) )*COS( ALPHAF )*SIN (BETA ) )- XL 3 AD( N1 ) * ( S I N (PS I ( Nl ) ) * 

2 COS(BETA) -CO S( P S I ( Nl ) ) *COS ( AL PHAF ) *S I N ( BET A) )) 

ZOVERH=-S IGMA*GAMMA*SIN( ALPH AF )* ( XLOAD ( Ml ) *COS ( P SI (Ml ) )- 

1 XLO ADI Nl )*COS( PSKNl )) ) 

GO TO 630 

622 ETA=ETA1-XL0AD(N1)*SIGMA*(S IN (PS I (Nl) )*COS(BET A) +COS( PSIINI ) )* 

1 COS( ALPHAR >*SIN(BETA))-SIGMA*SIN(BETA)*(LRR*COS( ALPHAS H- 

2 HRR*SIN( ALPHAS) ) 

ZETA = ZETAl/( i.O-S I GMA^GA MMA *ZE TA1 *( XLOAD(Nl I *COS (PS I ( Nl ) ) * 

1 SI N ( AL PH AR)+LRR*SIN( ALPHAS )-HRR*COS( ALPHAS)) ) 

X0VERH=SIGM4*G4MM4*( XLOAOI Ml ) * (COS ( PS I ( Ml ) ) *COS ( ALPH AF ) *COS 

1 (BET A) -SIN (PS I (Ml ) )*SIN(BETA) )-XLOAO(N 1 ) *( COS( PS I ( Nl ) ) * 

2 COS(ALPHAR)*COSIBETA)-SIN(PSI( Nl ) ) *S IN ( BET A) )-COS ( BET A ) * 

3 ( LRR*COS ( ALPHAS) +HRR*S IN ( ALPHAB ) ) ) 

Y0VERH=SIGM4*G4MM4*( XLC AC ( Ml ) * ( S I N ( PS I (Ml ) ) *COS ( BET A ) + 

1 COS (PS I (Ml ) ) *COS (ALPFAF)*SIN(BETA) )-XLO ADI Nl )*( S IN( PS I ( Nl ) ) 

2 *C0 S (BET A) +COS( PSI (Nl) )*COS( ALPHAR) *SIN (BETA)) -SIN (BETA)* 

3 (LRR*CQS ( ALPHAB ) *HRR *SI N ( ALPHAB ) ) ) 

ZOVERH=-SIGMA*GAMMA*( XLOAOI M1)*C0S( PS I ( Ml ) ) *SI N( ALPHAF)- 

1 XLO AC ( Nl )*COS( PSKNl )) *S IN (ALPHAR) -LRR*SIN( ALPHAS) > 

2 HRR*COS( ALPHAB)) 

GO TO 630 

623 ETA=ETA1-XLDAD(N1)*SIGMA*(SIN(PSI(N1))*C0S(BETA) +C0 S( PSKNl)) 

1 *COS( ALPHAR )*S IN (BETA) )-SIGMA*SIN( BETA) *(LRR*COS( ALPHAB) 


( L 

56) 

(L 

57) 

( L 

58) 

( L 

59) 

(L 

60) 

(L 

61 ) 

(L 

62) 

(L 

63) 

( L 

64) 

(L 

65) 

<L 

66) 

C t 

67) 

(L 

68) 

(L 

69) 

(L 

70) 

(L 

71) 

(L 

72) 

(L 

73) 

( L 

74) 

(L 

75) 

( L 

76) 

<L 

77) 

(L 

78) 

( L 

79) 

i L 

80 ) 

(L 

81 ) 

(L 

82) 

( L 

83) 

(L 

84) 

( L 

85) 

( L 

86) 

( L 

87) 

(L 

88) 

< L 

89) 

(L 

90 ) 

(L 

91) 

( L 

92 ) 

(L 

93) 

( L 

94) 

( L 

95) 

(L 

96) 

( L 

97) 

(L 

98) 

(L 

99) 

(L 

100) 

(L 

101) 

(L 

102) 

(L 

103) 

<L 

104) 

(L 

105) 

(L 

106) 

(L 

107) 

(L 

108) 

(L 

109 ) 

a 

110 ) 

(L 

111 ) 

(L 

112 ) 

(L 

113) 

(L 

114) 

( L 

115) 

(L 

116) 


43 



Appendix L - Continued 


2 +HRR*S IN ( ALPF AB ) ) 

ZETA=ZETAl/( 1.0-SIGMA*GAMM4*ZETA1*< XLOADINl ) *COS (PS I (Nl l»*S IN 
1 ( ALPHAR) +LRR*S INI ALPHAS >-HRR*C3S( ALPHAB I ) ) 

XOVERH= S I GM A*G AMM A* ( XLO ADI Ml ) * ( COS ( PS I ( MI > )*COS( ALPHAR )*C0S 

1 (BETA>-SIN(PSI(M1))*SIN(BETA> )-XLOAD(Nl ) *(CO SIP SI (Nl ) » * 

2 COS (ALPHAR ) *COS( BETA)-SIN(PSI(N1))*SIN(BETA) I) 
YCVERH=SIGMA*GAMMA*( XLOAOI Ml) * I S I Nl PS I ( Ml ) ) *COS < BET A ) +COS 

1 I PS 1 1 Ml ))*COS( ALPHAR >*SIN(BETA) I -XLO AO I N 1 1 *( SIN I P S 1 1 N 1 ) ) 

2 *COS(BETA) +CCSIPSI (Nl ) I *COS( ALPHAR )*S IN (BETA )) > 
ZCV£RH=-SIGMA*GAMMA*SIN( ALPHAR I* (XLOAOIMl) *COS( PSIIM1 ) )- 

1 Xl04D(Ni)*COS(PSI(Nim 

GO TO 630 

624 £TA=ETA1-XL0AD(N1)*SIGM4*(SIN(PSI INI )>*COS(BETA> -*COS I PS I (Nl i » * 

1 COSIALPHAF )*SIN(BETA)I 

2ETA=ZETAl/( l.O-XLOACINl )*S IGM A*GAMMA*Z ETA l*COS ( PS I ( Nl I)* 

1 S IN ( ALPHAF ) l 

XOV£RH= SIGMA*GAMM A*( XLO AO( Ml ) *( COS ( PS I ( Ml )J *COS ( ALPHAR)*C3S 

1 (BETA )-SIN(PSI(Ml))*SIN(BETA) I *COS ( BET A ) *( LR R*CO S ( ALPHABI 

2 *HRR*SIN( ALPHABI )-XLOAO(Nl » * ( COS (PS I ( Nl ) )*COS( ALPHAF)* 

3 COS(BETA)— SIN(PSI(N1))*SIN(BET4))) 

YOVERH= SI GMA*GAMM A* ( XLO ADI M 1 ) *( S INI PS I ( Ml ) ) *C0 S ( BETA) +C0SI P SI 

1 (Ml ))*COS( ALPHAR ) *S INI BETA) ) +S I N( BET A ) * ( LRR*COS (ALPHABI* 

2 HRR*S IN ( ALPH AB ) l-XLCADI Nl)*(SIM(PSI(Nl) )*C0S(BETA) + 

3 COSIPSI (Nl) ) *COS (ALP HAFI*SIN( BETA) I ) 

ZOVERH=-SIGMA*GAMMA*(XLOAD( Ml) *C0S( PS II Ml I ) *SI N ( ALPHAR ) * 

1 LRR*SIN( ALPHAB ) -HRR*COS ( ALPHAB) -XLO AOI Nl )*COS(PSI (Nl ) )* 

2 SIN( ALPHAF) ) 

630 CALL DLTAS (C(K)I 

********************* SEE APPENDIX Q FOR SUBROUTINE DLTAS *************** 
DO 805 Ll=l , 28 

805 XDELT At LI )=XOELTA( L 1 )*OELTA( LI J 
802 CONTINUE 
801 CONTINUE 

00 807 L3 -l > 28 

8C7 CELTAIL 3 )=XDELTA( L3 ) *SUML*C0NST1 
DO 500 LL = 1 *28 

500 DELTA! LL )=D£LT A ( LL )*PI *( SI GM4**2 1 *GAMMA/ 4. 0 
WRITE (6,149) C ( K ) 

(DELTA!! I ,1=1,25,4) 

(DELTA! I ) , 1 = 2, 26,4) 

(OELTA(l) ,1=3,27,4) 

(OELTAI I ), 1=4,28,4) 


WRITE 
WRITE 
WRITE 
WRITE 
00 814 


( 6, 15C ) 
(6,151 ) 

( 6,152 ) 
(6,153 » 
L4=l ,28 


814 X DELTA! L4 1=0. 

41 CONTINUE 

42 CONTINUE 
GO TO 1 

100 FORMAT (1H1/////////59X*START NEW CASE*) 

149 FORMAT (1X*CHI =* F7.3/) 

150 FORMAT ( 3 X5H ( W , l ) 71 FI 7.4 ) ) 

151 FORMAT (3X5H(U,L )7(F17 .4 ) ) 

152 FORMAT ( 3X5HI W , D ) 7( FI 7 .4 ) ) 

153 FORMAT (3X5H(U , D )7( F17 .4 )// ) 

210 FORMAT (1X13K1H-)) 

211 FORMAT (1X1H111X1HI31X61HC3RRECT10N FACTORS FOR CORRECTING FROM A 
1WIND TUNNEL WHICH IS25X1HI) 

212 FORMAT ( 1 XI H 1 1 IX 1H 1 1 17 ( 1H- 1 1HI ) 

213 FORMAT ( 1X1HI11X1HU6X1HI5X6HCL0SE05X1HI 16X1HI 2X12HCLOSED FL00R2X1 
1H I6X4H0PEN6 X 1H I 16X IH I 5X 6HCL0SED4XH I) 

214 FORMAT ( l X1HI3X5HDELT A3X 1HI 5X6HCL0SED5X1HI 4X9H0N 80TT0M3X1H I 6X4H0P 
1EN6X1HI 6X4H0NLY6X1H 1 5X5HFLOOR6X 1HI 5X6HCL0SE D5X1 H I3X9H0N B0TT0M3X IH 


(L 

1 17) 

CL 

118) 

(L 

119) 

(L 

120) 

(L 

121 ) 
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21) (L 177) 

215 FORMAT I 1X1H 1 1 1 X1HI 1 6X1 HI 6X4H0NL Y6X 1HI 16X1 8H I ( GRCUND EFFECT) I6X4H IL 178) 

10NL Y6X 1HI16X1HI6X4H0NLY5X1HI ) IL 179) 

216 FORMAT I 1X1H11 1X1HI 84 ( 1H-) 1 HI3 2I 1H- ) 1HI I (L 180) 

217 FORMAT < 1 XI H 11 1 X1H 1 36X1 1 HTO FREE A I R 37X1HI 8X1 6HTC GROUND EFFECT 8X1 (L 181) 

1HI) (L 182) 

218 FCRMAT I IX 13 1 1 1H- ) / ) <L 183) 

701 FORMAT l 1H1//41X*AVERAGE WIND-TUNNEL INTERFERENCE OVER TANDEM ROTO IL 184) 

IRS*/) (L 185) 

702 FORMAT (49X*EFFECT OF FRCNT ROTOR ON FRONT ROTOR*//) IL 186) 

703 FORMAT I49X*EFFECT OF REAR ROTOR ON FRONT ROTOR*//) IL 187) 

704 FORMAT I48X*EFFECT OF REAR ROTOR ON REAR ROTOR*//) IL 188) 

705 FORMAT I49X*EFFECT OF FRCNT ROTOR ON REAR ROTOR*//) IL 189) 

706 FORMAT I58XA10* LO ADI NG*// 6X*G AMM A =*F7. 3 , l OX*S I GMA =*F6.3, IL 190) 

110X*L I RR ) /R =*F6 .3 * 1 OX* ALPHA IF R ) =*F8. 3, 10X*ALPH At BODY) =* IL 191) 

2F8.3//6X*ZETA =*F7 .3, 10X*E TA =*F6. 3, 10X*H( RR ) /R =*F6 .3 , IL 192) 

310X*ALPHAI RR) =*F8 .3 , 14X*BET A =*F8.3//> IL 193) 

707 FORMAT I26X*F0R SIGMA = C, THE EFFECTS OF THE FRCNT ROTOR ON THE F IL 194) 
1RCNT ROTOR, THE REAR R0T0R*/26X *0N THE FRONT ROT OR, THE REAR ROTOR IL 195) 

2 ON THE REAR ROTOR, AND THE FRCNT ROTOR ON THE*/ 39X*RE AR ROTOR ARE IL 196) 

3 ALL IDENTICAL, AND ARE GIVEN AS FOLLOWS,*//) IL 197) 

900 FORMAT I 1 1 , F 9.3 , 6F 10 . 3/3F10 . 3 ) IL 198) 

999 STOP IL 199) 

END IL 200) 
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APPENDIX M 


FORTRAN PROGRAM FOR CALCULATING THE AVERAGE WIND-TJNNEL 
INTERFERENCE CVER UNLO ADED-ROTOR CONFIGURATIONS 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 OECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
APPENDIX 0. 


INPUT WILL BE FOUND AT AODRESS l (TWO CARDS PER CASE! IN FORMAT 900. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE CENTER OF THE ROTOR. THE RE- 
QUIRED INPUT VARIABLES FOR THE ROTOR, ON THE FIRST CARO, ARE 


LIR ROTOR 0 ISK-LOAD-DI S TR I BUT I ON INDICATOR, LIR=1 FOR UNIFORM LOAD- 

ING, L I R=2 FOR TRIANGULAR LOAOING 

ZETA1 SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HANO WALL DIVIDED BY TUNNEL 

SEMIWIDTH 

GAMMA WIOTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMAR RATIO OF ROTOR OIAMETER TO TUNNEL WIDTH 

ALPHAR ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE, DEG 


THE REQUIRED INPUT VARIABLES FOR THE WING, GIVEN ON THE SECOND CARD, ARE 


LIW WING SPAN-LOAD-DISTRIBUTION INDICATOR, LIW=1 FOR UNIFORM LOADING, 

L I W=2 FOR ELLIPTIC LOADING 

LW CISTANCE OF WING APEX BEHIND ORIGIN, NONDI MENS ION AL I ZE 0 WITH RE- 

SPECT TO ROTOR RADIUS 

HW 01 STANCE OF WING APEX ABOVE ORIGIN, NONOI MENS IONAL I ZEO WITH RE- 

SPECT TO ROTOR RADIUS 

SIGMAW RATIO OF WING SPAN TO TUNNEL WIDTH 
LAMBDA WING SWEEP ANGLE, OEG 

ALPHAS ANGLE OF ATTACK OF BODY CARRYING WING AND TAIL, OEG 
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Appendix M — Continued 


THIS PROGRAM REJECTS CASES OF ZERO SPAN. SINCE THE EQUATIONS ARE FORMED 
IN TERMS OF ROTOR RADIUS, SUCH CASES REPRESENT INPUT ERRORS. 


THIS PROGRAM COMPUTES INDEPENDENTLY THE INTERFERENCE ON THE WING AND ROTOR 
DUE TO THEIR OWN PRESENCE, AS WELL AS THE INTERFERENCE ON EACH DUE TO THE PRES- 
ENCE OF THE OTHER. 


PROGRAM WINDTUNI INPUT, OUTPUT ,T APE 5= 1 NPUT , T A PE6=0 LTPUT > (M 1) 

COMMON ZETA, ETA, GAMMA, XOVERH.YOVERH.ZOVERH, DELTA (28) (M 2) 

DIMENSION XDELTA128) , PS 1(20 ) ,RLOAO< 20) .RUNIFI 20) .RTRIA120) , (M 3) 

1 XLEI1PI»XLOAD(20) ,C(8) <M4) 

REAL LAMBDA, LW (M 5) 

DATA (RUNIFI II ,1=1, 20/4*0.2981 ,8*0.6255,8*0.8921/ (M 6) 

DATA (RTRIAI1), 1=1, 201/4*0. 4386, 8*0. 7296, 8*0. 9262/ (M7) 

DATA (C( 1), 1*1, 81/20. ,30., AO. ,50. ,60. ,70. ,80. ,90./ (M 8) 

P 1=3. 141 59265358979 IM 91 

RAD=0. 0174532925199 (M 10) 

DO 803 L 1 = 1 , 2 8 (Mil) 

803 X DELT A ( L 1 ) =0 . <M 12) 

PS1 (l)=(PI/4.» (M 13) 

PS I ( 2 )=3.*PSI (1 ) ( M 14) 

PSI(3)=5.*PSIl 1 > (M 15) 

PSI(4) = 7.*PSI( 1) (M 16) 

PSI(5)=PSI(13)=(PI/8.) (M 17) 

PSI(6)=PS I( 14 )=3 . *PS I (5 ) (M 18) 

PSI (7>=PSI(15)=5.*PS1(5) (M 19) 

PSI(8)=PSI(16)=7.*PSI(5) (M 20) 

PSI (9) =PS I( 17)=9.*PSI (5) (M 21) 

PSI(10)=PSI( 18)=il.*PSI( 5) (M 22) 

PS I ( 11 )=PSI(19) = 13.*PSI(5) ( M 23) 

PSI(12)=PSI(20 )=15.*PSI (5) (M 24) 

XLE(1)=XLE110)=0.43579 IM 25) 

XLE( 2 ) = XLE( 9 ) =0.71422 (M 26) 

XLE(3) = XLE( 8) =C. 86603 (M 27) 

XLE ( 4) =XLE ( 7) =C. 95394 (M 28) 

XLE ( 5 ) = XL El 6 ) =0.99499 (M 29) 

1 READ ( 5 , 9CC ) LIR.ZETA1.ETA1, GAMMA, S I GMAR, AL PHAR ,L I W,L W,HW, (M 30) 

1 SIGMAW, LAMBDA, ALPHAS IM 31) 

IF (EOF, 5) 999,47 (M 32) 

47 AALPR= ALPHAR ( M 33 ) 

AALPB=ALPHAB (M 34) 

ALAM=L AMBD A (M 35) 

LAMBDA=LAMBDA*RAC (M 36) 

*LPHAR= ALPHAR4R AD (M 37) 

ALPH AB= ALPH AB*R AO (M 38) 

WRITE l 6,100 ( M 39) 

IF (LIR.EC.l) GO TO 804 (M 40) 

I ALPHA=10HTRI ANGUL AR ( M 41 ) 

DO 808 M2=l , 20 (M 421 

808 RLOADl M2 )=RTRIAIM2) (M 43) 

GO TO 806 (M 44) 

804 I ALP HA= 10 HUN I FORM ( M 45 ) 

DO 809 M2=l , 20 ( M 46 ) 

809 RL0ADIM2) =RTR I A( M2 ) (M 47) 

806 IF (LIW.EQ.l) GO TO 852 (M 48) 

I BET A=8FELLIPTIC I M 49) 

DO 851 M3 = l ,10 (M 50) 

851 XLOAD(M3)=XLE(M3) (M 51) 
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Appendix M - Continued 


GO TO 850 

852 I BET A=8 HUN IFORM 
00 853 M3 = l ,10 

853 XL0AD(M3)=1 .0 

850 IF (SIGMAR.NE.O.) GO TO 855 
WRITE < 6, 701 i 

WRITE ( 6,706 ) I ALPH A, I BET A, $ I GMAR , Z ET A l , LW , A 4LPR , SI GM AW , ETA 1 , 
1 HW,44LP8,G4MF4,4L4M 

WRITE (6,2101 
WRITE (6,211) 

WRITE (6,2121 
WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE ( 6,217 ) 

WRITE (6,218) 

WRITE ( 6,707 ) 

GO TO 1 

855 DO 42 I EL EM = 1 ,4 
WRITE (6,701) 

GO TO (601,602,603,604), I EL EM 

601 WRITE (6,702) 

GO TO 61C 

602 WRITE ( 6, 703 ) 

GO TO 610 

603 WRITE (6,704) 

GO TO 61C 

604 WRITE (6,705) 

610 WRITE (6,706) I ALPH A , I BETA, S IGMAR , Z ET A1 f LW, AAL PR ,SI GMA W , ETA 1 , 

1 HW, AALPS, GAMMA, AL AM 

WRITE (6,210) 

WRITE (6,211) 

WRITE (6,212) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218) 

CO 41 K=l,8 

M7=N7= 20 

GO TO (611,612,613,614), I EL EM 
RCTOR ON ROTOR 

611 SUMI=0 *0025 

IF (ETA1.EQ.1.) $UML=C.005 

GO TO 812 

WING ON ROTOR 

612 SUML =0.00630 52 

IF (SIGMAW.EQ.O.) GO TO 615 
IF (LIW.EQ.1) SUML=0* OC 5 
N7=10 

IF (ET41.NE.1 . ) GO TO 812 
SUML=0. 0126104 
IF (LIW.EQ.1) SUML=0.0 10 
N7=5 

GO TO 812 


( M 52 ) 
(M 53) 
(M 54) 
(M 55) 
( M 56) 
(M 57) 
( M 58) 
( M 59) 
( M 60) 
( M 61 ) 
(M 62) 
(M 63) 
( M 64) 
(M 65) 
( M 66) 
( M 67 ) 
(M 68) 
(M 69) 
(M 70) 
(M 71) 
(M 72) 
(M 73) 
( M 74) 
< M 75) 

( M 76) 
(M 77) 
(M 78) 
(M 79) 
( M 80 ) 
( M 81) 
(M 82) 
( M 83) 
(M 84) 
( M 85) 
(M 86) 
(M 87) 
( M 88) 
(M 89) 
(M 90) 
(M 91) 
( M 92) 
(M 93) 
(M 94) 
(M 95) 
(M 96) 
(M 97) 
(M 98) 
( M 99) 
(M 100) 
(M 101) 
(M 102) 
(M 103) 
(M 104) 
(M 105) 
(M 106) 
(M 107) 
(M 108) 
(M 109) 
(M 11") 
(M 111) 
(M 112) 
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Appendix M - Continued 


615 SUML-0 ♦ 05 CM 

N7=l (M 

XLO A0( 1 ) = 1.0 CM 

IF (ETAi.EQ.l.) SUML^O.IO CM 

GO TO 812 (M 

(M 

WING ON WING (M 

(M 

613 SUML =0*0126104 CM 

IF (SIGMAW.EQ.O*) GO TO 616 CM 

IF (LIW.E0.1I SUML= 0*010 CM 

M7=10 (M 

N7=10 CM 

IF IETAl.NE.1.1 GO TO 812 CM 

SUML*0 *0252208 CM 

IF ILIM.EQ.il SUML=0.020 CM 

N7=5 CM 

GC TO 812 CM 

616 SUML=1. CM 

XLOAOC 1 1=1 *0 (M 

M7=*N7*i (M 

GO TO 812 CM 

CM 

ROTOR ON WING CM 

(M 

614 SUML=0 .0050 CM 

IF CETAl.EO.i.) SUML-C.Olt) CM 

M7=10 CM 

IF ( S l GMAW • NE *0 • ) GO TO 812 CM 

M7=l (M 

SUML*0 • 05 CM 

IF (ETAi.EQ.l.) SUML=0. 10 CM 

812 CO 801 Ml-1 ,M7 CM 

00 802 Nl-1 *N7 CM 

XST AR= ( U.-2.*FLOAT(M1H/10. CM 

YST AR-( ll.-2.*FL0ATI Nl) 1/10. CM 

811 GC TO <621,622,623,6241, IELEM CM 

( M 

ROTOR ON ROTOR CM 

CM 

621 IF (ETA1.NE.1.) GO TO 625 CM 

IF (PSI(Nl) .GT.PI) GO TO 802 CM 

625 ET4=ET41-RL04D(N1)*SIGMAR*S INC PS I (Nl ) I CM 

Z ET A=ZET Al/ ( 1 .O-RLQ AC (Nl )*S I GM AR*GAMMA*Z ETA 1 *COS IP SI (Nil) CM 

1 *S INI ALPHAR I I CM 

XCVERH* S IGMAR*GAMM A*CQS( ALPHAR )*(RLQAD(Ml) *COS( P SI (Ml )) CM 

1 •RLOAC(Nn*COS(PSI(Nim CM 

YOVERH* S IGMAR*GAMMA*( RLOADC Ml ) *SI N ( PS I <W 1 ) ) -RL3 A0( Nl 1* SIN ( P SI (Nl U (M 
1 ) (M 

Z0VERH*-S IGMAR*GAMMA*SIN (ALPHAR )* (RLOAOC Ml >*COS< PSI (Ml )) ( M 

1 -RLO AC( Nl ) *COS (PSI(Nll) ) (M 

XLOAOC Nl 1 = 1 .0 (M 

GO TO 630 (M 

(M 

WING ON ROTOR (M 

(M 

622 IF (ETA1.NE.1..0R.SI GMAW.NE *0 • ) G3 TO 627 CM 

IF (PS {(Ml ). GT.PI) GO TO 80 2 (M 

627 ETA=ETAi- YST4R*$ IGMAW (M 

Z ET A*Z ETA 1 / ( 1 .0-S I GM AR*G AMMA*ZET A1*(A0S(YSTAR I * ( SI GMA W/ S I GMAR ) CM 

1 *TAN(LAMBDA)*SIN(ALPHAB)+LW*SIN( ALPHAB ) -HW*CO$( ALPHAB) )) (M 


113) 

114) 

115) 

116) 

117) 

118) 

119) 

120 ) 
121 ) 
122 ) 

123) 

124) 

125) 

126) 

127) 

128) 

129) 

130) 

131) 

132) 

133) 

134) 

135) 

136) 

137) 

138) 

139) 

140) 
141 ) 

142) 

143) 

144) 

145) 

146) 

147) 

148) 

149) 

150) 
151 ) 

152) 

153) 

154) 

155) 

156) 

157) 

158) 

159) 

160) 
161) 
162) 
163) 
164 ) 
165) 
166 ) 

167) 

168) 

169) 

170) 
171 ) 

172) 

173) 
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Appendix M - Continued 


XCVERH= S I GM AR* GAMMA* ( RLO AD( Ml ) *C 0S( AL PH AR >* COS < P S I ( Ml )) 

1 -ABS(YSTAR)*(S IGM AW/ S I GMAR )*T AN ( L AMBD A ) ♦CO S ( ALPHAB) 

2 -LW*COS ( ALPHAB )-HW*SIN( ALPHAB) ) 
YCVERH=SIGMAR*GANMA*(RL0AD(M1)*SIN(PSI(M1) )-YSTAR*( SIGMAW/ 

1 SIGMAR)) 

ZOVERH=-$ I GMAR *GAMMA*( RLOAD ( Ml )*SIN ( ALPHAR ) *COS( PS I (Ml ) ) 

1 -ABS(YSTAR)*(S IGMAW/ SIGMAR )*T AN ( L AM BD A ) *S I N ( AL P HAB ) 

2 -LW*S INI ALPHAB )+HW*COS ( ALPHAB ) ) 

GO TO 630 

WING ON WING 

623 ETA=ETA1-SIGMAW*YSTAR 

ZETA=ZETAl/( 1 •0-SIGMAR*GAMMA*ZETA!*(ABS(YSTAR)*( S IGM AW / SI GM AR ) 
1 *T AN ( L AMSO A ) *S I N ( ALPHAS ) +L W * S I N ( AL PH A B ) -HW * C OS ( ALPHAB) ) ) 
XOVERH=SIGMAW*GAMMA*TAN (L AM BDA ) *COS ( AL PHAB ) *( AB S (X STAR ) 

1 -ABSIYSTAR)) 

YCVER H=0 *2 *S IGM AW*GAMMA* ( FLOAT (Nl)-FLDAT(Ml) ) 

ZOVERH=-S IGMAW*GAMMA*TAN(LAMBDA)*SIN( AL PHAB ) * ( AB S ( XST AR) 
l -ABS(YSTAR)) 

GC TO 630 

ROTOR ON WING 


624 


IF 

IF 


(ETA1.NE.1.) GO 
(PSI(Nl).GT.PI ) 


TO 626 
GO TO 802 

626 ETA=ETA1-RL3AD(N1)*SIGMAR*SIN(PSI(NI) ) 

Z ETA=Z ET A1 / ( 1.0-RLQAC(Ni )*S I GM AR*GAMMA*Z ET A 1 *COS ( PS I ( N1 ) ) 

1 *SIN( ALPHAR) ) 

XCVERH* S I GMAR + GAMM A* ( ABS (XST AR ) * ( S I GMAW/ S I GM AR ) * TAN ( L AM BOA ) * 

1 COS( ALPHAB )+LW*C0S(ALPH4B)+HW*SIN{ALPHAB)-RLG4D(Nl )* 

2 COS (ALPHAR )*C0S(PSI(N1) ) ) 

YOVERH=S I GMAR ♦GAMMA* (XSTAR*(SIGMAW/SlGMAR)-RLOAD (Nl ) *SIN (PS I (N1 ) ) ) 
ZOVERH=-S I GM AR ♦GAMMA ♦( Ab S( X STA R ) * ( S I GMAW/ S I GMAR ) *TAN (L AMBD A ) * 

1 SI N( ALPHAB )+LW*S IN(ALPHAB)-HW*COS( AL PH AB ) -RL 0 ACM N 1 ) * 

2 SIN( ALPHAR) ♦COSt PS I (Nl) ) ) 

XLGAD(N1>=1.0 

630 CALL OLTAS (C(K)) 


(M 
(M 
(M 
(M 
( M 
(M 
(M 
( M 
(M 
(M 
( M 
(M 
( M 
(M 
(M 
( M 
(M 
(M 
( M 
(M 
(M 
(M 
( M 
(M 
(M 
(M 
( M 
(M 
(M 
(M 
(M 
(M 
(M 
( M 
(M 
( M 
(M 
(M 


174) 

175) 

176) 

177) 

178) 
179 ) 
180) 
181 ) 
182) 

183 ) 

184) 

185) 

186 ) 

187) 

188) 
189) 
190 ) 

191) 

192) 

193) 

194) 

195) 

196) 

197) 

198) 

199) 

200 ) 
201 ) 
202 ) 
2 03 ) 

204) 

205) 
2 06 ) 
207 ) 
208) 
209) 
210 ) 
211 ) 


*****44************** SEE APPENDIX 0 FOR SUBROUTINE DLTAS ********************** 



DC 8C5 Ll=l,28 

(M 

212 ) 

805 

XDELT A ( L 1 )=XDELTA(L1 M-DELT41 LI )*XL04D(N1 ) 

(M 

213) 

802 

CONTINUE 

(M 

214) 

801 

CONTINUE 

(M 

215) 


DO 807 L3=I ,28 

( M 

216) 

8C7 

DELTA (L 3)=XDELTA(L3)*SUML 

(M 

217 ) 


WRITE C 6 * 149 ) C(K) 

(M 

218) 


WRITE (6,1501 (DELTA! I) ,1=1 ,25,4) 

(M 

219) 


WRITE (6,151) (DELTA! I ), 1 = 2,26,4) 

IM 

220) 


WRITE (6,152) (DELTA! I ), 1 = 3, 27,4) 

( M 

221) 


WRITE (6,153) (DELTA! I ), 1=4,28,4) 

(M 

222 ) 


DO 814 L4=l ,28 

(M 

223) 

814 

XDELT A ( L4 ) =0 . 

(M 

224) 

41 

CONT INUE 

(M 

225) 

42 

CONTINUE 

(M 

226 ) 


GO TO 1 

(M 

227) 

100 

FORMAT (1H1/////////59X*START NEW CASE*) 

(M 

228) 

149 

FORMAT (1X*CHI =* F7.3/) 

(M 

2 29) 

150 

FORMAT <3X5H(W,L )7(F17.4) ) 

(M 

230) 

151 

FORMAT <3X5H(U,L )7(F17.4>) 

(M 

231 ) 

152 

FORMAT (3X5H(W,D)7(F17.4 )) 

(M 

232) 

153 

FORMAT (3X5H(U,D)7(F17.4)//) 

(M 

233) 
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210 FORMAT I1X131I1H-)) ( M 

211 FORMAT ( 1X1HI11X1HI31X61HCORRECTION FACTORS FOR CORRECTING FROM A (M 

1WIND TUNNEL WHICH IS25X1HII IM 

212 FORMAT ( 1X1H II 1X1H 111 7( 1H- ) 1HI ) IM 

213 FORMAT (1X1HI11X1HI1 6X1HI5X6HCL0SED5X 1HI 16X 1HI 2 X12HC LOSEO FL00R2X1 IM 

lFi I6X4HQPEN6X1H I 16X 1HI 5X6HCL0SED4X 1HI I IM 

214 FORMAT 1 1 X1H I3X5H0ELT A3X 1HI 5X6HCL0SED5X1HI4X9H0N B0TT0M3X1 H 1 6X4 HOP (M 

1EN6X1HI6X4H0NLY6X1HI 5X5HFL00R6X1HI 5X6HCL0SED5X1H I3X9H0N B0TT0M3 X1H IM 
21) IM 

215 FORMAT I 1X1H 11 1X1H11 6X1HI6X4H0NL Y6X IHI 16X1 8H II GR CUND EFFECT ) ISX4H IM 

10NLY6X1HI16X1HI 6X4H0NLY5 X1H I ) IM 

216 FORMAT 1 1 X1H II 1 X 1H I 841 1H- ) 1H 1 3 21 1H- ) IHI ) IM 

217 FORMAT tlXIHIl 1X1HI36X11HT0 FREE AIR37X1HI8X16HT0 GROUND EFFECT8X1 IM 

IHI ) IM 

218 FORMAT 1 1X131 1 1H- ) /I IM 

701 FORMAT 1 1H1/ /42X*AVERAG£ INTERFERENCE OVER AN UNLOADED ROTOR MODEL IM 

1*/) IM 

702 FORMAT !54X*EFFECT OF RCTOR ON ROTOR*//) (M 

703 FORMAT I55X*EFFECT OF WING ON ROTOR*//) IM 

704 FORMAT !55X*EFFECT OF WING ON WING*//) IM 

705 FORMAT (55X*EFFECT OF ROTOR ON WING*//) IM 

706 FORMAT I34X.A10* ROTOR L0ADING*19X , A8* WING LOADING*// IM 

119X* SIGMA! ROTOR) =*F6.3 , 10X*ZETA =*F6.3 ,10X*LW/R =*F6.3,13X IM 

2*ALPHA! ROTOR) =*F7.3//19X* SIGMA! WING) = *F 6 .3, 1 OX*£TA =*F6.3, IM 

310X*HW/R =*F6. 3 » 10X* ALPH AI BODY ) =*F7. 3//39X*GAMMA =*F6.3, (M 

427X*LAMB0A =*F7.3//) IM 

707 FORMAT I 40X*S IGMAI ROTOR ) EQUALS ZERO, THIS PROGRAM IS NOT SUITABLE IM 

1 FOR USE WITH SUCH CASES.*//) IM 

900 FORMAT I I1»F9. 3, 4F 10.3/1 1*F9.3,4F10.3) IM 

999 STOP IM 

END IM 


234) 

235) 

236) 

237) 

238) 
2 39) 
240 ) 
241) 
242 ) 

243) 

244) 
245 ) 

246) 

247) 

248) 
249 > 
250) 
251 ) 

252) 

253) 
2 54) 

255 ) 

256 ) 
2 57) 

258) 

259) 

260) 
261) 
262) 
263 ) 
264) 
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FORTRAN PROGRAM FOR CALCULATING THE AVERAGE WIND-TUNNEL 
INTERFERENCE OVER A TAIL 
BEHIND AN UNLOADE O-ROTOR CONFIGURATION 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE DLTAS WHICH IS GIVEN IN 
APPENDIX Q. 


INPUT WILL BE FOUND AT ADORESS 1 (THREE CARDS PER CASE! IN FORMAT 900. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN TO BE AT THE CENTER OF THE ROTOR. THE RE- 
QUIRED INPUT VARIABLES FOR THE ROTOR, ON THE FIRST CARO, ARE 


LIR ROTOR DISK— LOAD -DISTRIBUTION INDICATOR, LIR=1 FOR UNIFORM LOAD- 

ING, L IR = 2 FOR TRIANGULAR LOADING 

ZETA1 SEMIHEIGHT CF TUNNEL DIVIDED BY HEIGHT OF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FROM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWIDTH 

C-AMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMAR RATIO OF ROTOR DIAMETER TO TUNNEL WIDTH 

ALPHAR ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE, DEG 


THE REQUIRED INPUT VARIABLES FOR THE WING, GIVEN ON THE SECOND CARD, ARE 


LIW WING SPAN-LOAD-DISTRIBUTION INDICATOR, LIW=1 FOR UNIFORM LOADING, 

L I W = 2 FOR ELLIPTIC LOADING 

LW DISTANCE OF WING APEX BEHIND ORIGIN, NOND I MENS I ON AL I Z E D WITH RE- 

SPECT TO ROTOR RADIUS 

HW DISTANCE OF WING APEX ABOVE ORIGIN, NONOI MENSION AL I ZED WITH RE- 

SPECT TO ROTOR RADIUS 

SIGMAW RATIO OF WING SPAN TO TUNNEL WIDTH 
LAMBDA WING SWEEP ANGLE, DEG 
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Appendix N - Continued 


ALPHAS ANGLE OF ATTACK OF BODY CARRYING MING AND TAIL, OEG 

THE REQUIREO INPUT VARIABLES FOR THE TAIL, ON THE THIRD CARO, ARE 
SIGHAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL TAIL LENGTH BEHINO ORIGIN AT ZERO ANGLE OF ATTACK, NOND I ME NS I ON- 

ALIZEO WITH RESPECT TO ROTOR RADIUS 

TH TAIL HEIGHT ABOVE ORIGIN AT ZERO ANGLE OF ATTACK, NONDIMENSION- 

ALIZED WITH RESPECT TO ROTOR RADIUS 

THIS PROGRAM REJECTS CASES OF ZERO SPAN. SINCE THE EQUATIONS ARE FORMED 
IN TERMS OF ROTOR RADIUS, SUCH CASES REPRESENT INPUT ERRORS. 

THIS PROGRAM COMPUTES INDEPENDENTLY THE INTERFERENCE AT THE TAIL DUE TO 
THE PRESENCE OF BOTH THE WING AND THE ROTOR. 


PROGRAM WINDTUNI INPUT, OUTPUT ,T APE5= I NPUT , T A PE 6=0LTPUT) (N l) 

COMMON ZETA, ETA, GAMMA, XC VERH ,YOVERH , ZOVERH, DELTA (28) IN 2) 

DIMENSION XDEL TA128I ,PSI(20l,RLOAD(2O), RUN I F { 2 0 1 ,R TR I A( 20 ) » IN 3) 

I XLEUOI.XLOADI 20) ,CI 81 IN A) 

REAL LAMBDA, LW IN 5) 

DATA IRUNIFII) ,1=1, 20>/4*0. 2981, 8*0. 6255, 8*0. 8921/ IN 6) 

DATA IRTRIAII), 1=1, 20)/4*0. 4386, 8*0. 7296, 8*0. 9262/ IN 7) 

DATA I C I I ) , I=l,8)/20. ,30. ,40. ,50. ,60. ,70. ,80. ,90./ IN 8) 

PI=3. 14159265358979 IN 9) 

RAD*0.C174532925199 IN 1") 

CO 803 L 1=1 , 28 IN 11) 

803 X0ELTAIL1 )=0. IN 12) 

PSI II ) = I P 1/4 . ) IN 13) 

PSI I2)=3.*PS III ) IN 14) 

PSII3)=5.*PSII 1) IN 15) 

PSI(4S=7. *PS III) I N 16) 

PSI I5)=PSII13)=(PI/8.) IN 17) 

PSII6) = PSII14)=3 . *PS 1 1 5 ) IN 18) 

PSI I 7 )=PS II 15 )=5 . *PS I 15) IN 19) 

PSI IS ) =PS II 16 >=7.*PS II 5) IN 20) 

PSI19>=PSII17)=9.*PSII5) IN 21) 

PSH10)=PSI! 18 ) = 11.*PSI 15) IN 22) 

PSI 111) =PSI! 19 )=13.*PSI1 5) IN 23) 

PSII12)=PSI120)=15.*PSl 15) IN 24) 

XL El 1)=XLE110)=0. 43579 I N 25 ) 

XLE I2)=XLEI9) =0.71422 IN 26) 

XLEI3 l=XL£l 8) =0.86603 IN 27) 

XLEI4) =XL El 7 ) =0.95394 IN 28) 

XLEI 5I = XLE 1 6 ) =0.99499 IN 29) 

1 READ 15,900) L IR, ZETA1 , ETA1, GAMMA, S IGMAR, ALPHAR ,LI W,L W ,HW, 1 N 30) 

1 SIGMAW, LAMBDA, AL PH AB , SI GM A T , TL ,TH IN 31) 

IF 1 EOF , 5 I 999,47 IN 32) 

47 AALPR= ALPHAR I N 33 ) 

A ALP B= ALP HA B IN 34) 

ALAM=L AMBDA IN 35) 

L AMBDA=LAMBDA*RAD IN 36) 

ALPHAR= AL PH AR*R AD IN 37) 
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ALPHAB= ALPH4B*R4D 


( N 

38) 


WRITE (6,100) 


( N 

39 ) 


IF (LIR.E0.1I GO TO 

804 

(N 

40) 


I ALPHA= 10HTR I ANGULAR 


( N 

41 ) 


CC 808 M2 = l * 20 


(N 

42) 

808 

RLOAD (M2)=RTRI4(M2) 


(N 

43) 


GC TO 806 


( N 

44) 

8C4 

I ALPH A= 10 HUN IF ORM 


(N 

45 ) 


DO 809 H2=l , 20 


(N 

46) 

809 

RLQA0(M2)=RTRIA(M2) 


( N 

47) 

806 

IF (LIW.EQ.l) GO TO 

852 

(N 

48) 


IBETA=8FELLIPTIC 


< N 

49) 


DO 851 M 3= 1 , 1 0 


(N 

50 ) 

851 

XL0A0(M3)=XLE(M3) 


(N 

51) 


GO TO 85C 


( N 

52 ) 

852 

I BETA=8HUN IFORM 


( N 

53) 


DO 853 M3 = 1 ,10 


(N 

54) 

853 

XLOAD* M3) = l .0 


( N 

55) 

850 

IF (SIGMAR.NE.O. > GO 

TO 855 

(N 

56) 


WRITE (6,731) 


(N 

57) 


WRITE ( 6,706) I ALPH A , 

I BETA, SIGMAR,ZETA1 ,LW, A ALPR , SI GMA W , ET A1 , 

( N 

58) 


1 HW, A4LP8,S IGM4T,TL ,TH, GAMMA, 4L AM 

( N 

59) 


WRITE (6,210) 


(N 

60) 


WRITE ( 6,211 ) 


(N 

61 ) 


WRITE (6,212 ) 


t N 

62) 


WRITE (6,213) 


(N 

63) 


WRITE ( 6,214) 


( H 

64) 


WRITE (6,215) 


(N 

65) 


WR PTE (6,216) 


< N 

66) 


WRITE (6,217) 


( N 

67) 


WRITE (6,218) 


(N 

68 ) 


WRITE (6,707) 


(N 

69) 


GO TO 1 


( N 

70) 

855 

CO 42 IEL EM = l ,2 


( N 

71) 


WRITE (6,731) 


(N 

72) 


GO TO (601,602 ), IELE* 

(N 

73) 

601 

WRITE (6*702) 


(N 

74) 


GO TO 61 C 


( N 

75) 

602 

WRITE (6,703 > 


( N 

76) 

610 

WRITE 16,706) I ALPHA , 

IBETA,S IGMAR,ZETA1,LW, AALPR , S I GMAW , E TA 1 , 

( N 

77) 


I HW, AALPB, SIGMA T,TL , TH, GAMMA, ALA M 

(N 

78) 


WRITE (6,213) 


(N 

79) 


WRITE (6,211) 


( N 

80) 


WRITE (6,212) 


(N 

81 ) 


WRITE (6,213) 


(N 

82 ) 


WRITE (6,214) 


(N 

83) 


WRITE (6,215) 


< N 

84) 


WRITE (6,216) 


(N 

85) 


WRITE (6,217) 


( N 

86) 


WRITE (6,218) 


( H 

87) 


DO 41 K = 1 , 8 


(N 

88) 


N7=20 


( N 

89) 


*7=4 


< N 

90) 


GO TO (611,613), IELEM 

(N 

91 ) 




(N 

92) 


EFFECT OF ROTOR 


( N 

93) 




( N 

94) 

611 

SUML=0 .0125 


(N 

95) 


IF (ETA1.EQ.1.) SUML 

=0.025 

( H 

96) 


IF (SIGMAT.NE.O.) GO 

TO 812 

( N 

97) 


SUML=C *C50 


( N 

98) 
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IF (ETA1.E0.1.) SUML =C. 103 



( N 

1 99 ) 



M7= 1 



( N 

100) 



GO TO 812 



(N 

101 ) 

c 





(N 

102) 

c 


EFFECT OF WING 



(N 

103) 

c 





(N 

104) 


613 

SUML =0.031526 



<N 

105 ) 



IF « SIGMAW.EQ.O. .OR . S IGMAT . EQ.O . ) GO 

TO 616 

CN 

106) 



IF (LIW.EQ.II SUML=0.025 



(N 

107) 



M7=4 



(N 

108 ) 



N7=10 



(N 

109) 



IF (ET41.NE.1.I GO TO 812 



(N 

110) 



SUML =0.063052 



CN 

til ) 



IF (LIW.EQ.U SUML=0.050 



(N 

112 ) 



N7=5 



CN 

113) 



GO TO 812 



(N 

114) 


616 

IF (SIGMAW.EO.O.. AND. SI GMAT .NE .0 . ) GO 

TO 

612 

(N 

115) 



IF (SIGMAW.NE.O. .ANO.SIGMAT.EO.O. ) GO 

TO 

614 

(N 

116) 



SUML=l • 0 



CN 

117) 



XLOAOC 1 »-1.3 



CN 

118) 



M7=N7=1 



(N 

119) 



GO TO 812 



CN 

120) 


612 

SUML =0.25 



(N 

121 ) 



H7=4 



(N 

122) 



N7=l 



(N 

123) 



XLQADC 1 ) =1 .0 



CN 

124) 



IF (ETAl.NE.l.) GO TO 812 



CN 

125) 



M7=2 



( N 

126) 



$UML=0 .50 



CN 

127) 



GO TO 812 



CN 

128) 


614 

SUML=0. 126104 



CN 

129) 



IF CLIW.EQ.i) SUML= 0. ICO 



CN 

1 30) 



N7=10 



CN 

131) 



H7=l 



CN 

132) 



IF ( ETAl .NE • 1 • ) GO TO 812 



(N 

133) 



SUML=0. 25208 



(N 

134) 



IF CLIW.EQ.il SUMI=0. 200 



CN 

135) 



N7= 5 



(N 

136) 


612 

00 801 Ml=l i M7 



CN 

137) 



CO 802 Nl = l ?N7 



CN 

138) 



XST AR= ! 11 .-2 .*FLOATC Ml ) ) /10 . 



CN 

139 ) 



YST AR-t ll.-2.*FL0ATCNl> 1/10. 



CN 

140) 



GO TO C 621t 622 ), I EL EM 



(N 

141 I 

c 





CN 

142 ) 

c 


EFFECT OF ROTOR 



(N 

143) 

c 





CN 

144) 


621 

ETA=ETA1-RL3ADCN1 )*SIGMAR*SINCPSI CN1) ) 



CN 

145) 



Z£TA=ZETA1/(1 ,0-RLO AD(N 1 )*S I GM AR*GAMMA*ZETA1 *COS (PS I (Nl))* 

CN 

1 46) 



1 SIN( ALPHAR ) ) 



(N 

147) 



XOVERH=SIGMAR*GAMMA*(TL*COS( ALPHAS »4-TH*S IN ( ALPHAS) -RLO A0( Ml ) * 

CN 

148) 



1 CCS( ALPHAR)*COS(PSI(Nim 



CN 

149) 



YOV ERH=S IGMAR*GAMMA*(0.25*(5.0-2.0*FL3AT(Ml ) )*( S IGMAT/SI GMAR ) - 

CN 

150) 



1 RL0AD(N1)*SIN(PS KN1 )>) 



CN 

151) 



ZOVERH=S IGMAR*GAMMA*(TH*COS( ALPHAB ) -TL*S I N( ALPHAB) +RLDAD* 

CN 

152 ) 



1 SIN(ALPHARl*COS(PSI(Nim 



(N 

153) 



XLOAOC Nl 1 = 1.0 



CN 

154) 



GO TO 630 



(N 

155) 

c 





CN 

156 I 

c 


EFFECT OF WING 



CN 

157) 

c 





(N 

158) 


622 

ZETA=Z ETA1/ ( 1 .0-S IGMAR*GAMMA*ZET Al* ( ASS (YST AR >*( SIGMAW/SIGMAR) 

CN 

159 I 
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1 *T AN ( LAMBDA I *S IN ( ALPHAS I +LW S IN ( AL PHAB) -HW*COS( ALPHAS ) ) ) 

ET4=ETA1-YSTAR*SIGMAW 

XOVERH=S IGM AR*GAMMA* l (TL-LW)*COS( ALPHAS )♦( TH-HWI *SIN( ALPHAS ) - 
1 ABS(Y ST AR)*(SIGM AW/SI GM4RI*TAN( LAMBDA I *COS ( ALPH AB) I 

YOVERH=GAMMA*( 0 .25 *( 5 .0-2. 3*FLO AT (Ml ) > *S IGM AT-YSTAR*S IGM AW ) 
ZGVERH=-S IGM AR *GAMMA*( ( TL-LW ) *S IN I ALPH AB) ♦ ( TH-HW )*COS( ALPHAS) - 
1 ABS( YSTAR)*( SIGM AW/S IGM AR)*TAN( LAMBDA) *S IN (ALPHAS) ) 

630 CALL OLTAS (C(K)I 

* 4 * 4 * **************** SEE APPENDIX 0 FOR SUBROUTINE DLTAS *********** 
DO 805 Ll = l » 28 

805 XOELTA( LI )=XOELTA( LI )+DELTA(Ll )*XLOAD( N1 ) 

802 CONTINUE 
801 CONTINUE 

DO 807 L3=l,28 

807 DELTA! L3 )=XDELTA(L3 )*SUML 
WRITE (6,149) C ( K ) 

(OELTA(I), 1=1,25, 4) 

=2,26, A) 


WRITE 
WRITE 
WRITE 
WRITE 
DG 814 


(6,150 ) 
(6,151 ) 

(6.152) 

(6.153) 

L 4=1 » 28 


( DEL TA ( I ) , I’ 

(DELTA( I ) , 1=3,27,4) 
( DELTA( I ) , 1=4,28,4) 


814 XDELTAI L4 ) = 0. 

41 CONTINUE 

42 CONTINUE 
GO TO 1 

100 FORMAT (1H1/////////59X*START NEW CASE* ) 

149 FORMAT (1X*CHI =* F7.3/) 

150 FORMAT ( 3X5H( W , L )7( F17.4 ) ) 

151 FORMAT (3X5H(U,L )7( F17.4 )) 

152 FORMAT ( 3X5H( W, D) 7( F 17.4 ) ) 

153 FORMAT ( 3X 5H( U , D )7( F 1 7 .4 )// ) 

210 FORMAT (1X13K1F-) ) 

211 FORMAT (1X1HI11X1H131 X6 1HC3RREC T I ON FACTORS FOR CORRECTING FROM A 
1WIND TUNNEL WHICH IS25X1HI) 

212 FORMAT (1X1HI 11X 1H 1117(1 1 — ) 1 H I ) 

213 FORMAT ( 1X1HI11XIHI16X1HI5X6HCLOSED5X1HI 16X1HI 2X12HCLOSED FL00R2X1 
IF I6X4H0PEN6X1H 1 16X1H 1 5 X6HCL0SED4X 1H I ) 

214 FORMAT (1X1HI3X5FDELTA3X1HI 5X6HCL0SED5 XI HI 4X9H0N BOTTOMS XI H I 6X4 HOP 
1EN6X1H I 6X4H0NL Y6X1H 1 5X5HFL00R6X 1H I 5X6HCL0SE05XIH 13 X9HQN B0TT0M3X1H 
21) 

215 FORMAT ( 1X1HI11X1HI1 6X1H I6X4H0NL Y6X1HI 16X18HI (GRCUND EFFECT) I6X4H 
1GNLY6X1HI16X1H 16 X4H0NL Y 5X1H I > 

216 FORMAT ( 1X1H II 1 X1H I 84( 1H-) 1H 13 2( 1H- ) 1HI ) 

217 FORMAT ( 1X1HI11X1HI 36X11HT0 FREE AIR37X1HI8X16HT0 GROUNO EFFECT 8X1 
1HI) 

218 FORMAT ( 1X131 ( 1 H- )/ ) 

701 FORMAT (1H1//35X*AVERAGE INTERFERENCE OVER A TAIL BEHIND AN UNLOAD 
1E0 ROTOR MODEL*//) 

702 FORMAT (55X*EFFECT OF RCTOR CN TAIL*//) 

703 FORMAT (55X*EFFECT OF WING ON TAIL*//) 

706 FORMAT (34X.A10* ROTOR L0ADING*19X , A8* WING LOADING*// 

119X*SI GM A ( ROTOR ) =*F6.3 , 10X*2E TA =*F6. 3 , 10X*LW/R =*F6.3,10X 
2*ALPHA( ROTOR) =*F7. 3// 1 9X*S IGMA ( WING ) =*FS . 3 , 10 X*ET A =*F6.3, 

2 10X*HW/R =*F6.3,10X*ALPHA(B0DY) =*F7.3//19X*SIGMA( TAIL ) =* 

4F6.3,21X*TL/R =*F6.3,30X*TH/R =*F6. 3//39X*GAMMA =*F6.3,27X 
5*LAMB0A =*F7 .3// I 

707 FORMAT (40X*SIGMA( ROTOR) EQUALS ZERO, THIS PROGRAM IS NOT SUITABLE 
1 FOR USE WITH SUCH CASES.*//) 

900 FORMAT ( I1,F9.3,4F10.3/I 1,F9.3,4F 10. 3/3F10. 3) 

999 STOP 
END 


I N 

160) 

IN 

161 ) 

IN 

16?) 

IN 

163) 

(N 

164) 

IN 

165) 

IN 

166) 

( N 

16?) 

E******* 

IN 

168) 

(N 

169) 

IN 

170) 

IN 

171) 

(N 

172 ) 

(N 

173) 

(N 

174 ) 

IN 

175) 

IN 

176) 

(N 

177) 

IN 

178) 

(N 

179) 

IN 

180) 

IN 

181 ) 

CN 

182) 

IN 

183) 

(N 

184) 

IN 

185) 

IN 

186) 

(N 

187) 

(N 

188 ) 

IN 

189 ) 

(N 

190 ) 

IN 

191 ) 

IN 

192) 

IN 

193) 

IN 

1 94) 

IN 

195) 

IN 

196) 

IN 

197) 

IN 

198) 

IN 

199) 

IN 

200 ) 

IN 

201 ) 

IN 

202 ) 

IN 

203) 

IN 

204) 

IN 

205) 

IN 

206) 

CN 

207) 

(N 

208) 

IN 

209 ) 

IN 

210 ) 

IN 

211 ) 

IN 

212 ) 

IN 

213) 

IN 

214) 

IN 

215) 

IN 

216) 

IN 

217) 

IN 

218) 

IN 

219) 
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APPENDIX 0 


FORTRAN PROGRAM FOR CALCULATING THE AVERAGE WIND-TUNNEL 
INTERFERENCE OVER SIOE-BY-SIDE R0T3R CONFIGURATIONS 


THIS PROGRAM WAS WRITTEN IN COC FORTRAN, VERSION 2.1, TO RUN ON CDC 5000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENT IONE 0 COMPUTERS WHICH CARRY 
TFE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE OLTAS WHICH IS GIVEN IN 
APPENCIX Q. 


INPUT WILL BE FOUND AT ACCRESS 1 ITHREE CARDS PER CASE) IN FORMAT 930. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN AT THE POINT MIDWAY BETWEEN THE TWO ROTORS. 
THE REQUIRED INPUT VARIABLES FOR THE ROTORS , GIVEN ON THE FIRST CARD, ARE 


LIR ROTOR D IS K- LOAD -01 STRI8UTIQN INDICATOR, LIR=1 FOR UNIFORM L3 AD- 

ING, L IR=2 FOR TRIANGULAR LOADING 

ZETA1 SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT CF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FRCM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMAR RATIO OF ROTOR CIAMETER TO TUNNEL WIDTH 

SR DISTANCE BETWEEN ROTOR CENTERS DIVIDED BY ROTOR DIAMETER 

ALPHAR ANGLE OF ATTACK OF TIP-PATH PLANE OF ROTORS, DEG 

THE REQUIRED INPUT VARIABLES FCR THE WING, GIVEN ON THE SECOND CARO, ARE 


LIW WING SPAN-LOAD-OISTRIBUTION INDICATOR, LIW=1 FOR UNIFORM LOADING, 

L IW=2 FOR ELLIPTIC LOADING 

LW DISTANCE OF WING APEX BEHIND ORIGIN, NO NDIMENSIONALI ZED WITH RE- 

SPECT TO ROTOR RADIUS 

HW OISTANCE OF WING APEX ABOVE ORIGIN, NONDI MENSI ON AL I ZED WITH RE- 

SPECT TO ROTOR RADIUS 

SIGMAW RATIO OF WING SPAN TO TUNNEL WIDTH 
LAMBDA WING SWEEP ANGLE, DEG 
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AUMMil ANGLE OF ATT ACK OF BODY CARRYING WING AND TAIL, DiG 


THIS PROGRAM REJECTS CASES OF ZFRO SPAN. SINCE THE EQUATIONS ARC FORMED 
IN TERMS JF ROTOR RADIUS, SUCH CASES REPRESENT INPUT ERRORS. THIS PROGRAM ALSC 
DETERMINES AND REJECTS CASES IN WHICH TOTAL ROTOR SPAN EXCEEDS THE TUNNEL WIOTH. 


THIS PROGRAM COMPUTES INDEPENDENTLY THE INTERFERENCE AT EACH OF THE THREE 
ELEMENTS DUE TO ITS OWN PRESENCE , AS WELL AS THE INTERFERENCE ON EACH DUE TO TEE 
PRESENCE OF THE OTHER TWO ELEMENTS. IN SYMMETRICAL CASES, THE EFFECTS ON, AND 
CAUSED BY, EACH OF THE ROTORS IS IDENTICAL. CONSEQUENTLY, ONLT THE INTERFERENCES 
RELATED TO ONE OF THE ROTORS IS CALCULATED. 


NOTE THAT THIS PROGRAM IS ALSO SUITABLE FOR TW I N- PROPF L LER TILT-WING MOD- 
ELS, TILT-ROTOR MODELS, FAN-IN-WING MODELS, AND OTHER SIMILAR TYPES. JUDICIOUS 
CHOICE OF INPUT VARIABLES WILL SATISFY THE REQUIREMENTS OF THESE, AND MANY 
OTHER MODELS. 


PROGRAM WINDTUNUNPUT, OUTPUT, TA PE 5 = I NPUT , T APF6 =OUT PUT ) (0 1) 

COMMON ZET A , 5 T A , GA mm A , XO VERH , YO VER H, Z OVE RH , DE L T A ( 2 8 ) ( 0 2) 

DIMENSION X DE L T A ( 2 8 I , P S I ( 2 P I , R L 0 A D ( 20 ) , R UN I F ( 2 0 1 , R TR I A ( 20 ) » (0 3 ) 

I XLE( l n ) »XLDA0(20),C( 9) (0 4) 

REAL LAMBDA, LW ( C 51 

DATA (RUNIFI I) , 1 = 1 ,201/4*0. 2981 ,8*0.6255,8*0.8921/ ( 0 6) 

DATA (RTRIA(I)»I=1,?0)/4*C'.4386, 8*0.72 96,R*C , .9?6 2/ (C7I 

DATA (C( I) , 1 = 1 , 8)/20. ,30. ,43. ,53. ,63. , 7D., 80. , 90. / (C8) 

PI=3 .14159268358970 (0 9) 

R A 0=0. 31 74 5 32 92 SI 99 (0 10) 

DO 303 LL=L,28 (Oil) 

803 XDELT4IL1 > =">. (0 121 

PS If 1 )= (P 1/4. ) (013) 

PS I( 2 > = 3 .*PS I ( 1 ) (0 14) 

PSH 3 ) = 5 . * P SI ( 1) (D 15) 

PS l(4) = 7.*PS(< 1 ) (0 16) 

PS I(S) = PSI (13 ) = ( P I / 8 . ) (0 17) 

PS I( 6) = PSI (14)=3.*PSf( E) (0 18) 

PSI(7)=PS 1(1 ? >=5.*PSI (5) (0 19) 

PS I< 8) =PSI ( 16) =7.*PSI t 5) (0 20) 

PS I( 9 )=PS I I 1 7 ) = 9. *PSI ( 5 ) (0 21) 

PS I l l))=PS I( 18 )=1 1 .*PSI (S ) ( 0 22) 

PS I( 11 ) = PSI ( 19 )= 13.*PSI ( S) (0 231 

PS I( 12 ) = PS 1(23 I = 1 5 . * P S I (5) (0 24) 

XLE( 1)=XLE( IP 1=0.43579 (0 25) 

XL E( 2 I = XLE ( 9 ) =0.71422 ( 0 26) 

XL E( 3)=XLE< 8) =*. 86633 (0 27) 

XL E( 4 l = XL E ( 7 ) =0. 95304 ( 0 28) 

X L E ( 5 ) = XLE ( 6 ) =0.99499 (0 29) 

1 READ (5,90*)) LIR, 7.ETA1,ETA1, GAMMA, SIGMAR, SR, AL PH AR , L I W , L W , HW , (0 30) 

1 SIGMAW, LAMBDA, ALPHA!) (0 31) 

I F ( EOF, 5 I 099,47 ( 0 32 ) 

4 7 A ALP R= AL PHAR ( 0 33) 

AALP3=ALPHA3 (0 34) 

AL AM = L AMBDA (0 35) 

LAMBDA = L AM90A*RAD (0 36) 

ALPHAR= ALPH4R*RA0 (0 37) 

ALPHA3=ALPHAR*RAD (0 38) 

WRITE ( 6 , 100 ) (0 39 ) 

IF ( L IR .EG .1 ) GO TP R04 (0 40) 
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I ALPHA=10HTRI ANGUL AR 

00 80 B M?=U20 

80 8 RL 0 AD ( M 2 ) =RT R 1 A( M2 ) 

GO TO 806 

80 4 I AL 3 HA= 10HUN I FORM 
DO 8 09 M2 = L ,20 
800 RLOAO( M2 >=RUNIF (M2) 

806 IF (LIW.EQ.I ) GO TO 35? 

1 f i ETA = 3HEL L I P T I C 
DO 851 M 3= 1,10 

851 X L 0 A 0 ( M 3 ) = X L F ( M3 ) 

GO TO 350 

852 I 6 tT A = 8 HUM t FOR M 
DO 853 M 3 = 1 ,10 

P53 XLOAOt M3 )=1 .o 
8 5C WIDTH-S IGMAR*( l.o+SR) 

IF (WIDTH. IT. 1 .. AND.SIGMAP .NE.O.) GO TO 855 
m ITE (6,701 ) 

WR ITE (6,711) I ALPHA, I 3 E T A , S I QM AR , Z £T A1 , LW , A A L PR , S IGM A W , E T A1 , 

1 HW, A ALP 3, GAMMA f SR , ALA M 

WRITE ( 6 » 2 1 n ) 

WR IT E (6,211) 

WRITE (6,212) 

WR ITE (6,213 ) 

WRITE (6,214) 

WRITE (6,215) 

WR ITE ( 6,216 ) 

WRITE (6,217) 

WR ITE (6,213 ) 

IF ( S I GMAR • ED # n • ) WRITE (6,712) 

IF (WIDTH. GE.l .) WRITE (6,717) 

GO TO 1 

85 5 DO 4 2 I FL EM =1 ,9 

IF (ETAl.EQ.l.) GO TO ( 6P0 f *nf), 6^0 , 42 , 42 , 42 , 60* f 6C-'3 ,42 I , IELEM 
60 C WR ITE (6,701 ) 

IF ( ETA1.EQ.1. ) GO TO ( 5 n 1, 50 2, 50 3, 42 , 42 *42 ,6^7 ,504,42 ) , IELEM 
GO TO (601,6- 2,603, 6° 4 f 6C5, 6 06,607, 6^8 ,6^9) , IELEM 
6^ l WR ITE ( 6, 7 32 ) 

GO TO 610 

60 2 WR I TE ( 6, 703) 

GO TO 610 

603 WRITE ( 6,704) 

GO TO 610 

604 WRITE (6,705 ) 

GO TO 510 

60 5 WR ITE (6,705 ) 

GO TO 6 10 

606 WRITE (6,707) 

GO TO 610 

607 WRITE (6,7 >8 ) 

GO TO 610 

50 8 WRITE (6, 709) 

GO TO 6 1 9 
609 WRITE ( 6,710 ) 

GO TO 610 
50 1 WR ITE ( 6,713 ) 

GO TO 619 

502 WRITE (5,714) 

GO TO 610 

503 WRITE (6,715 ) 

GO TO 610 


(0 41 ) 

( 0 42 ) 
(0 43 ) 
(0 44 ) 
(0 45 ) 
(0 46 ) 
10 47 ) 
(0 48 ) 
(0 49 ) 
(0 50 ) 
(0 51 ) 
(0 5 ?) 
( 0 53 ) 
(0 54 ) 
( 0 55 ) 
( 0 56 ) 
(C 57 ) 
( 0 58 ) 
(0 59 ) 
(0 6 ° ) 
( 061 ) 
( 0 62 ) 
(0 63 ) 
(0 64 ) 
(0 65 ) 
(0 66 ) 
( 0 67 ) 
( 0 68 > 
( 0 69 ) 
( 0 70 ) 

(0 71 ) 
( 0 72 ) 
(0 73 ) 
(0 74 ) 
( 0 75 ) 
(0 76 ) 
(0 77 ) 
(0 78 ) 
(0 79 ) 

(0 39 ) 

(0 31 ) 
(0 82 ) 
(0 R 3 ) 
(0 84 ) 
(0 85 ) 
( 0 8 6 ) 
(0 37 ) 
(0 38 ) 
(0 89 ) 
(0 90 ) 
(0 91 ) 

(0 92 ) 
(0 9 3 ) 
(0 94 ) 
(0 95 ) 
(O 96 ) 
(0 97 ) 
(0 98 ) 
(0 99 ) 
(0 100 ) 
(0 101 ) 
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5^4 WRITE ( 6,71b ) 

610 WRITE (6,711) I M PH A , I BE T A , S I GMAR , 7 R T A i , L W , A AL P R , S I GM AW , ET A l , 
1 HW, AALPO, GAMMA, SR , AL AM 

WR ITS (6,21" ) 

WR IT F (6,211) 

WRITE (6,212 ) 

WRITE (6,213) 

WRITE (6,214) 

WRIT E (6,215 ) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,213) 

DO 41 <=1,3 

M7 =N 7 = 2 0 

GO TO (611,611,612,611,611,612,613,614,614), IELEM 
ROT'JR ON ROTOR 

611 SU ML =0.002 5 
GO TO 812 

WING ON ROTOR 

612 SO ML = 0 • 6 30 5 7 

IF ( SIGMAW .EQ.O. ) GO TO 6 1^ 

IF (LIW.EQ.i) SUML=0 .O n 5 
N7 =1 0 
GO TO 312 

615 S’J ML =0 . ° 5 
XLOAD( 1 ) = l . 7 
N 7 =1 

GO TO 812 

WING ON WING 

61? S U ML = 0 * 0 1 2 6 1 0 4 

IF ( SIGMAW .EQ.O. ) GO TO 6 16 
IF (LIW.EQ.l ) SUML = r 1° 

M7 =10 
N7 =1 0 

IF (ETA1.N6.1 . ) GO TO 312 
SO ML =0 •'■'2 5220 8 
IF ( LIW.EQ.l ) SUML = A • G 20 
N 7 =5 

GO TO 312 

616 SO ML = 1 « 

XL OAD( l ) = 1 
M 7 =N 7 = 1 
GO TO 812 

ROTOR ON WING 

614 SU ML =3. "0 59 
M 7 = 1 3 

IF ( SIGMAW.NF.^.I GO TO 812 
M7 =1 

S U ML = 3 .05 

81 2 00 831 *1=1 , M 7 

DO 302 Nl =1 , N 7 
XS TAR=( 11.-2 .AFLOAT ( Ml ) ) /IP. 

YSTAR= ( 11. -?.*FL0AT(N1) )/10. 


(0 10 2) 
(0 103) 
(0 1^4) 
(0 1° 5 ) 
(0 1 0 6 ) 
(0 107) 
(0 108) 
(o io<n 
(o no 
(0 111 ) 
(0 112) 
(0 113) 
(0 114) 
(0 115) 
(0 116) 
(0 117) 
(0 118) 
<0 119) 
(0 120) 
(0 121 ) 
(0 12?) 
(0 123) 
(0 124) 
(0 125) 
( 0 126) 
(0 127) 
(0 128) 
(0 12 0 ) 
(0 130) 

( 0 131 ) 
(0 132) 
(0 138) 
(0 134) 
(0 135) 
(0 136) 
(0 137) 
(0 133) 
(0 1 7 9 ) 
(0 140) 
(0 141) 
(0 14? ) 
(0 143) 
(0 144) 
(0 145) 
(0 146) 
(0 147) 
(0 148) 
(0 149) 
(0 15^) 
(0 151) 
(0 152) 
(0 153) 
(0 154) 

(0 155) 
(0 156) 
(0 157) 
( 0 158) 
(0 159) 
(0 16^) 
(0 161) 
(0 16?) 
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Appendix O - Continued 


811 GU TO (621 ,622 ,6?3, 624, 625, 626, 627, 628, 629) , IELEM 

RIGHT R0T3R IN RIGHT ROTOR 

621 ET A=ET41 -RL0ADIN1 )*SI GMAR*SIN( »SI( Nl )» -SR*S I GM AR 

YOVEI<H=S IG8AR*GAMMA* (RLOAOI Ml) *S IN ( PS I ( 81 >) -RLO AO INI ) *S18(P SI 
1 (Nl) ) ) 

GO TO 631 

LEFT ROTOR ON RIGHT ROTOR 

62 2 ET A= ET A 1 -RLOAD I Nl )$SIGMAR*SIN( PSI (Nl) )+SR*SlGMAR 

Y0VERH=SIG8AR*GAMMA*(R10AD(M1)*S IN I PS I ( 81 ) ) -RL 0 AO ( N 1 ) * S l Ml P S I 
l (Nl) J«-2.0*SR) 

GO TO 631 

LEFT ROTOR IN LEFT ROTOR 

62 6 ET A= ET A1 -R LOAD (Nl )*SIGMAR*SIN< PS I( Nl ) )+SR*S!GMAR 

Y3VERH=S IGMAR*GAMMA*( RLOAOIM1J *S IN (PSI (81 ) ) -RLO A OI N 1 ) *S I N I P S I 
1 (Nl) ) ) 

GO TO 631 

RIGHT ROTOR ON LEFT ROTOR 

62 5 ET A=ETA1-RLUADIN1 ) *S I GMAR*=SI N( PSI ( Nl ) )-SR*S IGMAR 

YO VfcRH= S IG MAR *GAMMA*( RIO AO ( M 1 ) *S I N ( PS I ( 81 ) ) -RL 3 A D( N 1 ) * S IN ( P S I 
1 (Nl ) )-2 .)*SR ) 

631 ZETA=7ETA1 /(l.O-RLOAD(Nl) *S I GM AR*G A MM A*ZE T A1 *COS ( P SI ( Nl ) ) * 

1 SIN(ALPHAR)) 

XO VERH= SI G8AR* GAMMA *COS( ALPHAR)* (RL)AO(Ml )* COS (PS 1(81 ) )- 
1 RLOAOINl )*COS ( PS I (Nl ) l ) 

ZOVtRH=-S IGMAR *GAMMA*SIN( AL PHAR )*( R LOATH Ml I *COS ( P SI (M l ) ) - 
1 RL 0 A 0 ( N 1 ) *CO S t P S I ( N 1 ) ) ) 

XLOAOI Nl ) =1 .0 
GO TO 640 

WING ON RIGHT ROTOR 

62 3 YO VERH=SI3M4R*CA8MA* (RL04D(M1)*SIN(PSI(81))*-SR-YSTAR* 

1 ( S IGM4W/S IGMAR ) ) 

GO TO 632 

WING ON LEFT ROTOR 

626 YOVERH=S IG8AR *GAMMA*( RLOAD(Ml) *S IN(PS I (81 ) )-SR-YSTAR* 

1 ( SIG.84W/S IGMAR )) 

632 ET 4=ETA1-YSTAR*S!GMAW 

ZET4=ZETAl /( 1 . C-S l GM AR*G AMM A* Z ET 41 * ( ASSIYSTAR )*( SI GM A WZ S I GM AR )* 
1 T AN ( L AMrtO A ) * S IN ( ALPHAS ) +LW* SI N I ALPHAS ) -HW*COS (ALPHAS) ) ) 

ZOVERH=-SIGMAP *GAMMA*(RLO AM Ml ) *S IN( AL PH AR ) *C0 S( P S I( M 1 ) )- 

1 ASS (YS TAR )*( SIGMA W/SI GMAR ) *T AN ( LA MBO A ) *S I N ( AL PH AS) - 

2 LW*S IN ( ALPHAS ) 4-HW*C0S( ALPHAS )) 

XO VERH=S IGMAR ♦GAMMA*! RLO AD(M1)*C0S I ALPHAS )*COS(P St ( Ml ) )- 

1 ASSIYSTAR ) *( S IG8AW/SIGMAR ) *TAN(LA8B0A)*CDS( ALPHAS ) - 

2 LW*COS ( AL P F A P )-HW*S IN( ALPHAS ) ) 

GO TO 640 

WING ON WING 

627 ET A=ETA1-YSTAR*SIGMAW 


(0 

163) 

<n 

16 4) 

( 0 

165) 

( 0 

166) 

(0 

167) 

i 0 

16 R I 

(0 

169) 

(0 

170) 

( 0 

171 ) 

(0 

172) 

i 0 

173) 

<Q 

174) 

( 0 

175) 

CO 

176 ) 

(0 

177) 

(0 

178) 

( 0 

179) 

<0 

180) 

(0 

181 ) 

(0 

182) 

<0 

18 3) 

(0 

184) 

(0 

185) 

(0 

186) 

( n 

187) 

n 

188 ) 

( G 

189) 

( o 

190 ) 

(G 

191 ) 

( G 

192) 

(0 

193) 

(0 

194) 

( n 

195 ) 

<n 

1%) 

(0 

197) 

(0 

198 ) 

( n 

199) 

( G 

2V) 

(G 

2 r 1 ) 

to 

2C2) 

( 0 

203) 

( G 

2^4) 

(0 

205) 

( 0 

2^6) 

(G 

20 7 ) 

(0 

2^8) 

( n 

20 9) 

(0 

210) 

( 0 

211 ) 

(□ 

212) 

i 0 

213) 

in 

214 ) 

(0 

215) 

( 0 

216) 

(0 

217) 

(G 

218) 

(G 

219) 

(G 

22 n ) 

(G 

221 ) 

in 

222 ) 

(0 

223 ) 
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ZETA = ZETAl/( l .'•-SIGMAR*GAMMA*ZETAI*IYSTAR*( S I GM AN/ S l GM AR ) * 

1 TAN ( LAMBDA ) *S IN ( 4L PHAR ) + LW*S IN ( ALPHAS > -HW*C()S< ALPHAB) ) ) 
XOVERH=SIGMAW*GAMMA*T AN (L AMBDA)*Cf)S( ALPHAS)*! ABSIXSTAR)- 
l AOSIYSFAT ) > 

YO VERH=l ,2*SIG«AW*GAMMA*( FLOAT! Nl) -FLOAT (Ml) ) 
ZnV£RH=-SIG'AAWX‘GAA ,M A*T AN (l AMBOA )* S IN I AL PH AS )* I AB S I X ST AR ) - 
l A3S I YS T AR ) ) 

GO TO 6 40 

RIGHT ROTOR ON WING 

62 8 ET A=ET A1 -R LO AD *S I CM AR *S IN (PS I INI ) >-SR*S IGMAR 

Y0VERH=SIGMAR*GAMMA*(XSTAR*(SI3MAW/STGMAR I-RLOAOINl )*SINI PSI 
1 IN1U-SR) 

GO TO STB 

LEFT ROTOR ON WING 

62 0 ET 4=ET AI-RLO AD*S I CM AR *S IN I PS 1 1 N1 ) )+SR *S l G'* AR 

YOVERH=SlSMAR*GAMMA*(XSTAR*(SI GMAW / S l G MA R ) -RL 0 AD (N l ) * S I N ( PS I 
1 (NI))+SR) 

63 3 ZETA=ZETAI/( 1 . O-SI GM AR*G AMMA* Z E T A I *RLOAD ( NL ) *S I N ( AlPHAR)* 

1 COSIPS I INI ) ) ) 

XOVEKH=SIG‘HAR*GAM M A*( ABS(XSTAR)*(SIGMAN/SIGMAR)*TAN(LAMBDA)* 

1 COS I AL PH AB ) + 1 W*COS ( AL P HAB ) > HW * S I N I A L PH AB ) -PLO ADI Nl ) * 

2 COSI ALPHAR )*CCS(PSI INI m 

ZOVERH=-SIGMAP*GAMMA* (ABS(XSTAR)*! SIGMAW/SIGMAR) *T AN ( L AMBOA I * 

1 S INI ALPHAB )+LW*S IN 1AL PHAB)-HW*C )S I A L PH AB ) -RL D ADI N 1 ) * 

2 SINI ALPHAR )*CPS(PSI(N1 >) > 

XL DAD INI )=1.H 

64 n CALL DLTAS ( C IK ) ) 


(0 

22 4 ) 

(0 

225) 

( a 

22 6) 

(0 

227) 

(0 

22 8) 

(G 

229) 

( □ 

2 39 ) 

(0 

231) 

(G 

232) 

( 0 

233) 

CO 

234) 

{ G 

2 3 5) 

(0 

236 ) 

( 0 

237) 

( 0 

238 ) 

(0 

239) 

(0 

24 r ) 

(0 

241 ) 

<0 

242) 

( G 

243) 

(0 

244) 

(0 

245 ) 

( 0 

246) 

( 0 

247) 

(0 

248 ) 

(0 

249) 

(0 

250 ) 

( 0 

251 ) 

(0 

25?) 

(0 

253) 

( G 

254) 


**** **** ************* 4PPFN0IX Q FOR SUBROUTINE DLTAS ********************** 



OG 305 

L 1 = 1 *23 

(0 

255) 

80 5 

X CELT A I LI ) =XOELTA (LI ) +DEL T A I L 1 ) *X LO ADI Nil 

(0 

256 ) 

80 2 

CONT INUE 

(0 

257) 

30 1 

CONTINUE 

(0 

258) 


on 80 7 

13-1*28 

(0 

259 ) 

30 7 

0 E L T MLB) = X DEL T A U 3 ) *SUML 

(0 

260 


WR ITE 

(*6,149) C(K) 

(0 

261 ) 


WR ITE 

< 6 * 1 5 r * ) t DELTM I ) , 1 = 1 , 25 ,4) 

(0 

262 ) 


WR ITE 

(6*151) < DE LT A ( I ) ,1=2*26,4) 

(0 

26 3) 


WR ITE 

(6,15?) { DELTA( I) , 1=3* ?7 f 4) 

( 0 

264) 


WRITE 

(6,153 ) (DELTA! I ) , 1=4,23,4 ) 

(0 

265) 


00 814 

L4= 1 » 28 

(0 

266) 

PI 4 

XDELTM L4 ) =0 * 

( 0 

267) 

41 

CONTINUE 

(0 

268 ) 

42 

CONTINUE 

(0 

269) 


GO TO 

1 

( 0 

270) 

100 

FORMAT 

(1H1/////////59X*START NEW CASE*) 

CO 

271 ) 

149 

FORMAT 

( 1 X *C H I = * F7.3/) 

( 0 

272) 

1 5 C 

FORM AT 

( 3X 5 H( W f L ) 7 ( FI 7.4 ) ) 

(0 

273 ) 

151 

FORMAT 

(3X5H(U»L ) 7 ( F 1 7 • 4 ) ) 

<0 

274) 

15? 

FORMAT 

(3X5H(W,D)7(F17.4) ) 

CO 

275 ) 

153 

FORM AT 

(3X5H(U,D)7(F17.4)//) 

(0 

276 ) 

210 

FORMAT 

(1X131 C1H-) ) 

( 0 

277) 

21 1 

FORMAT 

( IX 1 HI 1 IX1HI31 X61HC0RRECTinN FACTORS FOR COPRCCTING FROM A 

(0 

27 8 ) 

1WIN0 TJNNEl WHICH IS 2 5X1 H l ) 

(0 

279) 

21 2 

FORMAT 

I 1X1HI 11XIHU17I 1H-) IHI ) 

(0 

280 ) 


213 FORMAT I1X1HI11X1HII6X1HI5X6HCL0SEDSX1HI16X1HI2X12HCL0SED FLOORRXl (0 231) 

IHI 6X4H0PEN6XIHI 16X1HI 5X 6HCI.0 SE04 XI H I ) (0 282) 

214 FORMAT ( 1 X 1 H I 3* 5H0EL T A3 X1H I 5 X 6HCL OSED 6X 1H I 4X 9H0N B CT TOM 3 XU I 6 X4 HOP (0 283) 
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l E N 6X 1 H I 6 X4 HONL Y6 X 1 H I 5X 5 HFL00R6 X 1 H I 5X6HCL 3 S F 05 X H I 3X9H0N BOTTOM 3 XIH ( 0 


21 ) (0 

215 FGR'IM ( 1X1HM IX 1 H I 1 6X1HI 6X4H0NL Y6 X 1HI 16X 1 8HI ( GROUND EFFECT) I6X4M (0 

10NLY6X1HI 16X1HI 6X4H0NL Y5X1HI ) (0 

216 E0RM4T (lXIHm X1HI«4( 1H-UH132C 1H-)1HI) (0 

217 FORMAT < 1 X 1 H 1 1 1 X 1 H I 3 6X 1 1 HTC1 FREE A I R37X 1H I «X 1 6HTO GROUND EE EEC T 8 X 1 (0 

1 H I ) (0 

218 FORMAT (1X131(1 H-)/) (Q 

7C1 FORMAT ( 1H1//35X*AVFRAGE INTERFERENCE OVER SIDE-BY-SIDE AND/OR TIL (0 

IT-ROTOR MODEL S */ ) (0 

702 FORMAT (48X*EFFECT OF RIGHT ROTOR ON RIGHT ROTOR*//) (0 

703 FORMAT <43X*EFFECT OF LEFT ROTOR ON RIGMT ROTOR*//) (0 

7 04 FORMAT <5lX*EFEECT OF WING ON RIGHT ROTOR*//) (0 

705 FORMAT (49X*EEFECT OF LEFT ROTOR ON LEFT ROTOR*//) (0 

70 6 FORMAT (4SX*EPFECT OF RIGHT ROTOR ON LFFT ROTOR*//) (0 

70 7 FORMAT ( 52 X *t E F ECT OF WING ON LEFT ROTOR*//) (0 

708 FORMAT (55X*EFFECT OF WING ON WING*//) (0 

709 FORMAT (51X*EFFECT OF RIGHT ROTOR ON WING*//) (0 

71 C FORMAT <52X*EFFECT OF LEFT ROTOR ON WING*//) (0 

711 FORMAT ( 34 X * A1 0* ROTOR LO AD I N G*1 9 X , A R * W I NG LOADING*// (0 

119X*SIGMA( ROTOR) = *F6 . 7 , IP X*ZF TA = *F6* 7, IPX* LW/R =*F6.3,1PX (0 

2*ALPHA( ROT OR ) =*F7 . 3/ / 1 ^X*S I GM A ( WING) =*F 6 . 3 , 1 ^ X* E TA =*F6.3 , ( Q 

310X*HW/R = *F 6 • 3 1 1 ° X* Al PH A ( BODY ) = *F 7 . 3 / / 1 Q X*G AM M A =*F6.3, (0 

4 28 X* S R/R =*F6.3,27X*t AMBDA =*F7*3//> (0 

712 FORMAT ( 40 X* S I GM A < ROTOR ) EQUALS ZERO, TMIS PROGRAM IS NOT SUITABLE (0 

1 FOR USE WITH SUCH CASES**// ) (0 

713 FORMAT (53X*FFFECT OF EITHER ROTOR ON ITSELF*//) (0 

714 FORMAT (47X*EFFECT OF OTHER ROTOR ON EITHER ROTOR*//) (0 

715 FORMAT (51X*EFFECT OF WING ON EITHER ROTOR*//) (0 

716 FORMAT (51X*EEFECT OF EITHER ROTOR ON WING*//) (0 

717 FORMAT (4DX*R0T0R SYSTEM IS TOO WIDE FOR THE WIND TUNNEL*//) (0 

900 FORMAT ( 1 1 , F9 . 3 , 5E 1 r \. 3 / I 1 , E 9 . 3 , 4F IP . 7 ) (0 

999 STOP (0 

END (0 


28 4) 

285 ) 
28 6) 
287 ) 
2 8 8 ) 

289) 

290) 

291 ) 

292 ) 

293) 

294) 

295) 

29 6) 

297) 

298) 

29 9 ) 

300 ) 

301 ) 

30 2 ) 
3C 3 ) 
3 r 4 ) 
39 5 ) 
306) 

7 ) 

? r 8) 
399 ) 

310) 

311) 

31 2 ) 

313) 

314) 
31 5 ) 

316) 

317) 
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APPENDIX P 


FORTRAN PROGRAM 


BEHIND A 


FOR CALCULATING THE AVERAGE WIND-TUNNEL 
INTERFERENCE OVER A TAIL 
SI OE-BY-SI DE ROTOR CONFIGURATION 


THIS PROGRAM WAS WRITTEN IN COC FORTRAN, VERSION 2.1, TO RJN ON COC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MCDIF ICATIONS MAY BE REQUIREO PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE ECU IV AL ENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS PROGRAM REQUIRES THE USE OF SUBROUTINE OLTAS WHICH IS GIVEN IN 
APPENCIX Q. 


INPUT WILL BE FOUND AT ADDRESS 1 (THREE CARDS PER CASE) IN FORMAT 900. NOTE 
THAT THE REFERENCE ORIGIN IS CHOSEN AT THE POINT MIDWAY BETWEEN THE TWO ROTORS. 
THE REQUIRED INPUT VARIABLES FOR THE ROTORS, GIVEN ON THE FIRST CARO, ARE 


LIR ROTOR DISK- LOAD -0 IS TR I8UTI0N INDICATOR, LIR=I FOR UNIFORM LOAD- 

ING, L IR=2 FOR TRIANGULAR LOAOING 

ZETA1 SEMIHEIGHT OF TUNNEL DIVIDED BY HEIGHT OF ORIGIN ABOVE FLOOR 

ETA1 DISTANCE FRCM ORIGIN TO RIGHT-HAND WALL DIVIDED BY TUNNEL 

SEMIWICTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMAR RATIO OF ROTOR DIAMETER TO TUNNEL WIDTH 

SR DISTANCE BETWEEN ROTOR CENTERS DIVIOED BY ROTOR DIAMETER 

ALPHAR ANGLE OF ATTACK OF TIP-PATH PLANE OF ROTORS, DEG 

THE required input variables for the wing, given on the second card, are 

LIW WING SPAN-LOAD-DISTRIBUTION INDICATOR, LIW=1 FOR UNIFORM LOADING, 

L I W=2 FOR ELLIPTIC LOADING 

LW DISTANCE OF WING APEX BEHIND ORIGIN, NONOI MENS I ONAL I ZE 0 WITH RE- 

SPECT TO ROTOR RADIUS 

HW DISTANCE OF WING APEX ABOVE ORIGIN, NON DI MENSI ON AL I ZED WITH RE- 

SPECT TO ROTOR RADIUS 

SIGMAW RATIO OF WING SPAN TO TUNNEL WIDTH 
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Appendix P - Continued 


LAMBDA WING SWEEP ANGLE. DEG 

ALPHAS ANGLE OF ATTACK OF BODY CARRYING WING AND TAIL, DEG 
THE REQUIRED INPLT VARIABLES FOR THE TAIL, ON THE THIRD CARD, ARE 
SIGMAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL TAIL LENGTH BEHIND ORIGIN AT ZERO ANGLE OF ATTACK, NONDI MENS ION- 

ALIZED WITH RESPECT TO ROTOR RADIUS 

TH TAIL HEIGHT ABCVE ORIGIN AT ZERO ANGLE OF ATTACK, NOND I MENSI ON- 

ALIZED WITH RESPECT TO ROTOR RADIUS 


THIS PROGRAM COMPUTES INDEPENDENTLY THE INTERFERENCE AT THE TAIL CAUSED 
BY THE PRESENCE OF THE THREE LIFTING ELEMENTS. IN SYMMETRICAL CASES, THE EFFECT 
OF THE TWO ROTORS IS IDENTICAL. CONSEQUENTLY, ONLY THE INTERFERENCE CAUSED BY 
CKE ROTOR IS CALCULATED. 


THIS PROGRAM REJECTS CASES OF ZERO SPAN. SINCE THE EQUATIONS ARE FORMED 
IN TERMS OF ROTOR RADIUS, SUCH CASES REPRESENT INPUT ERRORS. THIS PROGRAM ALSO 
DETERMINES AND REJECTS CASES IN WHICH TOTAL ROTOR SPAN EXCEEDS THE TUNNEL WIDTH. 


NOTE THAT THIS PROGRAM IS ALSO SUITABLE FOR TWI N- PROPELLER TILT-WING MOD- 
ELS, TILT-ROTOR MODELS, FAN-IN-WING MODELS, AND OTHER SIMILAR TYPES. JUDICIOUS 
CHOICE OF INPUT VARIABLES WILL SATISFY THE REQUIREMENTS OF THESE, AND MANY 
OTHER MODELS. 


PROGRAM WINOTUNI INPUT .OUTPUT ,TAPE5= I NPUT, TAPE6=0UTPUT ) 
COMMON ZETA,ET A.GAMMA.XOVERH .YOVERH, ZOVERH, DELTA (28 ) 
DIMENSION XDELTAI28 ), PS 1(20 > ,RLOAD(20) ,RUNIF< 20) ,RTRI A( 20) , 
1 XLE( 10) ,XL0AD(2C) ,C( 8) 

REAL LAMBDA, LW 

DATA (RUN IF ( I) , 1*1 , 2 0 ) /4*0. 296 1 ,8*0. 6255 , 8* 0. 8921 / 

DATA (RTRIA(I) »l=l»20)/4*0.4386» 8*0.7296,8*0.9262/ 

CAT A (C( 1), 1 = 1*81/20. ,30. ,90. ,50. ,60 . ,70. , 80. , 90 ./ 

PI=3. 14159265358979 
RAD=0. 01 74532925 199 
DO 803 Ll=l,28 
803 XDELTAILl ) *0 . 

PSI(1)=(PI/4.I 
PS I ( 2 ) = 3 .*PS 1(1 ) 

PSI (3) =5.*PSI( 1 ) 

PSI(4 )=7.*PSI( l> 

PSI( 5 ) = P$ I( 1 3 ) — ( PI/8. ) 

PSI(6)=PSl( 14>=3.*PSI( 5) 

PSI(7)=PSI(15)=5.*PSI(5) 

PSU8)=PSI( 16> = 7.*PSI (5) 

PSI(9) = PSI(17)=9.*PSI(5) 

PSI( 10)=PSI( 18)=11.*PSI (5) 

PSI(11)=PSI( 19)=13.*PSI(5> 

PSI(12)=PSI(20)=15.*PSI(5) 

XLE< 1 )= XL E(10)=0. 43579 
XLE ( 2)=XLE(9) =0.71422 
XLE ( 3 ) =XLE ( 8 ) =0.86603 


(P 

1) 

( p 

2) 

(P 

3) 

<p 

4) 

( p 

5) 

< p 

6) 

( p 

7) 

( p 

8) 

(P 

9) 

( p 

10) 

< p 

11) 

( p 

12) 

(P 

13) 

( p 

14) 

(P 

15) 

< p 

16) 

(p 

17) 

(P 

18) 

( p 

19) 

( p 

20 ) 

< p 

21) 

( p 

22) 

< p 

23) 

{P 

24) 

( p 

25 ) 

(P 

26) 

(p 

27) 
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Appendix P - Continued 



XLE ( 4 ) = XL E ( 7 ) =0.95394 

( P 

28) 


XLE ( 5 ) =XL E ( 6 ) =0.99499 

(P 

29) 

1 

READ (5,900 ) L I R , Z ET A 1 , £T A 1 , GAMMA , S IGM A R , SR , AL P HAR ,LIW,LW,HW, 

( P 

30) 


1 S IGM AW ♦ LAMBDA , AL PH AB , SI GM A T , TL , TH 

< P 

31 ) 


IF (EOF, 5 ) 999,47 

(P 

32) 

47 

A AL PR= AL PHAR 

( P 

33) 


aalpb=alpfab 

( P 

34) 


AL AM=L AMEC A 

( P 

35) 


L4MBDA=LAM8DA*R4D 

(P 

36 ) 


ALPHAR= ALPHAR*RAD 

( P 

37) 


ALPH AB = ALPH AB*RAD 

l P 

38 ) 


WRITE (6, nOJ 

(P 

39) 


IF (LIR.EQ.l ) GO TO 804 

( P 

40) 


I AL PHA= 10 HTR I ANGULAR 

( P 

41 ) 


00 808 M2 = l ,20 

(P 

42) 

808 

RLOAO(«2 l=RTRIA( M2 ) 

( P 

43) 


GO TO 806 

( P 

44) 

804 

I ALPHA = 10HUNIF ORM 

(P 

45 ) 


00 809 M2= 1 , 20 

( P 

46 ) 

8C9 

RLO AO (M2) =RUN I F ( t*Z ) 

(P 

47) 

806 

IF (LIW.EG.l ) GO TO 852 

( P 

48) 


IBETA=8FELLIPTIC 

( P 

49) 


CO 851 M 3= 1 , 10 

(P 

50) 

851 

XL0AD(M3 )=XLE(M3) 

(P 

51) 


GO TO 850 

( P 

52) 

852 

I BET A=8 HUN I FORM 

(P 

53) 


CO 853 M 3= 1 ,10 

( P 

54) 

853 

XLO A0< M 3 ) = 1 .0 

( P 

55 > 

850 

WlDTH=SlGMAR*( 1.0+SR) 

( P 

56) 


IF (WIDTH. LT.1..ANO.SIGMAR.NE.O.) GO TO 855 

( P 

57) 


WRITE (6,701 ) 

( P 

58) 


WRITE (6,706) I ALPHA , I BET A, S IGMAR , Z E TA1 ,LW, AALPR ,S I GMAW , ET A1 , 

( P 

59) 


1 HW, AALPB,S IGM AT, TL ,TH,SR , GAMMA, ALAM 

(P 

60) 


WRITE (6,210) 

(P 

61 ) 


WRITE (6,211 ) 

( P 

62) 


WRITE (6,212 ) 

( P 

63 ) 


WRITE (6,213 ) 

(P 

64) 


WRITE (6,214) 

(P 

65) 


WRITE (6,215 ) 

( P 

66 ) 


WRITE (6,216) 

( P 

67) 


WRITE (6,217) 

(P 

68) 


WRITE (6,218) 

( P 

69) 


IF (SIGMAR.EO.O. ) WRITE (6,707) 

( P 

70 ) 


IF (WIDTH. GE.l .) WRITE (6, 708) 

(P 

71) 


GO TO 1 

( P 

72 ) 

855 

CO 42 I EL EM =1 , 3 

( P 

73) 


IF (ETAl.EQ.I.) GO TO (600,42,600), I EL EM 

(P 

74) 

600 

WRITE (6,701 ) 

( P 

75) 


IF (ETA1.EQ.1.) GO TO (604,42,603), I ELEM 

(P 

76) 


GO TO ( 601,602,603), I EL EM 

(P 

77) 

601 

WRITE ( 6,702 ) 

( P 

78) 


GO TO 610 

( P 

79) 

602 

WRITE (6,703) 

( P 

80 ) 


GO TO 610 

( P 

81) 

6C3 

WRITE (6,704) 

(P 

82) 


GO TO 610 

( P 

83) 

604 

WRITE (6,705) 

( P 

84) 

610 

WRITE (6,706 ) I ALPHA , I BET A, S IGMAR , ZE T A 1 , LW , AALPR , SI GMAW, ETA 1 , 

( P 

85) 


1 HW, 44LP3, SIGMATtTL ,TH, SR , GAMMA, ALAM 

( P 

86) 


WRITE (6,210) 

( P 

87) 


WRITE (6,211) 

(P 

88) 
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Appendix P - Continued 


WRITE ( 6,212 ) 

WRITE 16,213) 

WRITE 16,214) 

WRITE ( 6,215) 

WRI TE l 6,216) 

WRITE (6,217) 

WRITE (6,218) 

CO 41 K. = 1 » 8 
N 7 = 2 C 
M 7=4 

GC TO (611,611,612 ), IELEM 
EFFECT CF ROTOR 

611 SUML=O.Cl 2 5 
GO TO 812 

EFFECT OF WING 

612 SUML =0 .0315 26 

IF (SIGMAW.EQ.O* .OR.SIGMAT. EG. 9. I GO TO 615 
IF (L IW .60.1 ) SUML-0 .025 
N7= 10 

IF ( ETA1 .NE.l . ) GO TO 8 12 
SUML -0*0126104 
IF CLIW.EO.il SURL =0 *05C 
N7= 5 

GO TO 812 

615 IF (SIGMAW.EO.O. *ANO*SIGMAT *NE*0. ) GO TO 613 
IF CS1GMAW.NE.0. .AND .SI GMAT . EO.O. ) GO TO 616 
SUML= 1 • G 

X LO AD ( 1 1=1.0 
R7 = N7 « 1 
GO TO 812 

613 S UML=0 * 2 5 
P 7 = 4 

N7= 1 

X LO AC ( 11 = 1.0 

IF ( ETA 1 *NE . 1 * I GO TO 812 
SuML=0 .50 
M 7 = 2 

GO TO 812 

616 SUML = C • 1261 C4 

I F (L IW .EO.l ) SUML-O. IOC 
N 7= 1 0 
M7= 1 

IF CETA1.NE.1.I GO TO 812 
SUML=C. 25208 

IF CLIW.EQ.l) SUML=0.20C 
N7=5 

812 DO 801 M 1 = 1 ,M7 
DO 802 N l =1 , N7 
XST AR=( 11 .-2 .♦FLOAT { Ml 11/10. 

Y ST AR= ( 11.-2. ♦FLOAT (N1 > ) / 1 0 • 

GO TO ( 622,623,621 ), I EL E ^ 

EFFECT OF RIGHT ROTOR 

6 22 ETA=ETA1-RLQA0(N1 )*SIGMAR*SIN( PS I (Nl) )-SR*S IGWAR 

YOVERH=-SIG N 1AR*GAMMA*(0. 25* ( 2. 0*FL0 A T ( M 1 ) - 5 . 0 1 * ( SI GMA T / S I G M AR ) 
1 *RLOAO( Mi )*S IN( PSI (M) ) + SR ) 


( P 89 ) 
(P 90 
(P 91) 
( P 92 ) 
(P 93) 
( P 94) 

( P 95) 
(P 96) 
(P 97) 
( P 98) 
( P 99 ) 
(P K' ? ) 
( P 101 ) 
(P 102) 
(P ir-3) 
(P 1 n 4 1 
(P 1*5) 
(P 106) 
(P 1071 
(P 1 °8 ) 
(P 1*9) 
(P 11?) 
(P 111) 
(P 112) 
(P 113) 
(P 114) 
( P 115) 
(P 116) 
(P 117) 
(P 118) 
(P 119) 
(P 120) 
(P 121) 
(P 1 22) 
(P 123) 
(P 124) 
(P 125) 
(P 126) 
(P 127) 
( P 128) 
(P 120) 
(P 130) 
(P 131 ) 
(P 132) 
(P 133) 
(P 134) 
(P 135) 
(P 136) 
(P 137) 
(P 138) 
(P 139) 
(P 140) 
(P 141) 
(P 142) 
(P 143) 
(P 144) 
(P 145) 
(P 146) 
(P 147) 
(P 148) 
(P 149) 
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Appendix P - Continued 


GO TO 625 

EFFECT OF LEFT ROTOR 

6 23 ETA=ETA1-RL0A0IN1 > *S IGM4R*SINI PS 1 1 N l)> +SR*S I GM4R 

Y0VERH=-SIGMAR*GAMMA*I 0.25*1 2.0*FL3ATIM1)-5.0>*I SIGMAT/SIGM AR> 

1 ♦RL0ADIN1>*SIN(PSIIN1>)-SR» 

GO TO 625 

EFFECT OF WING 

62i YST4R = I 11.0-2. 0*FL0AT ( N1 D/10.0 
ETA=ETA1-YSTAR*S IGMAW 

ZETA=ZETA1/I 1.0-SIGMAR*GAMMA*ZETA1*I A8SCYST AR ) *1 SIGMAW/SIGM4R > 

1 *T ANI L AMBDA ) *S INI ALPHA8 > *LW*SIN l AL PHABI -HW*CCS< ALPHAS) ) ) 

X0VERH= SIGM AR*GAMM A* ( ( TL-LW > *C0S I ALPHAS )*{ TH-HW, *SINI ALPHAB ) 

1 -ABSIYSTAR >*l S IGMAW/ SIGM AR »*T AN (LAMBDA ) *COS I ALPHAB) ) 

YCVERH=GAMM A*I 0.25*1 5.0-2. 3*FL0ATIM1» ) *SIGMAT-YSTAR*SIGMAW) 
Z0VERH=-SIGMAR*G4MM4*I ITL-LW,*S IN, A LPHAB » -I TH-HW )*COS ( ALPHAB) 

1 -ABSIYSTAR ,*IS IGMAW/ SIGN AR »*T AN I LAMBDA) *S IN I ALPHAB) ) 

GO TO 6 30 

62 5 ZETA=ZETA1/,1.0-RL0A0,N1 » *S I GM AR*GAMM4*ZET A1 *C3S IP S 1 1 N1 ))* 

1 S INI ALPHAR ) ) 

XOVERH* SIGM AR*G AMM 4* ITL* COS I AL PHAB 1 *TH*SINt ALPHA B)- 
1 RLO AO INI) *C0S I ALPHAR ,*COS I PS UNI) )) 

Z0VERH=SIGM4R*G4MM4*,TH*C0S< ALPHAS, -TL*S INI ALPHAB) +RL 040 INI ) 

1 *SINI ALPHAR) *COS IPSIIN1 )) ) 

XLOAOI NI )=1 .0 
630 CALL OLTAS I C I K I » 

*+********** 4 ******** see appendix o for subroutine dltas ********** 

00 805 Ll = i * 28 

805 XOELTAILl )*XDELTAILI ) +DELTAI Lll *XLO ADI N 1 > 

802 CONTINUE 
801 CONTINUE 

DO 807 L3=i *28 

807 CELT A I L3 )=XDELTAIL3) *SUML 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
DO 81 A L4=1.28 
814 XOELTAI L4 )=0. 

41 CONTINUE 

42 CONTINUE 
GO TO 1 

100 FORMAT IIH/////////59X*START NEW CASE*) 

149 FORMAT 1 1 X *CHI = * F7.3/) 

150 FORMAT I 3X5HI W ,L I7IF17.4 ) ) 

151 FORMAT (3X5HIU.L J7IF17.4)) 

152 FORMAT I 3X5H ( W , 0) 71 F17 .4 ) ) 

153 FORMAT I3X5HIU.0I7IF17.4 )// ) 

210 FORMAT UX131I1H-I) 

211 FORMAT I 1X1HI11X1H 131X6 1HC3RRECTI0N FACTORS FOR CORRECTING FROM A 
1WIND TUNNEL WHICH IS25X1HI) 

212 FORMAT I 1X1HI 1 IX 1H 1 11 71 1 H- ) 1HI ) 

213 FORMAT I lXIHIllXIHI 16X1HI5X6HCL0SED5X1HI 16XIHI 2X12HCLOSEO FL00R2X1 
1H I6X4H0PEN6X1HI 1 6X1HI5X6HCL0SED4X1HI ) 

214 FORMAT 11X1HI3X5HOELTA3X1HI 5X6HCL0SE05X1HI 4X9H0N 80TT0M3X1H I 6X4 HOP 
1EN6X1HI 6X4HCNL Y6X1HI 5X5HFL00R6 X 1H 1 5X6HCL 0SE05X1H I3X9H0N B0TT0M3 XI H 
21 ) 


( 6 , 

1491 

C ( K 1 





16, 

150 ) 

(DELTA* 

n. 

r 1 = 1 i 

,25 

,4) 

(6, 

151 > 

(DELTA ( 

i), 

>1 = 2, 

> 26 

,4) 

(6, 

152 > 

(DELTA! 

i ) i 

>1=3, 

>27 

,4) 

(6, 

153) 

(DELTA! 

i ) i 

> 1=4, 

>28 

,4) 


(P 

150) 

f P 

1511 

(P 

152) 

(P 

153) 

(P 

154) 

(P 

155) 

(P 

156) 

(P 

157) 

(P 

158) 

(P 

159) 

( P 

160) 

(P 

161 ) 

(P 

162) 

(P 

163) 

(P 

164) 

(P 

165) 

(P 

166) 

<P 

167) 

(P 

168) 

(P 

169) 

(P 

170) 

(P 

171) 

(P 

172) 

(P 

173) 

<P 

174) 

(P 

175) 

(P 

176) 

(P 

177) 

(P 

178) 

t******* 

(P 

179) 

(P 

180) 

(P 

181 ) 

( P 

182) 

(P 

183) 

(P 

184) 

CP 

185) 

(P 

186) 

( P 

187) 

(P 

188) 

(P 

189) 

(P 

190) 

( P 

191) 

<P 

192 ) 

CP 

193) 

(p 

194) 

(P 

195) 

<p 

196 ) 

(P 

197) 

(p 

198) 

(p 

199) 

(p 

200) 

(p 

201 ) 

(P 

202) 

(p 

203) 

(P 

204) 

(P 

2 05) 

(P 

206) 

(p 

207) 

CP 

208) 

(P 

209) 
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Appendix P - Concluded 


215 FORMAT ( 1 XI H II IX 1H I 16X1H I6X4H0NLY6X1H I 16X1 8HI (GROUND EFFECT) I6X4H 
10NLY6X1HI16X1H I6X4H0NL Y5 X1H I ) 

216 FORMAT ( IX IF 11 1 XIH 1 84( 1H- >1 H 13 2 ( 1H- ) 1HI ) 

217 FORMAT ( 1 XI H 1 1 IX 1H I 36X1 1 HTO FREE AI R37X 1 HI 8X1 6HT C GROUND EFFECT 8X1 
IF I ) 

218 FORMAT 
701 FORMAT 


I 1X13K 1 H- ) / ) 
(1H1//28X*AVERAGE 


INTERFERENCE OVER A TAIL BEHIND SIDE-BY-S 


1IDE AND/OR TILT ROTOR MODELS*/) 

702 FORMAT <55X*EFFECT OF RIGHT ROTOR*//) 

703 FORMAT !56X*EFFECT OF LEFT ROTOR*//) 

704 FORMAT <59X*EFFECT OF WING*//) 

705 FORMAT (55X*EFFECT OF EITHER ROTOR*//) 

706 FORMAT (34X.A10* ROTOR LOAOI NG*19X , A8* WING LOADING*// 

119 X* SIGMA (ROTOR) =*F6.3, 10X*ZETA =*F6. 3 , 10X*LW/R =*F6.3»10X 


2*ALPHA( ROTOR) =*F7.3//19X*SIGMA( WING) =*F6. 3. 10X*ETA =*F6.3, 

31 OX *HW/R =*F6 . 3 » 10X* ALPHA! BOOY ) =*F7.3//1 9X*SIGMA( TAIL ) =*F6.3, 

410X*TL/R =*F6.3,10X*TH/P =*F6. 3, 10X*SR/R =*F7.3//39X 

5*GAMM A =*F6.3,29X*LAMBDA =*F7.3//) 

707 FORMAT (40X*SI GMA (ROTOR ) EOUALS ZERO, THIS PROGRAM IS NOT SUITABLE 
1 FOR USE WITH SUCH CASES.*//) 

708 FORMAT (40X*R0T0R SYSTEM IS TOO WIDE FOR WINO TUNNEL*//) 

900 FORMAT ( 1 1 , F9. 3 , 5F1 0 .3 /1 1 ,F9. 3 ,4F 10 . 3/ 3F 10. 3) 

999 STOP 
END 


CP 

210) 

(P 

2 1 i ) 

CP 

212) 

CP 

213) 

t p 

214) 

(P 

215) 

(P 

216 ) 

CP 

217) 

(P 

218) 

CP 

219) 

CP 

220) 

CP 

221) 

(P 

222 ) 

CP 

223) 

CP 

224) 

CP 

225) 

CP 

226) 

CP 

227) 

CP 

228) 

(P 

229) 

CP 

230) 

CP 

231 ) 

(p 

232) 

CP 

233) 
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APPENDIX Q 


SUBROUTINE DLTAS 


THIS SUBROUTINE WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS SUBROUTINE 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS SUBROUTINE IS REQUIRED BY THE PROGRAMS OF ALL PRECEEDING APPENDICES. 


SUBROUTINE DLTAS IANGL) (Q 1) 

COMMON ZETA, ETA, GAMMA, XOVERH , YOVERH .ZOVERH , OELTA ( 2 8 ) (Q 2) 

DIMENSION V(3,9) ,ADEL( 28 ) IQ 3) 

SC*S INI 4NGL*0 .01 74532925199 ) IQ 4) 

CC=C0SUNGL*0. 0174532925199) IQ 5) 

Z6=ZETA*Z0VERH+1. IQ 6) 

Z8=-Z 6 (Q 7) 

Z7=Z8-1. IQ 8) 

DO 8 J 1=1 , 28 (Q 9) 

8 DELTA! J1 )=C. I Q 10 ) 

CO 10 M = 1 » 7 IQ ID 

DO 10 N=1 *7 (Q 12) 

IF IN.EQ.4.AND.M.E0.4) GO TO 10 IQ 13) 

DO 11 J 1=1 , 3 (Q 14) 

DO 11 J2= 1 » 9 (Q 15) 

11 VI J1 , J2 )=0. IQ 16) 

DO 12 J 1=1 , 28 IQ 17) 

12 A DEL I J1 )=0. IQ 18) 

AM=M-4 (Q 19) 

AN=N-4 IQ 20) 

X = ZET A* XOVERH IQ 21) 

Y=ZETA*(Y0VERH-2.*AM*GAMM4fGAMMA*I l.-ETA»*( l.-I-l. )**M) ) (Q 22) 

Z=ZET A* ( ZOVERH-4 .*AN I IQ 23) 

4=SQRT< X*X+Y*Y+Z*Z» 10 24) 

e=A+Z*CC-X*SC IQ 25) 

VII, 1 )=! I X*X*Y*Y )/!B*A*A*A) )-IlZ*A*CC)/I B*4 I )**2 IQ 26) 

V(2*l ) =-I X*Z )/ I B*A*A*A )- IZ+A*CC ) *( X-A*SC )/ I 8*B* A*A ) IQ 27) 

VI3,1 )=I I Y*Y*Z*Z)/18*4*4*4I )-{ ! X-4*SC )/IB*4 ) >**2 IQ 28) 

Z=-Z-2. IQ 29) 

A=SQRT( X*X+Y*Y+Z*Z ) (Q 30) 

B=A+Z*CC-X*SC IQ 31) 

V(1,3) = I IX*X+Y*Y )/CB*A*A*A) )-! I Z* A*CC ) / I B*A ) )**2 (Q 32) 

VI2 ,3 ) =-( X*Z) / I B* 4* 4*41-1 Z+A*CC ) *( X- A*SC ) / I B*B* A*A ) (Q 33) 

V 13 ,3 )= ! I Y*Y+Z*Z >/ I8*A*A*A) )-I IX- A* SC )/ l B*A ) ) **2 IQ 34) 

IF (ANGL. £0.90.0) GO TO 13 IQ 35) 

X=X-ISC/CC) IQ 36) 

Z=-Z-l. IQ 37) 

A=SQRT { X*X*Y*Y *1 *Z ) IQ 38) 

B=A»Z*CC-X*SC 10 39) 

V(1,2)=HX*X+Y*Y)/IB*A*A*A))-(IZ+A*CC)/(B*A))**2 (Q 401 

V(2,2)=-IX*Z)/IB*A*A*A)-IZ+A*CC)*<X-A*SC)/(B*B*A*A> IQ 41) 
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Appendix Q - Continued 


V! 3,2)=! t Y*Y+Z*Z)/(B*4*A*4> ) - ( ! X- A*SC ) / ( B*A I > **2 
B= A— X 

V!l.5)=( !X*X+Y*Y)/!B*A*A*A) ) -< Z/ < B* A )) **2 
V!2,5)=Z/! A*A*A) 

VI3,5)=X/!A*A*A) 

Z=-Z 

B = A+Z *CC-X*SC 

V!i»4)=lCX4X*-Y*Y)/l B*A*A*AI )-((ZH*C:i/l8*MI**2 
V!2,4)=-!X*Z)/ !B*A*A*A)-!Z+A*CC )*!X~A*SC )/(B*B*A*A) 

V! 3,4 ) = ! ( Y*Y+Z*Z)/!B*4*4*4> )-! f X- A*SC ) / ! B*4) ) **2 

13 A DEL ( 1 1 =V ( 1 , 1 ) -V ( 1 . 2 ) -V ( 1 , 3 > +V ( 1 , 4 > 

ADEL(2)=V!2,i)-V(2,2)+V! 2,3 )-V(2 ,4) 

ADEL!3)=V!2, ll-V(2,21-V( 2,3) +V< 2, A) «-2.*V( 2,5) 

ACEH 4)=V!3»1)-VC3»2)+V(3»3)-V(3»4)+2.*V!3,5) 

4DEL!5)=< !-l.) **!M+N) )*4DELt 1 ) 

ADEL16>=! !-l . )**!M+N))*AOEL! 21 
ACEL ! 7) =! (-1 *) ♦*! M+N ) ) *ADEL (3) 

ADEL !8) = ! (-1.>**(M+NI)*ADEL(4) 

ADEL ! 9 ) = (I -1. )**«)*( V( 1 , 1 1 - Vt 1,2)+Vt 1, 3 ) -V ( 1 , 4 ) +2 , * V ( 1 ,5) ) 
ACEL(10)=C!-1, J<MMI*!V12,1)-V(2,2)-V! 2, 3) *V! 2, 4 ) +2 .*V!2 ,5 ) ) 
ADEL (11 )= ! 1 -1 • l**M)*(V(2,i)-V(2,2)+V( 2,3)-V(2,4) ) 

ADELC12 )=l!-l. J**M)*!V!3,i)-V< 3, 21- V! 3, 3 ) +V( 3 ,4) ) 

DO 14 J 1 = 1 , 12 

14 DELTAC J1)=DELTA( J1)+ADEL( Jl ) 

10 CONTINUE 

DO 15 Jl«l,8 

15 DELTA! J 1 +20 )=DELT A ! J1 ) 

X=ZET4*XCVERH 

Y = Z ET A*YOVERH 
Z = Z 7 

A=SQRT ( X*X+Y*Y+Z*Z ) 

B=A+Z*CC-X*SC 

V!1,7)=!!X*X+Y*Y)/!B*4*A*4) )-! <Z+A*C: )/! B*A))**2 
V!2»7)=-(X*Z)/!B*A*A*A)-( Z> A*CC I *( X-A*SC )/(B*B*A*A) 

V<3.7>=< ! Y*Y+Z*Z )/!B*A*4*4) ) - ! ( X- A* SC ) / ! 8*4 )) **2 
IF CANGL.EQ.90.0 ) GO TO 16 
X=X-! SC/CC ) 

Z=Z6 

A=SQRTC X*X*Y*Y*Z*Z ) 

B=A+Z*CC-X*SC 

V C l » 6 ) = ( !X*X*Y*Y )/ ( B*A*A*A) )- ( ( Z+A*CC ) / ( B*A ) I **2 

VI 2,6) = - 1 X*Z) /<B*4*A*4)-(Z+4*CC)*! X-4*SC )/(B*B*4*4) 

V ( 3,6 )= ( ( Y*Y*Z*Z )/!B*A*A*A) ) - ! ( X- A*SC )/ ( B*A ) ) **2 
B= A-X 

V!l»9)=!!X*X*Y*Y)/( B*A*A*A) ) - ( 1/ ( B*A ) ) **2 
V!2,9)=Z/!A*A*A) 

V(3,9J = X/ (A* A* A J 
Z = Z8 

B=A+Z*CC-X*SC 

V ( 1 , 8 )= ( ! X*X+Y*Y)/!B*A*A*A) )-! ( Z*A*CC)/C B*A) )**2 
V!2,8)=-!X*Z )/IB*A*A*4)-!Z+A*CC)*<X-A*SC )/(B*B*4*4» 

V (3 , 8 ) = ( I Y*Y >Z *Z )/ C 8 *AAA *A ) )-( ! X- A*SC) / ( B*A) )**2 

16 DELTA! 13 )=-V(l ,6)-V( 1,7 I ♦VC 1,8) 

DELTA! 14 ) = -V I 2 , 6 )*V< 2 ,7 ) -V! 2 ,8 ) 

CELT A! 15 )=-V 12 , 6 )-V ! 2 , 7 ) +V! 2 ,8 ) +2 . *V! 2 , 9) 

DELTA! 16)=-VI3,6)+V!3,7)-V(3,8)*2.*V(3,9) 

DELTA! 17 ) = -V (1,6)+V!1»7)-V( 1 , 8 ) +2 . *V< 1 , 9 ) 

DELTA! 18 ) =DE LT A < 1 5 ) 

DELTA! 19 ) = DEL T A! 14) 

DEL T A! 20 )=-V ! 3 , 6 )-V( 3, 7) +V! 3,8 1 
DO 17 J 1= 1 , 4 
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Appendix Q - Concluded 


17 CELT A! J1 )*DELTA! Jl ) + DELT A! J 1 + 12) (Q 1 03 » 

DO 18 Jl-5« 12 (Q 10A) 

18 DELTA I J1 )=DELTA!Jl)+DELTA!Jl+8) !Q 105) 

AMT=-2.*GAMMA*ZETA*ZETA/3. 14159265358979 (0 106) 

CO 19 J 1*1 ,28 (0 107) 

19 DELTA! J1 )=AMT*DELT A< Jl) !0 108) 

RETURN (Q 109) 

END (Q 110) 
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